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ABSTRACT 


' The same ratio correlation may be generated by many different combinations of relationships between 
' absolute measures, but a single set of absolute-measure statistics leads to one, and only one, correlation 
~ between any particular set of ratios formed from these absolute measures. The passage from ratio correlation 
© toinference about relations between absolute measures is ambiguous at best and often misleading. 

Algebraic statements exhibiting ratio correlation as a function of absolute-measure statistics are offered 
for types of ratios commonly used in petrography. These statements are all derived from Pearson’s general 
' formula for index correlation. They yield good approximations only if the fraction s/X¥ for each absolute 
> measure is suitably small. 

Several practical examples drawn from petrographic literature are described. In most of these cases 
» the ratios seem to have been used either to order the data or in the hope that they would throw some light on 
' telationships between absolute measures. The results are shown to be on the whole indecisive and ambigu- 
* ous and in a few cases decidedly misleading. 

The formation of ratios should be confined to those problems in which hypotheses being tested deal with 
' tatios. Absolute measures are always preferable when large numbers of observations must be recorded with- 
out benefit of satisfactory hypothesis. Ratios can always be drawn from tables of absolute measures; fre- 
' quently, absolute measures cannot be reclaimed from tables of ratios. 


INTRODUCTION a duty or a pleasure, our task is to 
search for order and even sometimes to 
invent it when we cannot find it; we 
cannot help preferring an orderly, sys- 
tematic result to a disorderly, unsystem- 
atic one. 

Unfortunately, the search for order 
easily gets out of hand, and this is par- 
ticularly the case when it is necessary to 
record and assimilate large numbers of 
observations concerning which our hy- 
potheses are either inadequate or inde- 
cisive. We do not understand the signif- 
icance of the observations, but it is still 
incumbent upon us to record all and re- 
member a good many of them. It is quite 
as impossible to record a chaos as to re- 
‘Manuscript received August 4, 1948. member one, and so we are obliged to 


In petrography, as in any natural sci- 
ence, one of the principal drives of the 
observer is to order or ‘‘make sense”’ of 
"his data. We all share the suspicion that 
| the truth, if we could discover it, would 
"be systematic; most of us would prob- 
'ably go even further and confess that 
when we have discovered a little order in 
| our data we feel that we have come one 
| Step closer to the truth. The procedure is 
"logically untenable, for history shows 
‘that false systems are possible; but to 
"question its psychological validity would 

be to question the whole basis of de- 
| Scrip tive science. Whether we regard it as 


239 





240 


impose system where we have not as yet 
been able to discover it. 

Where the observations consist of 
measurements of two or more variables, 
some of our commoner notions of order 
and system are generalized in the con- 
cept of correlation. If the measurements 
are uncorrelated, many ratios formed 
from them will generally exhibit correla- 
tion. Purely for the sake of efficient de- 
scription, the ratios seem more useful 
than the absolute values do. Their utility 
in one situation recommends their ap- 
plication in others, and we tend to over- 
look the fact that ratio correlations rare- 
ly permit unambiguous inferences about 
absolute values. 

Ratio correlation is an old subject in 
statistics—so old, in fact, that it is not 
even mentioned in some of the newer 
texts.2 The purpose of this paper is to 
restate, in petrographic context, some of 
the more important properties of ratios 
in correlation analysis. No revolutionary 
advance is likely to result from the at- 
tempt, but I hope it will intrigue other 
petrographers, as it has intrigued me, to 
discover that what we have learned 
piecemeal and with much pain has been 
available for the asking in the journal 
and text literature of a sister-science for 
over half a century. 


DEFINITIONS AND TERMS 


Definitions 1-3 are used to reach 
definitions 4 and 5; only 4 and 5 are re- 
quired in the discussion which follows. 


1. & is the mean of X,,, X;., X4;,... 
X;, Measurements, or 


i= LL > be 
i=l 


2It is treated briefly by Snedecor, Statistical 
Methods (4th ed., 1946), pp. 162-64, and at greater 
length by Yule and Kendall, An Introduction to the 
Theory of Statistics (13th ed., 1946), pp. 299-301. 


F. CHAYES 


2. %,; is the deviation of the ith X, from 
the mean of the X,’s, or 


ei = (Xai — Fe) - 


. 5, is the standard deviation, or square 
root of the variance, of the X,’s, or 


( ; , 
Ss, = ( —— sy. 
n— 1 Some i 

i=1 


. C, the coefficient of variation, is the 
ratio of standard deviation to mean, 
or 


. r, the product moment coefficient of 
correlation, is the ratio of the covari- 
ance of two variables to the geometric 
mean of their variances, or 


1 n 
——7 2 (41; %2;) 
i=l a p 
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The correlation coefficient is a re- 
markably useful index of the relation 
between two variables. It will always be 
true that —1 <r <1; if r = |1| the 
variables stand in perfect linear rela- 
tion to each other, while if r = o there is 
no linear correlation. If the variability 
be expressed as variance, r? gives the pro- 
portion of it attributable to linear regres- 
sion of either variable on the other; 
(1 — fr’) is the amount of “‘residual”’ vari- 
ation which cannot be described in this 
fashion. If r >o0, the quantities in- 
volved vary directly; if r < o their varia- 
tion is inverse. We use r in this note pure- 
ly as a sample description, so that in- 
formation about the distribution of the 
parent deviates is unnecessary. 


PEARSON’S GENERAL FORMULA FOR 
RATIO CORRELATION 


More than half a century ago, K. 
Pearson (1896-1897) showed that, if 
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ratios were formed from absolute measures, such that 


ie 
- 


the correlation between the ratios was implicit in the variances, covariances, 
and means of the absolute measures. His derivation assumed that the absolute 
measures were normally distributed and that in each case the coefficient of varia- 
tion was small enough that its third and higher powers could be ignored. Reed 
(1921) later arrived at the same result without specifying the parent-distribu- 
tion of the absolutes. For our purpose it is important to note that the remaining 
restriction, governing the size of the C values, is required only to reach a con- 
venient algebraic formulation; correlation between the ratios is a single-valued 
function of the variances and covariances of the absolute measures, whether 
or not the relationship may be conveniently stated. Pearson’s general formula, 
which is rarely given in elementary texts, is 

— Cg — Puig — PoC gt tay (2) 
(CRA C2 = 2 rio) 7? (C2 +4. C2— 2 reg) 


and (1) 





r 


It contains quantities of only two types, and these are symmetrically distributed. 
The denominator will always be positive, so the sign of r,, depends on the sign of 
the numerator. The quantity r,, will be zero only if 1; = tig = 123 = f24 = 0;' if 
one of these correlations is not zero, r,, will also differ from zero. But it may 
differ in sign from the nonzero r on the right, and it will almost certainly differ 
in size. A single combination of r and C values can lead to one and only one value 
for r,., but the same 7,, might result from many different combinations of r and C 
values; this is the foundation for the assertion that correlations between ratios 
do not permit unambiguous inferences about the variables from which the ratios 
are formed. 


SPECIAL CASES OF RATIO CORRELATION 


By manipulating the absolute measures it is possible to develop from equation (1) 
simpler formulas for ratio correlations covering a wide variety of practical problems. 
In this section the more important relations are listed, and in the succeeding section 
afew practical applications are described. 


CORRELATION BETWEEN TWO RATIOS WITH COMMON DENOMINATOR 


If, in equation (1), X, = X,, so that 
Y= . and Z= : : 


then C, = Cy, 24 = 1, tie = Try, 23 = 154, and equation (2) reduces to 


ad rigCiC3 — rr C2 — rxC XC 3 +C3 
MO (C24C2 — 2 nC iC2) 7? (C2 + C2 — 2 reg C2) 7?" 





r 


3Or if fhe terms in the numerator cancel completely—a very unlikely circumstance. 
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This type of ratio correlation is com- 
mon in many descriptive sciences, for it 
is often convenient to use one variable as 
a sort of scaling unit. If the absolute 
measures are uncorrelated (e.g., 7. = 
11; = 2; = ©), equation (4) reduces to 


C 


i : 
yz (C?+C?2) 1/2 (C?+C?) 1/2? (5) 





which is Pearson’s definition of “spuri- 
ous correlation’ and will always be 
positive. Even in this greatly simplified 
form, many combinations of C values 
might yield the same r,,. For the special 
case in which C, = C, = C,, ry, = 0.5; if 
C, is notably larger than C, and C,, 7,, 
may be much larger than o.5. If C, = 
C,; = C,/2, for instance, r,, = 0.67. If 
C, = C,/2 = C./3, fy, = 0.79. 


CORRELATION OF A RATIO WITH A 
THIRD VARIABLE 


In equation (1) let X, = 1, so that 


dice s and Z=X;. (6) 


Then equation (2) reduces to 


Ni3C 1 — Togo 


=" CFC —-Intey® §) 





If neither term of the ratio is corre- 
lated with Z, r,, = o. If there is positive 
correlation between X, and X,, the size 
of r,, is increased, but its sign is unaf- 
fected. If r,, = 7.3, the sign of r,., and to 
a large extent its size as well, will de- 
pend on the sizes of C, and C,. If the 
terms of the ratio are completely uncor- 
related (7,, = 0), ry, will nevertheless 
differ from zero, provided only that 
either r,, or r., does so. 








CORRELATION OF A RATIO WITH ITS 
OWN DENOMINATOR 


In equation (1) let X, = X, and X, = 
1, so that 


be 


Y =—! and 


X Z=X2. (8) 





to 


Then equation (2) reduces to 


ae —C»2 
(C24C? = 2 rCiC2) 





(9) 


Puget 


which will be positive only if r..C; > C, 
If r,, = 0, equation (g) reduces to 


—C2 
ty: = ($C) (10) 
so that the “spurious” correlation is here 
negative. If, further, C, = C., 1, = 
—o.71, and if C, > C, it will be even 
larger; if C, is only twice as large as C,, 
for instance, 7), = —o.89. 


CORRELATION OF A RATIO WITH ITS 
OWN NUMERATOR 


In equation (1) let X, = X, and X,= 
1, so that 


a , , 
Y¥=> and z2=X,. (i 

Then equation (2) reduces to 
. Cir rile _ (12) 





(CRF CE= 2 riz) 


If r,. = 0, equation (12) further reduces 
to 

Vy: = waded -_ jo { 13) 

(C3+-C2) 

the “spurious” correlation again being 
positive, as in Pearson’s type case. If, in 
addition, C, = C., m2. = 0.71; if C, > G, 
it will be larger. 
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CORRELATION BETWEEN TWO RATIOS WITH THE SAME NUMERATOR 


This case is complementary to the original definition of spurious correlation (eq. 5) 
and is reached from equation (2) in similar fashion. If, in equation (1), X¥, = X;, so 
that 


» a an 
J a and 2= 7 (14) 








equation (2) reduces to 


Ci — 7ygiCg— TC iCot+ ral, 





= , 15 
™ (CE +CZ — 2 nC iC2) (CE ACE — 2 nix) ane 

If ro = 24 = 134 = 0, equation (15) reduces to 
Siete (16) 


(CFC) RCE) 


As in equation (5), the spurious correlation here is positive. If, in equation (16), 
(,= C. = Cy, ry. = 0.5, exactly as in equation (5). 

For the sake of simplicity the order of subscripts set up in equation (1) has been 
followed throughout. Notice, however, that equation (16) is actually identical 
with equation (5). In each case the numerator is the square of the coef- 
ficient of variation of the common term, and the denominator is the product of 
the square roots of the sums of squares of each unique term C with the common 
term C. A similar term-by-term identity holds for equations (15) and (4). Thus the 
ratio correlation will be identical whether the common term is used as numerator 
or as denominator, provided only that it occupies the same position in Y as in Z. 


CORRELATION BETWEEN RATIOS WHEN THE DENOMINATOR OF ONE IS THE 
SAME AS THE NUMERATOR OF THE OTHER 
If, in equation (1), X, = X;, so that 
— * Xe ‘ 
Y= >- and Z=3, (17) 
Xe A4 
equation (2) reduces to 
12 
MWC C2 — ryCiCg+ ragCXy —Ce 





‘yz T7"4 > sy /> ” 17% .7* \1/ ’ 1 8 

a (C3 ACF — 2 ri iC2) /? (CZ +C? — 2 regCxCy) rae 
and, if 7,, = 14 = 74 = 0, equation (18) ures are uncorrelated, the ratio correla- 
reduces to tion in this case will be equal in size and 
-C? opposite in sign to that between ratios 





“ (CFC) FCC) (19) with common denominator or common 

: ‘ numerator. In the same way, equation 

If, in addition, C, = C, = C,, ry; < (18) will yield correlations equal in size 
-0.5; if C, is greater than C, and C,, and opposite in sign to those given by 


', < —o.5. Thus if the absolute meas- equations (15) and (4). 













SOME PETROGRAPHIC RATIO 
CORRELATIONS 


In this section some practical instances 
of petrographic ratio correlation are re- 
viewed. The use of ratios is very common 
in petrography, and, although we often 
neglect to compute the actual correla- 
tion, we almost invariably arrange the 
paired values in tables or graphs which 
accept all the hazards but share few of 
the advantages of results obtained by 
computation. The examples cited here 
are of uneven grade and have in common 
only the property that each has come to 
my attention during the course of the 
past few months. In some cases computa- 
tions have been carried out, in others 
not. In some instances the usual inter- 
pretation of results passes from correla- 
tions of ratios to inferences about. abso- 
lute measures, while in others due allow- 
ance is normally made for the algebraic 
consequences of forming the ratios. 


THE GENERAL CASE: TWO RATIOS WITHOUT 
COMMON TERM 

G. T. Prior suggested in 1916 that in 
the stony meteorites the ratio MgO/ 
Nonmetallic le was negatively correlated 
with the ratio Ni/Metallic Fe. He seems 
to have considered this negative ratio 
correlation evidence for his hypothesis 
that the oxidized iron of stony meteorites 
formed at the expense of metallic iron. 
If the oxidized iron formed in this 
fashion, there would under most circum- 
stances be a negative correlation between 
oxidized and metallic iron; hence I think 
it is fair to say that Prior’s argument 
passes directly from a correlation be- 
tween ratios to an inference about rela- 
tions between two of the absolute meas- 
ures from which the ratios are formed.‘ 


4 Perhaps it would be even fairer to say that 
what I am criticizing is my own—and I think the 
common—impression of Prior’s position. In 1913 
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The absence of a common term render 
the case rather complex, but factors in. 
fluencing the sign of the ratio correlation 
may be visualized by entering equation 
(2) with X¥, = MgO, X, = Nonmetallic 
Fe, X; = Ni, X, = Metallic Fe. From 
the fact that r,, is negative Prior infers 
that r., is also negative. From equation 
(2) it is clear that 7,, is only one of four 
terms (containing combinations of eight 
variables) governing the sign of r,.. Now 
there is, in fact, a very high negative 
correlation between the two forms of 
iron (r,, = —0.68 as computed from the 
table given by Brown and Patterson 
[1947a]), but equation (2) shows that r,, 
might be negative regardless of the value 
of r.,. Only if the sum of the first three 
terms was equal to or larger than zero 
could it be argued that 7,, < o required 
fa < & 

The passage from knowledge about r,, 
to inference about r,, would be on some- 
what firmer ground if we had reason for 
supposing that the sum of these first 
three terms was not negative; but the 
simplest unexceptionable way to deter- 
mine the relation between XY, and X, is 
to compute r,,. In a situation as complex 
as this, the sign of the ratio correlation 


he argues (p. 35) that except for variations in the 
amount of nickel-iron, chondrites are, on the whole, 
characterized by “remarkable similarity, if not specif- 
icidentity” in chemical and mineral composition. In 
1916 (p. 26) he points out that “... the less the 
amount of nickel-iron in chondritic stones, the richer 
it is in nickel and the richer in iron are the magne- 
sium silicates’; but in his classification “. . . the main 
principle ...is the ratio of MgO to FeO in the 
magnesium silicates and in its interrelation with the 
richness in nickel of the nickel-iron.”’ In 1919 (p. 35° 
he clearly restates his belief in the oxidation hy- 
pothesis as an explanation of the dual ratio correla- 
tion (which in the language of this paper amounts 
to inferring a relation between absolutes from 4 
relation between ratios), regards the actual amounts 
of iron only as “support to this view,” and remarks 
that the amount of nickel in the metallic portion 
of meteorites is “. . . the factor which is seen to de- 
termine many of their other characters.” 
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is not a safe index of the sign of any one 
of the absolute-measure correlations. 
Another instance of dual ratio correla- 
tion is Washington’s well-known com- 
parison of 
K,O 
K,0 + Na:O 


; MgO 
with MgO + FeO 
in analyses of igneous rocks (1915). 
Figure 1 is his graph of the relationship, 
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this interpretation is commonly dropped 
in favor of one suggesting that in 
“potash-rich rocks” the “magnesia 
ratio” is usually high, and conversely. 
For this purpose ‘“‘potash-rich” may 
mean either that the ratio of potash to 
soda is large, that the ratio of potash to 
the sum of soda and potash is large, or 
just that the rocks are high in potash; 
and “magnesia ratio” similarly may 
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Fic. 1 


The ratio of MgO to the sum of MgO and FeO as a function of the ratio of Na,O to the sum of 


Na,O and K,O in analyses of igneous rocks (from Washington). 


and from it he argues that “. . . in igne- 
ous magmas potassium and magnesium 
on the one hand, and sodium and iron on 
the other, tend to vary together.” This 
conclusion is open to a variety of inter- 
pretations. The first impression of most 
people seems to be that it asserts the 
existence of two positive correlations, 
one between potash and magnesia, one 
between soda and iron. Some readers feel 
that the correlations but not their signs 
are specified. After a little consideration 


mean either the ratio of MgO to FeO or 
the ratio of MgO to the sum of MgO plus 
FeO. Washington’s own interpretation 
conforms to the first impression de- 
scribed above, viz., “...soda not un- 
commonly tends to vary with the iron 
oxides, while potash shows similar rela- 
tion to magnesia.” The correlation of 
these two ratios, however, is a function 
of at least ten absolute-measure correla- 
tion coefficients, only four of which are 
not potentially quite independent of one 
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another. It is certainly unreasonable to 
expect that valid inferences about rela- 
tions between any particular pair of 
absolute measures may be deduced from 
the mere fact that the ratio correlation 
is positive. If the purpose is to determine 
whether or not potash and magnesia are 
positively correlated, their correlation 
should be computed directly from the 
data. If it is desired to check the effect of 
soda, iron, or some function of the two 
on this relationship, partial correlation 
should be used. 

In any case, no one of the interpreta- 
tions advanced here may be regarded as 


TABLE 1 


VALUES OF ry: FOR 14 ANALYSES OF APATITE 
IRON ORES FROM NORTHERN SWEDEN 


ic 


- Mim Monn Com /Con | Nim/Nin 
P.O, SiO. +0.04 0.71 0.80 0.78 
P,O;... —0.04 0.68 0.84 0.85 


favored by figure 1. The ratio correlation 
means exactly what it says—that the 
ratio of potash to the sum of potash and 
soda varies directly with the ratio of 
magnesia to the sum of iron and mag- 
nesia. The graph suggests that, where 
sizable amounts of all four constituents 
are present, the correlation is, in fact, 
very high, and the relationship may be of 
considerable use if it is supplemented 
by information about relations between 
absolute measures; such information, 
however, may not be inferred from the 
ratio correlation but must be obtained 
directly from the data. 

As a final example of this case some 
calculations recently given by Lander- 
gren (1948) are shown in table 1. In the 
table the subscript m means “‘magnetic,”’ 
nm means “nonmagnetic.” These results 
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are useful for two reasons. They shoy 
the utility of the general relations, pro. 
vided that the restriction on the size of 
the C values is observed, and they pre. 
sent an opportunity for suggesting that 
this condition will sometimes be violated 
in practice. Disposing of this latter point 
first, minor constituents may vary wide. 
ly in individual analyses, but, by defini- 
tion, their averages are small. Therefore. 
C for a minor constituent may often be 
large enough that its third and higher 
powers cannot be ignored. The approxi- 
mations given in the preceding section 
will therefore be unsatisfactory. 

On the assumption that C values are 
small enough here that equation (4) is 
still applicable, it may be suggested that 
Landergren’s interpretation of table 1 is 
open to question. He concludes that 

. the quantity of Mn, Co and Ni 
entering the magnetic fraction and the 
nonmagnetic fraction respectively is de- 
pendent on the . . . content of P and the 
ratio P,O,/SiO,.’’ Now the correlations 
in the first line of table 1 conform ex- 
actly to equation (2) of the preceding 
section, while those of the second line 
fall under the case described by equation 
(7). Agreement between the two lines 
of the table is very good, so that one is 
led at once to wonder under what cir- 
cumstances equations (2) and (7) might 
lead to similar results. Equation (7) is 
reached from equation (2) by deleting 
all terms involving X,; table 1 would be 
explained if all correlations involving 
SiO, were negligibly small and the co- 
efficient of variation of silica were con- 
siderably smaller than that of phosphor- 
ous pentoxide. 

It thus seems quite possible that silica 
is not sensibly correlated with amy of the 
absolute values, and it is hard to see why 
the various m/n ratios of table 1 should 
be regarded as functions of silica unless 
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this can be established directly from the 
data. Unfortunately, Landergren has not 
recorded the basic measurements, and 
some or all of the C values may well be 
high enough to invalidate arguments 
based on equations (2) and (7); the case 
isincluded as an illustration of the power 
of the method under more auspicious cir- 
cumstances. 


CORRELATION BETWEEN TWO FRACTIONS 
WITH COMMON DENOMINATOR 

The relation which first attracted 
Pearson’s attention and on which he 
based the definition of spurious correla- 
tion was the correlation between two 
fractions with common denominator. In 
it fall most of the major parameters of 
the continental chemical-petrographic 
systems. As a hypothetical example, sup- 
pose that the quantities of olivine, py- 
roxene, and garnet are uncorrelated in a 
series of contact marbles. Letting olivine 
be X,, pyroxene X,, garnet X,, we have 
lo = Ti; = Y23 = 0. If two ratios, Y and 
Z, are formed from each set of measure- 
ments by using any one variable as de- 
nominator over the entire group, then 
r,, > 0, regardless of which variable is 
used as denominator. If C, = C, = C;, 
r,, = 0.5. If the object is to obtain as 
high a correlation as possible, the variable 
with largest C, often the one with smallest 
mean, should be used as common term. 

Incidentally, comparison of equations 
(5) and (16) indicates that this rule 
holds, regardless of whether the common 
term is numerator or denominator. Fur- 
ther, comparison of equation (4) with 
equation (15) shows that the nonspuri- 
ous ratio correlation (e.g., that holding 
when any or all of the absolute-measure 
correlations differ from zero) will also be 
identical whether the common term is used 
as numerator or as denominator (see dis- 
cussion following eq. [16]). 
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The Niggli system affords excellent 
examples of ratio correlation, for each of 
its major parameters is the quotient of a 
molecular amount divided by the sum of 
several molecular amounts. It is not in- 
frequently supposed that the relations 
between these parameters give a fair 
approximation of the relation between 
molecular amounts, molecular propor- 
tions, or weight percentages of the oxides 
after which they are named. 

Figures 2 and 3 illustrate a practical 
example. In each figure a graph of the 
Niggli symbols alk and gz for a group of 
analyses is shown on the left, and a graph 
of normative (or+ab) and normative Q 
for the same analyses appears on the 
right. Normative or and ab may be re- 
garded as multiples of the molecular 
amounts or weight percentages of the 
alkalis, whereas normative ( is an esti- 
mate of the molecular amount or weight 
percentage of uncombined silica. Rela- 
tions between the Niggli symbols are 
very different from those exhibited by 
the normative parameters. I have shown 
elsewhere (1948) that some approxima- 
tion of the relation between modal 
quartz and alkali-feldspar is just about 
the most important information to be 
extracted from these analyses; evidently, 
the result obtained will depend on the 
system used for reduction of the chemical 
analyses. In this particular case the 
Niggli parameters are misleading in the 
extreme, though the correlations between 
Niggli symbols are quite as “real’’ and 
“true” as the correlations between nor- 
mative molecules. 


CORRELATION OF A RATIO WITH A 
THIRD VARIABLE 
The correlation of a ratio with a third 
variable is not a particularly striking 
case, for the ratio correlation will be zero 
if neither term of the ratio is correlated 











with the third variable. It is to be noted 
from equation (7) that the coefficient of 
variation of the third variable does not 
directly affect the correlation between it 
and the ratio. The connection between 
equations (2) and (7) has already been 
described, with Landergren’s computa- 
tions as an example. 

If one is simply attempting to sys- 
tematize unruly data by the use of this 
type of ratio correlation, equation (7) 
shows that, in general, the highest ratio 
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correlation will be found if the ratio js 
formed from the two variables exhibiting 
the strongest direct variation. 
Nockolds (1947) has recently sug. 
gested that in aplites the ratio oy 
(or + ab) varies directly with normative 
quartz, and the same conclusion wa; 
reached by Vogt (1930) with regard to 
granites some years ago. One would sup. 
pose that in a case as simple as this it 
would be possible to pass directly from 
ratio correlation to inferences about 
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absolute measures. Equation (7) indi- 
cates that this is not so; for the ratio cor- 
relation to be positive, it is necessary 
only that 7,;C, — 723C. > o. Both ,; and 
r,,;may be positive or negative. If r,, > 0, 
we may not argue, for instance, that a 
highly siliceous granite will, in general, 
be rich in or or that a granite low in 
silica will, in general, be rich in ab+-an. 
I ought to add that neither Vogt nor 
Nockolds makes such inferences from the 
mere fact of the ratio correlation. Many 
of us, however, do make statements of 
the sort, usually based on direct knowl- 
edge of the absolute measures; and in un- 
wary moments nearly all of us imagine 
that the ratio correlation in some way 
confirms or even controls the observed 
relations between the absolute measures. 

In a comprehensive study of the Oden- 
wald diorites, Erdmannsdérffer has used 
a ratio correlation of this type in an at- 
tempt to discover the relation between 


‘hornblende, biotite, and quartz. The ob- 


served quantities of these three minerals 
vary widely and irregularly in over fifty 
modal analyses, but Erdmannsdorffer 
points out that the ratio 4b//bi decreases 
with increasing quartz. No computations 
are given, but a graph of the data leaves 
no doubt that a significant and probably 
large proportion of the variability of the 
ratio may be described as a linear 
function of quartz content. From this, 
Erdmannsdérffer concludes that ‘Man 
erkennt deutlich, wie mit steigendem 
Quarzgehalt, die Menge des Biotits 
zunimmt.”’ This should not be regarded 
as in any sense a necessary consequence 
of the ratio correlation, though it is, of 
course, one of many possible explana- 
tions. That it is not the correct explana- 
tion is suggested by the fact that Erd- 
mannsdorffer was led to use ratios in the 
first place by the “. . . scheinbare Regel- 
losigkeit . . .”” of the modal analyses. 
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CORRELATION OF A RATIO WITH ITS 
OWN DENOMINATOR 


A common case in petrography is the 
correlation of a ratio with its own de- 
nominator; and our interpretations of it 
are often correct for the reason that in 
making them we almost always rely 
heavily on information about the abso- 
lute measures. For some curious reason 
this information, which is the key to the 
situation, rarely receives the same promi- 
nence as does the ratio correlation. 

In his excellent study of the Tas- 
manian dolerites, A. B. Edwards (1942) 
has recently described the negative cor- 
relation between MgO and the ratio 
FeO/MgoO. It is by now evident that this 
correlation tells nothing of the relation 
between FeO and MgO. Equation (10) 
shows that, even if there is no correla- 
tion between FeO and MgO, there will be 
negative correlation between MgO and 
the ratio FeO/MgO and that, if the co- 
efficient of variation of MgO is equal to 
or larger than that of FeO, the absolute 
value of this correlation will be equal to 
or greater than 0.71, MgO and FeO still 
being uncorrelated. This seems to be the 
case with Edwards’ analyses. In each of 
the sills studied, MgO has a very much 
wider range and a somewhat smaller 
average value than FeO. The ratio cor- 
relation is thus compatible with the no- 
tion that there is no systematic relation 
whatever between MgO and FeO. 

What gives meaning to Edwards’ data 
is the persistence with which high MgO 
values occur in the lower, and low MgO 
values in the upper, portions of each sill. 
MgO has surely been extracted from the 
liquid during crystallization, but with 
little if any effect on the FeO content of 
the residue. One of the unfortunate by- 
products of our fascination with ratios is 
that we pay insufficient attention to the 
absolute measures. The most prominent 














feature of one of Edwards’ oxide profiles 
is the sharp increase in alumina, match- 
ing an equally sharp decrease in MgO, 
but apparently unrelated to silica; and 
throughout all the analyses there seems 
to be an inverse variation of alumina and 
magnesia. Yet the interesting possibility 
that alumina proxies for magnesia in the 
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Prior and has recently been revived by 
Brown and Patterson (19474, b; 1948), 
to whom we are indebted for the first 
reasonably complete compilation of 
stony meteorite analyses in many years, 
The earlier work was based on rather few 
analyses, but the Brown-Patterson table 
contains fifty-three in the critical range 
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Fic. 4.—The ratio —— asa function of total metal in stony meteorites (from Brown and Patter- 


Total Metal 
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pyroxenes rather early in the crystalliza- 
tion is not considered. Something of the 
sort would seem to be almost inevitable 
unless color index increases considerably 
and systematically with depth. 

Another example of this case is the 
correlation of Ni/Total Metal with Total 
Metal in stony meteorites. First dis- 
covered by Flight (1887), this relation- 
ship appears again as an introduction to 
the dual ratio correlation proposed by 





and makes it possible to establish the 
relationship as something more than a 
shrewd guess. These authors have com- 
puted average Ni/ Metal values for 2 per 
cent intervals of metal content and have 
plotted the average values of the ratio as 
a function of its denominator. Their 
graph of the relationship is shown here as 
figure 4. As most of the metal content is 
iron, one would expect a similar relation 
to hold between Ni/Metallic Fe and 
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Metallic Fe, and this proves to be the 
case. 

For the analyses listed in the Brown- 
Patterson table (1947a, p. 405), the co- 
eficient of variation of Ni is 0.49, and 
that of Fe is 0.60; under the circum- 
stances equation (g) shows that nothing 
can prevent the occurrence of negative 
correlation between Ni/Metallic Fe and 
Metallic Fe; for practical purposes the 
range of this ratio correlation is o 2 
ap Bee fe 

In forming ratios we are always care- 
ful to use quantities occurring in the 
same analysis; in forming, say, the 
eighth ratio of a group, we divide the Ni 
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cent Metallic Fe’ and repeated the 
averaging calculations previously carried 
through with the true sample pairs. The 
results are shown in table 2. The tend- 
ency for the average ratio value to de- 
crease with increase in Metallic Fe is 
quite as strong in the randomized pairs 
as in the true sample pairs. The explana- 
tion is not far to seek. The dispersion of 
Ni values is very much smaller than the 
dispersion of Metallic Fe. A little con- 
sideration will show that, when the range 
of the numerator is so much smaller than 
that of the denominator, the fraction is 
bound to decrease with increase in the 
denominator, regardless of the relation 


TABLE 2 


AVERAGE N1/METALLIC FE VALUES FOR 2 PER CENT INTERVALS OF FE 
TRUE AND RANDOM PAIRS 


I 3 
True sample pairs 0.122 0.136 
Randomized pairs 0.689 0.167 


value reported in the eighth analysis by 
the metallic Fe value in the same 
analysis, and the comparison is between 
the resulting fraction and its denomina- 
tor. What equation (10) tells us is that 
this precaution is quite unnecessary. The 
result would be the same, to the required 
precision, if we just as carefully avoided 
using the true sample pairs. As an illus- 
tration, I purposely randomized the Fe 
and Ni values in all the Brown-Patterson 
analyses containing less than 14 per 

5A C value of o.5 is really too high to permit 
rigorous application of the approximations given 
in the preceding section; C3 is 0.125, which may not 
be neglected if one wants an r value good even to the 
first digit. The demand for negative correlation is 
in this case so extreme, however, that the general 
outline of the argument is applicable, though the 
exact value of the correlation no doubt differs from 
that yielded by eq. (9). 


MID-POINT OF Fr CLAss 


5 } 7 9 Ir 13 
0.122 0.116 0.090 0.0907 0.086 
0.183 0.126 0.101 0.072 0.056 


between numerator and denominator. As 
a corollary, we ought not to expect to dis- 
cover much about this relationship from 
the ratio correlation. As a matter of fact 
there is very strong positive correlation 
between Ni and Metallic Fe in the 
analyses listed by Brown and Patterson. 


CORRELATION BETWEEN TWO RATIOS 
WITH COMMON NUMERATOR 
It has been shown that correlations 
between two ratios with common nu- 
merator will be the same as those found 


® By numbering the Ni values in order, assigning 
to each Metallic Fe value one of the first fifty-three 
nonidentical two-digit numbers drawn from a page 
of Tippett’s table of random numbers, ranking the 
results, and combining each Ni value of order i with 
the Metallic Fe of rank i, where i = 1 — 53, in- 
clusive. 
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when the common term is the denomina- 
tor; the discussion under “correlation 
between two fractions with common de- 
nominator” of the present section applies 
directly to this case. 


CORRELATION OF A RATIO WITH 
ITS NUMERATOR 
Correlations of a ratio with its 
numerator will be identical in size with 
those between ratio and denominator, 
but they will be of opposite sign. Be- 
tween equations (12) and (g) there is the 
same parallelism as has previously been 
noted between equations (18) and (s) 
(see discussion following eq. [19]). 
Correlation of a ratio with its numera- 
tor is quite commonly used in geochem- 
ical work, as, for instance, in the com- 
parison of Li/Mg with Li by Strock 
(1936) or of Ge/Si with Ge by Gold- 
schmidt and Peters (1933). Even if the 
ratio correlation were high, which is 
usually far from the case, it would not 
justify a conclusion that there was or was 
not some systematic relationship be- 
tween the absolute measures. If the ab- 
solute measures are uncorrelated, the 
ratio correlations will always be positive 
unless the numerator is constant. In 
geochemical work, particularly on minor 
elements, dispersion will often be large 
in relation to mean, so that the algebraic 
approximations given here will be unsatis- 
factory. This does not mean that the 
ratio correlation is independent of the 
absolute-measure statistics and is of 
course no reason for supposing that the 
latter are determined by the former. For 
the most part, satisfactory interpreta- 
tions of ratio correlations will still hinge 
on adequate information about relations 
between the variables from which the 
ratios are formed. 


CORRELATION OF TWO RATIOS WHEN THE Dp. 
NOMINATOR OF ONE IS THE SAME AS THE 
NUMERATOR OF THE OTHER 

With regard to location of the common 
term, equation (18) is just equation (16) 
or equation (5) with all signs changed. 
Thus, if Y = gz of the Niggli system and 


Z=alk, ty, = My = Ty If one 


were studying the relations between 
MgO, Metallic Fe, and Nonmetallic Fe in 
meteorites by means of ratios, he would 
find that this equality also held where, 
for instance, 


SUMMARY AND CONCLUSIONS 


When four variables are formed into 
two ratios, there will be no correlation 
between the ratios if there is no correla- 
tion between the variables, and there 
will, in general, be correlation between 
the ratios where there is correlation be- 
tween one or more pairs of the absolute 
measures if these pairs include a term 
from each ratio. But the size of the ratio 
correlation will, in general, differ from 
that of any one of the absolute-measure 
correlations, and its sign will provide no 
indication about the sign of any one of 
these. 

When three variables are formed into 
two ratios by using one variable as a 
common term, the size of the correlation 
will be the same whether the common 
term appears as the numerator in both 
ratios, as the denominator in both ratios, 
or as numerator in one and denominator 
in the other. In the first two cases the 
sign of the correlation will also be the 
same, but in the third it will be reversed. 

When three variables are compared 
by forming any two into a ratio, the cor- 
relation of this ratio with the third 
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variable can be zero only if each term of 
the ratio is uncorrelated with the third 
measure or if the terms of the ratio are 
merely coded values of each other. The 
ratio correlation is not a reliable criterion 
of the absolute-measure correlations. 

When a ratio is formed of two absolute 
measures and compared with one of the 
two, the sign of the correlation will de- 
pend on whether the common term ap- 
pears as numerator or as denominator, 
but its numerical value will be the same 
in either case. 

Whenever there is a common term in- 
volved, the ratio correlation cannot be 
zero if the absolute measures are uncor- 
related. In each specific case its sign is 
fixed; in all cases its size is determined 
by the relative variability of the absolute 
measures. 

There are very few exceptions to the 
general rule that no ratio correlation per- 
mits valid inference about the relation- 


ship between any two of the absolute 


measures from which the ratios are 
formed. On the other hand, a knowledge 
of the sample statistics of the absolute 
measures will permit valid deduction of 
the correlation between any set of ratios 
formed from them, provided only that 
the C values involved are sufficiently 
small. 

In descriptive petrography we com- 
monly resort to ratios when the absolute 
measures seem feebly correlated or un- 
correlated, and we frequently employ 
one or another of the common-term rela- 
tionships; so it is not surprising that we 
have turned up a number of remarkably 
high ratio correlations. A good deal of the 
significance that we attach to these cor- 
relations stems from the premise that in 
other samples or in other arrangements 
the same absolute measures, still uncor- 
related, might yield very much smaller 
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ratio correlations or even correlations of 
opposite sign. On occasion, most of us are 
not at all averse to passing from correla- 
tion between ratios to inferences about 
absolute measures. On all three counts 
we often lead ourselves, as well as 
others, astray. 

The misuse of ratios in specific prob- 
lems is easily corrected. The difficulties 
encountered in deciding how to record 
large numbers of measurements are not 
so readily eliminated; for it may be quite 
plausibly argued that the justification of 
the ratio is here purely that of a book- 
keeping device. Thus, in at least two in- 
stances (Niggli, 1923; Simonen, 1948) the 
Niggli symbols of large numbers of rock 
analyses have been recorded in place of 
the weight percentages or molar amounts. 
To recapture from these tables enough 
information for a study of the relation- 
ship between any pair of constituents 
would require almost as much work as 
the compilation of the table, and for 
some constituents (e.g., ferrous and 
ferric oxides) it would be impossible. 
Similarly, one can discover the composi- 
tion of any stony meteorite in the 
Brown-Patterson table, but only by un- 
doing a computation previously per- 
formed by Brown and Patterson. 

It is very easy to carry an argument of 
this sort too far. Sometimes hypotheses 
are framed in terms of ratios, and some 
phenomena are understandable only as 
ratios. I do not mean to offer a blanket 
condemnation of all ratio correlation, 
but it does seem to me that in most of 
the problems of descriptive petrography 
our primary concern is still with the 
simplest variables obtained by direct 
measurement or observation. There may 
be instances in which ratio correlations 
accomplish more for us than absolute- 
measure correlations do, but this will 
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usually occur only where our knowledge 
of the relations between the absolute 
measures is satisfactory. In any case, the 
ratio correlations are contained in the 
absolute measurements, whereas the con- 
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verse is not usually true. Even when at. 
tention is focused on ratios, the absolute 
measures ought to be recorded, and cer. 
tainly for purposes of tabulation the 
absolute measures are preferable. 
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SOME EXAMPLES OF THE APPLICATION OF THERMO- 
CHEMISTRY TO PETROLOGY' 


TH. G. SAHAMA’ AND D. R. TORGESON? 


ABSTRACT 


A brief summary is presented of measurements of heats of solution of minerals belonging to the forsterite- 
fayalite and to the enstatite-orthoferrosilite series and of artificial ilmenite and geikielite. 

“In the olivine and orthopyroxene series, the heat of solution in HF is found to be a linear function of the 
Mg: Fe ratio, indicating perfect isomorphism between the corresponding end-members. The importance and 
possibilities of applying calorimetry to the study of isomorphism are emphasized. On the basis of the calori- 
metric data available for the minerals in question, the heat and free energy are given for the reaction: olivine + 
quartz > 2 pyroxene. The difference of the stabilities of MgSiO, and FeSiO, is illustrated. 

The stability reaction of ilmenite and geikielite in the presence of olivine or orthopyroxene is calculated 
from the calorimetric data. The influence of the entropy of mixing isomorphic minerals upon the stability re- 


lation is emphasized. 


INTRODUCTION 

Besides the direct petrological obser- 
vations on natural rocks and minerals 
and experimental work on synthetic com- 
pounds, there are still other ways of in- 
vestigating the formation and stability 
relations of minerals. One of these meth- 
ods, indirect in itself, is offered by ther- 
mochemistry. Despite the fact that quan- 
titative thermochemistry, based on nu- 
merical data, apparently has a wide field 
of possible application to petrology, it 
has been used very little as yet for solv- 
ing petrological and mineralogical prob- 
lems. The reason for such a circumstance 
is revealed by the fact that calorimetric 
data for natural or artificial minerals are 
scanty. Therefore, it is to be hoped that, 
in the future, more calorimetric measure- 

'This is a report of work conducted by co- 
operation of the Geophysical Laboratory of the 
Carnegie Institution of Washington and the Pacific 
Experiment Station of the Bureau of Mines, United 


States Department of the Interior. Manuscript 
received September 30, 1948. 


*Professor of geochemistry, University of 
Helsinki, Finland, and at present visiting investi- 
gator of the Carnegie Institution of Washington, 
associated with the staff of the Geophysical Labora- 
tory. 


’ Formerly chemist, Pacific Experiment Station, 
Bureau of Mines. 


ments on minerals will be undertaken by 
scientists interested in problems dealing 
with the stability relations of minerals 
and mineral parageneses. A successful 
application of the thermodynamic prin- 
ciples in petrology evidently is not pos- 
sible without having available reason- 
ably comprehensive data obtained by 
experimental calorimetry and covering 
all the more important rock-forming 
minerals. 

From the viewpoint of the petrologist, 
the most important thermal quantities 
for minerals are, perhaps, the heats of 
formation, owing to the fact that, in a 
large number of reactions possible be- 
tween minerals, the largest contributor 
to the free energy usually is the heat of 
the reaction. The most suitable method 
of ascertaining heats of formation of sili- 
cates and of other rock-forming minerals 
is offered by solution calorimetry. This 
method is applicable to all minerals that 
are rapidly soluble in hydrofluoric acid 
or other solvents. If this method is used 
at elevated temperatures—say, near 75° 
C.—a considerable number of petrologi- 
cally important minerals can be investi- 
gated. 

By means of hydrofluoric acid solution 





255 











256 


calorimetry, the silicates of the olivine 
and orthopyroxene series were investi- 
gated. Previously, Mulert (1912) had 
found that pure magnesium ortho- and 
metasilicate did not dissolve rapidly 
enovgh ‘to allow a calorimetric heat-of- 
solution measurement. In the present in- 
vestigation it was found that, if the ma- 
terial is very finely powdered and if the 
experiment is made at an elevated tem- 
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cal a brief summary of the mineralogica| 
and petrological results of the work and 
to point out some possible applications of 
thermochemical methods and considera. 
tions to problems of petrology and min. 
eralogy. 








FORSTERITE-FAYALITE SERIES 





The samples of the minerals belonging 
to the forsterite-fayalite series, used for 



























TABLE 1 
HEATS OF SOLUTION OF THE FORSTERITE-FAYALITE 
SERIES IN 20.1 PER CENT HF AT 73°7 C. 
: Molecular AH* Uncertainty 
Sample 
Composition (Cal/ Mol) (Cal/Mol) 
. oa , = 
A. Forsterite, Dreis, Eifel, Germany. . a = | —94,250 160 
| 
B. Forsterite, Ultenthal, Tyrol, Austria = a —93,590 | 80 
C. Chrysolite, Marjalahti, Finland 2 —t | —93,930 | 140 
D. Hyalosiderite, Lake Leistilinjirvi, Nakkila, Fin- ee | 
land sas ‘0 50.0} —89,070 | 100 
; Fa 44.0f | as, 
FE. Hyalosiderite. Susimiki. V: a. Fink: Fo 54.2) | | 
». Hyalosiderite, Susimiki, Vampula, Finland. . Fa 45.8) | —88,820 | 160 
d = 4 
| | 
F. Ferrohortonolite, Tunaberg, Sweden... ” 5:9. | —82,080 | 160 
| a o4.1f | 
| | | 
— Fo 0.0 
“a. Synthetic fayalite. . -{ 5 
Fa. Synthetic fayalite | Fa 100.0 81,330 | 50 
“ alent L | 
* Throughout this paper the symbol AH corresponds with the Lewis and Randall usage. 


perature, the rate of the reaction can be 
increased sufficiently for an accurate de- 
termination. The determinations were 
made at 73°7 C. In addition, some heat- 
of-solution measurements of ilmenite and 
geikielite were made. The design of the 
calorimeter and the procedure employed 
has been given previously by Torgeson 
and Sahama (1948). The details of the 
results have also been reported (Sahama 
and Torgeson, 1948a and b). However, 
it seems pertinent to give in this periodi- 


the present investigation, are listed in 
table 1. In addition, the molecular com- 
positions, calculated on the basis of the 
chemical analyses as pure mixtures of 
Mg,SiO, and Fe,SiO,, and the measured 
heats of solution in 20.1 per cent HF 
(averages of the individual runs), ex 
pressed in calories per mol at 73°7 C., are 
given in the second and third columns, 
respectively. The fourth column indi- 
cates the uncertainty of the average. All 
heat-of-solution values are corrected for 








‘ical 
and 
1S of 
eTa- 
nin- 


ging 
| for 


od in 
com- 
f the 
es of 
sured 
| HE 

ex- 


imns, 
indi- 
e. All 











the impurities present in the samples. In 
designating the natural minerals, the 
nomenclature proposed by Deer and 
Wager (1939) was employed. 

The results of the heat-of-solution 
measurements listed in table 1 are plot- 
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ted against the composition in the upper 
part of figure 1. As is seen from the figure, 
the hdat of solution of the forsterite- 
fayalite series is, within the experimental 
errors, a linear function of the Mg:Fe 
ratio. The dashed line indicated in the 
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figure gives the best average and can be 
expressed by the following equation: 


AH = — 81,330 —140.5X,,, , 


in which X,,, is the molal percentage of 
Mg,SiO, in the corresponding isomor- 
phic mineral. An extrapolation to pure 
Mg.SiO, yields 

AH Mg,  adlians 95,380 cal/mol . 


The linearity of the relationship be- 
tween the heat of solution and the 
Mg:Fe ratio in the forsterite-fayalite 
series is mineralogically interesting. It 
indicates that the replacement Mg-Fe in 
the olivine structure is thermally perfect. 
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urements is that the perfectness of the 
isomorphism can be tested by calori- 
metric means. It shows the possibility of 
applying calorimetry to the study of iso. 
morphism of natural silicates and min- 
erals. It is to be expected that many 
other isomorphic series among the nat- 
ural minerals will show deviations from 
linearity, indicating changes in the struc- 
ture of the mixture as compared with the 
pure end-members. 


THE ENSTATITE-ORTHOFERROSILITE 
SERIES 


The samples of the enstatite-orthofer- 
rosilite series employed, with molecular 


TABLE 2 


HEATS OF SOLUTION OF THE ENSTATITE-ORTHOFERROSILITE 
SERIES IN 20.1 PER CENT HF AT 73°7 C. 


Sample 


G. Enstatite, Kjérstad, Bamle, Norway 
H. Eulite, Hirvisaari, Rantasalmi, Finland. 


K. Eulite, Tunaberg, Sweden 


In other words, the heat of isomorphic 
mixing is zero in this series. Structurally, 
it means that the isomorphic substitu- 
tion Mg-Fe in the olivine structure oc- 
curs without causing strain. 

The structural interpretation of the 
results of the heat-of-solution measure- 
ments expressed above contains in itself 
nothing new. All petrological observa- 
tions of natural rocks, supported by the 
experimental investigations carried out 
by Bowen and Schairer (1935), indicate 
that the isomorphism between forsterite 
and fayalite is complete without any de- 
tectable complications. The importance 
of the result of the heat-of-solution meas- 





Molecular All Uncertainty 
Composition (Cal/Mol) (Cal/Mol 
En 96.5 

> , — 62 ,67 8 

Fs 3.5 aon . 
En 27.1 

— 58,57 5 

Fs 72.9 5®,579 5° 
En 12.9 | 

> —57,9 2 

Fs 87.1f | ew ° 

| 


| | 
compositions, heats of solution in 20.1 
per cent HF (averages of the individual 
runs), and uncertainties of these aver- 
ages are given in table 2. The values are 
corrected for the impurities. The nomen- 
clature follows Poldevaart (1947). 

The heat-of-solution values given in 
table 2 are plotted against composition 
in the lower part of figure 1. The three 
points for the enstatite-orthoferrosilite 
series in the figure seem to indicate that 
in the orthorhombic pyroxene structure 
also the replacement Mg-Fe is perfect. 
The average shown by the dashed line in 
the figure is represented by the equation 

AH = — 57,010 —58.7X ng, 
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in which X,,, is the molal percentage of 
MgSiO, in the mixture. An extrapolation 
to the pure end-members yields 


— 62,880 cal/mol , 


AH gsio, 
— 57,010 cal/mol . 


AH FeSiO, 


In addition to the samples listed in 
table 2, a bronzite from Eerolamminvu- 
ori, Paistjarvi, Heinola, Finland, was in- 
vestigated (sample L). Corresponding 
to its mode of occurrence as a segregation 
ina labradorite rock, its alumina content 
is, according to the analysis published by 
Lokka (1943), unusually high (6.39 per 
cent Al,O,;). Its composition is thus out- 
side the enstatite-orthoferrosilite series 
proper. The mineral shows En:Fs = 
70.8:29.2 and gives a heat of solution 
(20.1 per cent HF at 73°7 C.) of AH = 
—62,500 + go cal/mol. This value is 
1,390 calories above the dashed line in 
figure 1. 

The deviation of the heat of solution of 
the Heinola bronzite from the heat-of- 
solution curve of the pure MgSiO,- 
FeSiO, series can be qualitatively ex- 
plained. As pointed out by Lokka (in his 
paper cited above), aluminum apparent- 
ly replaces both quadrivalent silicon and 
divalent iron and magnesium. Now, if 
trivalent aluminum enters into the struc- 
ture in place of quadrivalent silicon and 
divalent iron and magnesium, the elec- 
trostatic neutrality will, in general, be 
disturbed. Even if the structure as a 
whole remains neutral, individual oxygen 
ions might not all be equally balanced. 
The deviation from strict neutrality of a 
structure is a phenomenon that is very 
well known in some other minerals 
(sphene, andalusite, sillimanite, etc.) of 
comparatively simple chemical composi- 
tion. In minerals and mineral groups (am- 
phiboles, pyroxenes, micas, etc.) in which 
the isomorphic substitutions are high- 
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ly complicated, considerable deviation 
from neutrality might occur. Qualitative- 
ly, it is well understood that an overbal- 
ancing or underbalancing of the individ- 
ual oxygen atoms weakens the structure. 
Of course, the behavior of the Heinola 
bronzite represents just one singleexample 
of such a supposed weakening of the struc- 
ture causing a high heat of solution. There- 
fore, tog much weight should not be as- 
signed to it. However, it shows in prin- 
ciple the possibility of applying calori- 
metric methods to the study of the ener- 
gy balances connected with the isomor- 
phic substitutions in minerals. 


STABILITY OF ORTHOFERROSILITE 


As is well known through the investi- 
gations carried out by Bowen and Schai- 
rer (1932, 1935), pure orthoferrosilite is 
not stable at moderate and at high tem- 
peratures. In the lithophysae of some ob- 
sidians from different localities Bowen 
(1935) has found monoclinic pyroxenes 
with optical properties corresponding to 
those extrapolated for pure FeSiQ,. 
However, as pointed out by him, the 
natural occurrence of ferrosilite in ortho- 
rhombic or monoclinic form cannot yet 
be regarded as established. Whether or 
not FeSiO,, as a crystalline compound 
with pyroxene structure, is thermody- 
namically stable is a question that can- 
not be answered with certainty at this 
time. 

On the basis of the data for the ferro- 
magnesian ortho- and metasilicates given 
above, the stability relation of enstatite 
and orthoferrosilite can be illustrated. 
Since the heat of solution of quartz in 
20.1 per cent HF at 73°7C. has been 
measured by the authors (AH = —33,- 
ooo + 20 cal/mol), the heats of the fol- 
lowing reactions can be given: 

2. Fe.SiO, + SiO, > 2FeSiO; . 
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The entropies at 298°16 K. are known 
for Mg.SiO, and MgSiO, (Kelley, 1943), 
for SiO, (Kelley, 1941a), and for Fe,SiO, 
(Kelley, 19416). For FeSiO, the entropy 
can be roughly predicted as S},sio, = 
3 cal/deg/mol. Accordingly, the heats 
and free energies of the two reactions 
mentioned above are, at 298°16 K.., 


1. AH = —2,620 + 280 cal. , 
AF° = — 2,500 + 320; 


2. AH = —310 + 210, 
AF° = —460. 


Because of the great uncertainty in the 
entropy increment of reaction 2 and, on 


TABLE 3 


HEATS OF SOLUTION OF FETIO;, MGTIO,, FE,S1O,, 
AND OF THE EIFEL OLIVINE IN AN ACID MIX- 
TURE CONTAINING 10.05 PER CENT HF AND 
10.05 PER CENT HCL AT 73°7 C. 


as AH Uncertainty 

, (Cal/Mol) Cal/Mol) 
Synthetic ilmenite. . — 46,650 60 
Synthetic geikielite — 56,820 40 
Synthetic fayalite —75,370 30 


Forsterite, Dreis, Eifel, 


Germany —go,530 | 120 


the other hand, of the small heat of that 
reaction, no reliance can be placed in the 
sign of the free energy of that reaction. 
The difference in the stability of MgSiO, 
and FeSiO, is, however, evident. 


STABILITY RELATION OF ILMENITE 
AND GEIKIELITE 


The abundance relation of ilmenite 
and geikielite very strongly favors ilmen- 
ite. Geikielite is extremely rare and has 
been found so far only in the gems of 
Ceylon. 

To reveal the thermodynamical sta- 
bility relation of ilmenite and geikielite 
in the presence of olivine or orthopyrox 
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ene, the heat-of-solution measurements 
listed in table 3 were undertaken. Be. 
cause completion of the reaction of 
MgTiO; with pure HF proved to be too 
slow, an acid mixture containing 10.0; 
per cent HF and 10.05 per cent HC! was 
employed. The reported heat-of-solution 
values represent averages of the individ- 
ual runs. 

On the basis of the heat-of-solution 
values given in table 3 for synthetic 
fayalite and for Eifel olivine, the heat of 
solution of the forsterite-fayalite series in 
the acid mixture can be expressed by the 
following equation: 

AH = — 75,370 —165.2X nj. 


Extrapolation to pure Mg,SiO, yields 
AHmgsio, = — 91,890 cal/mol . 


‘Lhese heat-of-solution data yield, for 
the reaction Fe.SiO, + 2MgTiO, - 
Mg.SiO, + 2FeTiO,, the heat, AH = 
— 3,820 + 250 calories at 298716K. 
The entropies for Fe,SiO, and Mg,SiO,, 
as given by Kelley (19410, 1943), and for 
FeTiO, and MgTi0,, as given by Sho- 
mate (1946), yield AF = —4,510 cal- 
ories as the free energy of the reaction at 
298°16 K. This shows that the associa- 
tion Mg,SiO, + 2FeTiO, is thermody- 
namically stable at room temperature. 
The negative free energy persists also at 
higher temperature ranges. 

To approach the conditions prevailing 
in natural rocks, the following general re- 
actions were considered. 

1. n(Fe;,o0o—2-Mg,) 2SiO, + MgTiO; 
~>+ n(Fey,o0-pMg,) 2SiO, + FeTi0s, 


2. n(Fey.o0-.Mg,) SiO;+ MgTiO, 
—> & (Fey .o0 yMg,) SiO, -- Fel 10; : 


In these reaction equations, m and x are 
independent variables, and p depends 
stoichiometrically on m and x. 

Under the condition that the iso- 
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morphism of the forsterite-fayalite and 
the enstatite-orthoferrosilite series is 
thermally perfect, as indicated by the 
heat-of-solution data, the heats (AH) 
and the heat-capacity increments (AC,) 
of the two reactions are calculable and 
are independent of the numerical values 
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andC,(FeSiO,) =C,(Fe.SiO,) —C,(FeO). 
On the other hand, the entropy incre- 
ments (AS) of reactions 1 and 2 depend 
largely on the values of » and x, because 
of the entropy of mixing of the corre- 
sponding isomorphic mixtures of the 
olivines and orthopyroxenes. The en- 
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Fic. 2.—Stability relation of geikielite and ilmenite in the presence of olivine (1) or orthopyroxene (B) 
ina temperature range of 400°-1400° K.; x, mol fraction of magnesium in olivine (or pyroxene); », molecular 
ratio of olivine (or pyroxene) to geikielite (or ilmenite) in the rock. The dotted area indicates the stability 
lield of geikielite in the presence of olivine or pyroxene. 


of n and x. The following high-tempera- 
ture heat capacities were adopted: 
MgTiO, and FeTiO,, given by Naylor 
and Cook (1946); Fe,SiO, and MgSiO,, 
given by Kelley (personal communica- 
tion); Mg,SiO, and FeSiO,, estimated by 
means of the following approximations. 
C,(Mg.SiO,) = C,(MgSiO,) + C,(Mg0) 


tropy of mixing was calculated according 
to the relationship, ASyixing = —RN, In 
N, — RN, In N,, in which N, and N, are 
the mol fractions in the isomorphic min- 
erals in question. 

The free energies for the two general 
reactions 1 and 2 were computed from 
the data. The computations were made 
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according to the relationship d(AF/T) = 
—(AMH/T?)dT. The thermodynamic no- 
menclature follows Lewis and Randall 
(1923). The results of the computations 
are shown diagrammatically in figure 2, 
A and B. The stability fields of geikielite 
in association with olivine or orthopy- 
roxene, respectively, are indicated by the 
dotted areas on the individual sections 
corresponding to the different tempera- 
tures from 400° to 1,400° K. 

Some of the thermal data employed in 
the above computations are approxima- 
tions that may deviate more or less from 
the true values. Therefore, figure 2 can 
be applied only qualitatively or semi- 
quantitatively. However, it shows that 
the association geikielite + olivine (py- 
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roxene) represents a high-temperature 
paragenesis. In the presence of olivine or 
orthopyroxene, geikielite becomes stable 
only at higher temperatures and in rocks 
very rich in magnesium. This result indi- 
cates that geikielite may occur in those 
rocks more frequently than is known at 
present. Its proper identification in rocks 
is difficult and is possible only through 
separation and chemical tests. 

Independently of the particular sta- 
bility of geikielite, the present investiga- 
tion shows that the entropy of mixing of 
an isomorphic mixture is a factor of geo- 
logical importance and one which, in 
suitable instances, may greatly affect the 
stability relations of the minerals and 
mineral parageneses. 
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ISOSTASY—EXTENDED' 


ROSS GUNN 
Division of Physical Research, United States Weather Bureau 


ABSTRACT 


Basic principles are established for the equilibrium of a strong elastic lithosphere supported by a weak 
underlying magma. The lithosphere acts to spread out localized geological loads over a considerable area, 
and it is shown that isostatic equilibrium is achieved only when a load is reasonably uniform and distributed 
over a distance of some 330 km. This extended principle of isobaric equilibrium contains the isostatic prin- 
ciple as a special case. Basic mechanisms and methods for determining the deformation of the lithosphere 
under various kinds of loading are outlined. The geological significance of gravity anomalies is pointed out, 
and methods of forecasting their magnitude and distribution in terms of the deformation of the lithosphere 


are given. 


The extended theory of isostatic equilibrium that the author has called ‘‘isobaric equilibrium” finds 
application in all quantitative discussions of the deformation of the lithosphere. 


INTRODUCTION 

In a series of papers appearing for the 
most part in the Journal of the Franklin 
Institute and in Geophysics, the author 
has presented an extension of the prin- 
ciples of isostatic equilibrium (1937, 
1943@, b, 1944, 1945, 1947). This exten- 
sion deals particularly with the part 
played by a strong elastic lithosphere and 
is especially significant in areas exhibit- 
ing crustal deformation. Because these 
early papers were mathematical and not 
entirely understood by many readers, it 
is the object of the present paper to re- 
view the concepts and present them in a 
simpler and more usable form. For for- 
mal justification of some of the conclu- 
sions outlined herein, the reader may 
therefore find it desirable to refer to the 
original publications. 

Modern ideas concerning the equilibri- 
um of the lithosphere were introduced in 
1854, when J. H. Pratt (1855) prepared 
an extensive report on the figure of the 
earth and gave an explanation for the ap- 
parent deficiency in the gravitational at- 
traction of mountain masses as measured 
in closely adjacent areas. Shortly there- 
after, George B. Airy (1855), in discuss- 


‘Manuscript received July 8, 1948. 


263 


ing the discrepancy brought out by 
Pratt’s work, wrote: “It appears to me 
that the state of the earth’s crust lying 
upon the lava may be compared with 
perfect correctness to the state of a raft 
of timber floating upon water; in which, 
if we remark one log whose upper surface 
floats much higher than the upper sur- 
faces of the others, we are certain that 
its lower surface lies deeper in the water 
than the lower surfaces of the others.”’ 
The various ideas of Airy and Pratt 
have been examined in the light of many 
observations, and these have established 
the principle of isostasy, at least for 
large, relatively undeformed crustal 
areas. William Bowie (1931), one of the 
active workers in this field, defines iso- 
static equilibrium as 
the condition of rest which the outer materials 
of the earth tend to acquire. The equilibrium 
is between the crust and the subcrustal material 
and not within the crust itself [italics by the au- 
thor]. It is based on the idea that the mass of 
each unit section of the earth’s crust exerts the 
same pressure on the subcrustal material. Ac- 
cording to the principle of isostasy, if the earth’s 
crust were cut into prisms of the same cross-sec- 
tion by imaginary vertical planes, the prisms 
would have the same mass if the isostatic 
condition were perfect. 
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This formulation of the isostatic prin- 
ciple, which specifically ignores the 
strength of the lithosphere, has guided 
geological interpretation for thirty-five 
years. The principle, as stated, is an in- 
complete formulation of crustal equilib- 
rium, and the whole problem needs re- 
consideration. 

As a general rule, the differences be- 
tween the calculated and the observed 
values of the acceleration due to gravity, 
or the isostatic anomalies, in undeformed 
and level country and over most of the 
oceans are relatively low, being equiva- 
lent, on the average, to a layer of rock 
about 400 feet thick. However, in recent 
years, as methods for correcting for ele- 
vation, density, and topography have be- 
come more accurate and detailed data 
from the field have become available, it 
has become increasingly clear that many 
regions of the earth cannot possibly be 
in isostatic adjustment. Large gravity 
anomalies have been interpreted as indi- 
cating gross departures from normal den- 
sity in the vicinity of the observing sta- 
tions and frequently have led to mislead- 
ing statements suggesting that the anom- 
alies imply an unstable geological struc- 
ture, even though clear-cut observational 
evidence is available to show that the re- 
gion is entirely stable and that no dia- 
strophic activity is in progress. Such dis- 
crepancies in the interpretation of geo- 
logical structure in regions of gravity 
anomaly became so glaring, particularly 
in relation to theories of mountain-build- 
ing, that the author was led to review and 
extend isostatic theory (1937-1947) with 
the object of bringing the common 
gravity irregularities into the framework 
of natural physical law. This rationaliza- 
tion was accomplished by including an 
aspect of crustal stability specifically ig- 
nored in Bowie’s foregoing definition of 
isostasy by quantitatively formulating 
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the effect of finite strength and stresses 
within the crust itself. These neglected fac. 
tors are found to be of fundamental im- 
portance. 


IDEAL ISOSTASY 


Consider a newborn earth made up of 
rock material at such a temperature that 
its strength in shear is negligible. After 
circulation and mixing have acted for a 
few million years, the underlying layers 
of such an earth will approach the hydro- 
static state characterized by (a) a steadi- 
ly increasing temperature and density 
toward the center; (6) the coincidence of 
surfaces of equal density, equal pressure, 
and equal gravitational potential; and 
(c) a surface layer of cool and strong 
rock, whose density is notably less than 
the underlying hotter rocks. 

Suppose that great prisms of cool rock 
are ultimately formed at the surface by 
selective crystallization or some other 
means and that these aggregates are sys- 
tematically less dense than the underly- 
ing weak magma. Obviously, they will be 
supported according to Archimedes’ 
hydrostatic principle and will float in 
such a way that the weight of each prism 
will be exactly equal to the weight of the 
displaced underlying magma.? Consider 
a number of prisms supported in the 
above manner, all of different but com- 
parable density. Assume that each prism 
is self-supported and that no friction 
whatever exists between the lateral faces. 
Then, if a reference prism is selected and 
if the acceleration due to gravity is 
measured on it and if suitable corrections 
are made for the distance of the observ- 
ing station above some reference level 
and for the figure of the earth, it will be 
found that the measured values on all 

2 The term ‘“‘magma” is here used to denote any 


type of underlying rock having inappreciable 
strength in shear. 
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| prisms will not appreciably differ. This 
> js a consequence of the laws of gravita- 
’ tional attraction and potential theory, 
| which show that to a sufficient approxi- 


mation the vertical acceleration due to a 
horizontal plate of considerable extent 


> isgiven by the familiar relation 








(1) 


where Ag is the change in gravitational 
acceleration due to the plate, M is the 
mass per unit area of the plate, and y¥ is 
the gravitational constant. If the plate is 
not extensive, geometrical corrections 
must be applied, but equation (1) will 
represent the principal term. Since, by 
Archimedes’ principle, the mass per unit 
area of the floating block is, on the aver- 
age, exactly equal to the mass per unit 
area of the displaced liquid and since by 
equation (1) the gravitational defect de- 
pends on this quantity, it is clear that 
diferent thicknesses or densities of a 
floating block will not produce important 
changes in the measured value of the ac- 
celeration due to gravity. 

Now, if a layer of sediments be piled 
on top of one of the prisms by transfer 
from some other place, the prism will 
sink until it displaces a weight equal to 
that transferred to its surface; and, if 
sediments are washed off the prism, it 
will rise in a strictly analogous manner. 
Thus, as long as the prisms are perfectly 
free to move vertically in relation to one 
another, the basic ideas expressed by 
Airy apply. Because such a system of 
freely floating blocks will show a con- 
stant gravitational acceleration irrespec- 
tive of the sedimentary load placed there- 
on, we may describe the state as one of 
ideal isostasy. The principal term of the 
gravity anomaly for such configurations 
is zero. Although the commonly used 
statement that the mass per unit area 
above a given equipotential surface is 


Ag=2ryM , 
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constant is sufficiently accurate for most 
cases of ideal isostasy, yet occasionally 
the exact statement must be employed 
(Gunn, 1936), namely, that the mass per 
unit area multiplied by mean gravity in the 
prisms, is @ CONSTANT. 


DEPARTURES FROM IDEAL ISOSTASY 


In order to understand the basic mech- 
anisms, consider a simple experiment 
with the frictionless prisms mentioned 
above. It has been noted that adding or 
subtracting loads of sediments to the 
prisms will not influence the gravitation- 
al acceleration as measured at a standard 
level on these prisms. Suppose, now, that, 
instead of a load of sediments, a large 
force of unspecified origin acts down- 
ward, which will displace from under- 
neath the prism a certain mass of magma 
whose weight by Archimedes’ principle 
is exactly equal to the applied downward 
force. Because no material is transferred 
to the block and because mass has been 
forced out from under it, it is clear that a 
deficiency of mass exists in this prism in 
relation to the initial condition of ideal 
isostasy, and therefore a negative gravity 
anomaly will be measured. The magni- 
tude of this gravity anomaly is given by 
equation (1), where M is the mass per 
unit area of the magma displaced by the 
downward-applied force. In a similar 
way a vertically upward force might be 
applied to an adjacent prism, resulting in 
an inflow of excess mass below the prism 
in such a way that a positive gravity 
anomaly would be measured on this 
prism. Although we cannot conceive of 
such disconnected geological forces that 
will depress or lift frictionless prisms of 
the earth’s crust, it is entirely possible 
for vertical forces to be applied to each 
prism arising from the frictional contact 
with and motion of adjacent prisms. 
Such forces could be very important. It 
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is worth while, therefore, to go one step 
further and assume that, instead of the 
prisms being frictionless, they are actual- 
ly in elastic and cohesive contact with 
their neighbors, forming a great aggre- 
gate of prisms and constituting an elastic 
sheet having great internal strength. 
Frictional and cohesive forces between 
the vertical faces now prevent the indi- 
vidual adjustment of the prisms. For ex- 
ample, select a central prism and assume 
that sediments are piled on it alone. If 
this prism stood alone and was supported 
only by the weak asthenosphere, it would 
sink in the substratum and displace a 
weight of material exactly equal to that 
added. However, because of the assumed 
forces of friction or cohesion acting on 
the vertical faces of the prism, this can- 
not take place. The applied load at- 
tempts to force the prism down, but now, 
because it is mechanically associated 
with the neighboring prisms, the latter 
are forced down while the original prism 
is supported by shearing forces to exactly 
the same extent. Thus, under the selected 
prism not enough mass is forced out lat- 
erally to support it from the bottom, and 
part of the support is supplied by vertical 
shearing forces appearing on the vertical 
faces of the prism. Conversely, as a result 
of such forces, considerable material is 
forced out from under the adjacent 
prisms, even though no load of sediments 
is applied to them. Of course, considering 
the aggregate of prisms, just enough 
material is displaced to carry the super- 
posed load of sediments on the original 
prism, and a specialized condition of re- 
gional isostasy is achieved. 

The deficiency or excess of mass dis- 
placed from under any prism can general- 
ly be calculated in terms of the applied 
force or the vertical stress. For simplici- 
ty, consider a prism having frictionless 
walls and of unit cross-section supported 
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ideally by the substratum. Suppose this 
prism is pushed down by some arbitrary 
vertical force acting on unit area or by 
shearing forces arising from the finite 
strength of the lithosphere. Let this 
force per unit area be the vertical result. 
ant supporting stress, S». Such a stress, 
acting on a freely moving prism, will dis. 
place substratum material until a restor. 
ing stress, equal to the weight of the dis. 
placed liquid, restores equilibrium, or 


S,=M’'g, (2) 


where M’ is the mass per unit area of the 
displaced material and g is the accelera- 
tion due to gravity. It is clear, therefore, 
that a downward or positive resultant 
stress from any source whatever results 
in a deficiency of mass under the prism, 
while an upward stress results in an ex- 
cess of mass. If this change in mass is as- 
sumed to be in the form of a large flat 
sheet, the principal term of the accom- 
panying gravitational anomaly, Ag, in- 
creases with the applied stress; or, spe- 
cifically by equations (1) and (2), this is 
given by 


ghg=2ryS,. (3 


The frictional or cohesive forces acting 
on the vertical faces of the prisms neces- 
sarily produce downward (or upward 
forces in adjacent prisms, and_ these 
prisms are just as much in equilibrium as 
if they were supported ideally from the 
bottom. It follows from the above con- 
siderations that, if the resultant vertical 
stress distribution can be determined for 
any arbitrary load on the lithosphere 
either by measurement or by calculation, 
then it is possible to evaluate the gravity 
anomalies associated therewith. Con- 
versely, a knowledge of the space distri- 
bution of gravity anomaly and mass per- 
mits an estimate of the vertical super- 
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sed stresses and the resulting deforma- 
tions of the lithosphere. 

Because the integrated upward and 
downward stresses over a section of the 
lithosphere more than 400 km. long must 
add up to zero, it follows that the space 
integral of the principal term of the 
gravity anomaly over the same path 
must also approach zero. 


THE LITHOSPHERE AND ITS SUP- 
PORTING MAGMA 


The lithosphere may be defined as the 
cool, strong, rocky layer of the earth that 
is supported by the denser, weaker, and 
relatively hotter underlying layers. The 
upper few kilometers of the lithosphere 
are traversed by many shallow faults and 
fractures. As a unit, these upper layers 
are relatively weak, but at a depth of 
perhaps 5 km. and more, where the com- 
pressive stress due to the superincum- 
bent rocks exceeds the crushing strength, 
local fractures of the type regularly oc- 
curring at the surface itself become more 
or less self-healing, and this part of the 
lithosphere may be considered to be 
equivalent to a strong, uniform elastic 
sheet. By “strong,” the author means 
that the crustal material will support 
shearing stresses of 1 or 2 X 10° dynes/ 
cm? for periods of time comparable to the 
age of mountains. 

It is a convenient mathematical fiction 
to assume that the strong lithosphere ex- 
tends down to a definite level and that 
below this level the supporting magma is 
fluid in nature and has zero strength in 
shear. Actually, the transition from the 
strong lithosphere to the weak magma is 
continuous. The existence of high moun- 
tains and ocean deeps shows that rocks 
just below the surface are strong and will 
support large shearing stresses for long 
periods of time. But as the depth in- 
creases, the temperature rises steadily; 
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and, even though it is possible that the 
shear strength for short impulses is 
maintained, yet at sufficiently elevated 
temperatures progressive plastic de- 
formation will almost certainly be 
produced by stresses applied for a suffi- 
ciently long time. Too little is known of 
the mechanical properties of great sheets 
of rock, but there is accumulating evi- 
dence to show that the apparent shearing 
strength of geological materials decreases 
as the time of application of the stresses 
increases and that both the strength and 
the time interval necessary to induce 
appreciable creep decrease notably as the 
temperature increases. Thus the steady 
increase in temperature as one goes deep- 
er and deeper into the crust implies a 
continuous change in the apparent 
strength of the rocks. An examination of 
the continuous transition to the weak 
magma in comparison with a discontinu- 
ous one shows that, except for some un- 
certainty in the position of the neutral 
axis, an elastic outer surface of the earth 
so constituted behaves mechanically in a 
manner remarkably like a strong, uni- 
form, elastic lithosphere of perhaps 50- 
km. thickness, supported by and super- 
posed upon a weak magma. A constant 
that measures the ease of deformation 
of the lithosphere (characteristic con- 
stant = 0) is found to be by far the most 
important quantity describing its me- 
chanical characteristics, and this con- 
stant is related to an equivalent thick- 
ness that becomes the true thickness only 
when the transition from the lithosphere 
to the asthenosphere is assumed to be dis- 
continuous. It is sufficiently accurate, 
therefore, to assume that the lithosphere 
is a strong, elastic sheet of sufficient uni- 
formity that it obeys well-established 
laws of elasticity and that it is every- 
where supported on a weak underlying 
magma that, given adequate time, will 
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adjust itself to hydrostatic equilibrium. 
This time necessary for adjustment is 
short in most cases of geological interest 
and probably approximates 25,000 years. 


ISOBARY—A BASIC EQUILIBRIUM PRIN- 
CIPLE EMBRACING ISOSTASY AS 
A SPECIAL CASE 


Because it is assumed that the litho- 
sphere is everywhere supported on a sub- 
stratum having negligible strength in 
shear, it is possible to establish the cri- 
terion for vertical equilibrium in any se- 
lected prism. In the asthenosphere below 
the lowest section of the lithosphere, 
some surface of equal pressure or equal 
gravitational potential is adopted as a 
reference. If an applied vertical load on 
the lithosphere should increase the pres- 
sure at this selected level, hydrostatic 
laws show that transient horizontal pres- 
sure gradients will be set up, and these 
will force the weak rock outward more or 
Jess horizontally until the constancy of 
pressure over the selected equipotential 
surface is re-established. It is a property 
of the hydrostatic state that, when con- 
stancy of pressure is re-established, the 
system is in equilibrium. Now the pres- 
sure at this selected level is made up of 
three components, namely: (a) the 
weight per unit area of the column of 
weak rock extending from the reference 
level up to the bottom of the lithosphere; 
(b) the weight per unit area of the col- 
umn that extends from the bottom of 
the lithosphere to the actual surface; and 
(c) the vertical stress, S,, imposed on the 
column by the finite strength of the 
lithosphere. The total pressure arising 
from components a and b may be con- 
veniently written Mg, where M is the 
mass per unit area lying above the refer- 
ence equipotential surface and 7 is the 
average gravitational acceleration in any 
selected prism. 
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The criterion for isobaric equilibrium, F 


therefore, is 
P=Mg+S,=Constant, (4) 


: 4 
where P is the pressure at the selected 


equipotential or isobaric level and S, js 
the superposed vertical stress. This rela. 
tion differs from the analogous relation 
for isostasy by the inclusion of the ver. 
tical superposed stress, S,. The relation 
is generally applicable, even though each 
prism may be elastically jattached to its 
neighbors. The quantity S, is convenient. 
ly taken to be positive when it acts in the 
direction of g or downward. This exact 


result may be stated in words as follows: | 


For sections of the lithosphere to be in 
stable or isobaric equilibrium, it is neces- 
sary and sufficient, in prisms reaching 
from a selected surface of equal pressure 
in the asthenosphere to the actual sur- 
face, that the vertical force per unit area 
arising from the finite strength of the 
lithosphere, added to the product of the 
mass per unit area and the average ac- 
celeration due to gravity, shall be a con- 
stant. 

When the lithosphere is adjusted to 
this condition, one may conveniently say 
it is in “isobaric equilibrium,” since the 
word “isobaric” implies equal pressure. 
This extended formulation of isostatic 
equilibrium has been called “isobaric 
equilibrium” following a suggestion of 
C. E. Dutton, who originally discarded 
this more descriptive term in favor o 
“isostasy.” Isobary, therefore, corre- 
sponds to a state of vertical mechanical 
equilibrium, even though it permits real 
anomalies of gravity to exist. In level re- 
gions where vertical stresses are negli- 
gible, isobary obviously reduces to the 
more restricted equilibrium of isostasy. 
It is convenient to pick some standard 
column as a reference, which should pref- 
erably be near the center of a large unde- 
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formed region. At such a place one may 
be assured that the superposed vertical 
stresses usually due to deformation are 
small and the regional isostatic adjust- 
ment approaches the ideal, quite irre- 
spective of the finite strength of the 
lithosphere. Now, if gravitational meas- 
urements at other places are compared 
with those of such a reference point, 
making appropriate corrections for its 
position in relation to the geoid and the 
figure of the earth, it is clear that gravity 
anomalies will be an indication of vertical 
stress. Equation (3) may be combined 
with equation (4), and the principle of 
isobaric equilibrium may be written in 
the approximate but convenient form, 


P=Mg- BAS _ Constant, (5) 


ary 

where Ag is the principal term of the 
gravity anomaly. The exact anomaly 
may be calculated if the geometrical dis- 
tribution of stress, mass, and topography 
is available in any given case. The prin- 
ciple of isobaric equilibrium as here ex- 
pressed, together with observed numeri- 
cal values for the gravity anomalies, 
provides a valuable tool for the study 
and determination of lithospheric defor- 
mations, particularly in special geological 
structures that are definite exceptions to 
the principle of isostasy as usually de- 
fined. 


DEFORMATION OF A HYDROSTATI- 
CALLY SUPPORTED ELASTIC 
LITHOSPHERE 

In order to achieve a basic familiarity 
and understanding of the types of equi- 
librium that exist when a mountain is 
supported by isostatic forces and when it 
is supported by isobaric forces, it may be 
illuminating to introduce an analogy at 
this point. Suppose a small boy walks out 
on a deep pond covered with thin ice. If 
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the ice is too thin, he breaks through, 
and, according to Archimedes’ principle, 
he is then supported in the water only to 
the extent of the weight of water that he 
displaces by his body. Quite obviously, 
the support is in accordance with hydro- 
static principles, and, assuming he has a 
normal density, complete support is not 
achieved until he is almost immersed. 
Now, if the thickness of the ice increases, 
a thickness is finally reached which will 
support the boy. This thickness corre- 
sponds to the condition that the weight 
of the boy is distributed over a big 
enough area that the actual depression of 
the ice below its normal level is sufficient 
to displace a weight of water exactly 
equal to the weight of the boy. It might 
be said that, when the ice is strong 
enough, its downward depression by the 
boy produces a very shallow circular boat 
which displaces enough water to support 
him effectively. This downward deflec- 
tion of relatively thin ice is more or less 
apparent to every skater, and it is clear 
that the strength and elasticity of the ice 
plays an important part in the support of 
the skater. Now, as the thickness of the 
ice increases, its flexural rigidity also in- 
creases, and the support for the weight of 
the boy is spread out over a greater and 
greater area. Simultaneously, the mag- 
nitude of the systematic downward defor- 
mation of the ice due to the weight of the 
boy is so reduced that the weight of dis- 
placed water is maintained exactly equal 
to the weight of the boy. It should be 
clear, further, that the weight of the 
skater is never passed back to the sup- 
porting banks of the pond but is always 
passed downward into the supporting 
and underlying fluid. Whether the boy is 
on the ice or in the water is immaterial in 
relation to the amount of displaced un- 
derlying fluid, but it is clear that the dis- 
tribution of the displacement in the two 
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cases is radically different. It should be 
evident that the strength of the ice plays 
an important role in modifying the ge- 
ometry of any given system of loads. The 
analogy is such that the boy floating in 
the water corresponds to the floating 
equilibrium type of support envisaged by 
isostasy, whereas the distribution of dis- 
placement of the supporting water and 
the stress in the ice sheet correspond to 
isobaric equilibrium. In the language of 
isostasy, the “compensation” for the 
boy’s head, which is above water, is to be 
found immediately below in the deficien- 
cy of density of the boy’s body. On the 
other hand, when the boy is on top of the 
ice sheet, the compensation for his weight 
is spread out over a considerable area by 
the strength of the ice and is to be found 
in the water displaced and pushed aside 
by the downwardly deformed ice sheet. 
Attention is drawn to the fact that the 
ice even some distance from the boy is 
bent downward well below its equilibri- 
um position, and the mass of water in a 
typical unit area is accordingly less than 
normal. Thus, if a finite strength is as- 
signed to the lithosphere, most of the 
compensation for a given load is not 
under the load but to the sides. This is in 
sharp contrast to the basic assumptions 
of isostasy. 

Analysis shows that, with knowledge 
of the thickness and elastic modulus for 
the ice and the density of the underlying 
fluid, it is possible to calculate, using 
straightforward engineering principles, 
the magnitude of the downward defor- 
mation of the ice for all points immedi- 
ately under the boy and in his vicinity. 
It is found that the volume of the depres- 
sion caused by the boy, multiplied by the 
density of the water, is exactly equal to 
the mass of the boy. This is true for any 
thickness of ice. The dimensions of the 
depression and the effective radius over 








which the boy’s weight is distributed can 
also be determined. 

Although the above basic ideas under- 
lying isobary are very simple, it is evi- 
dent, in order to calculate the vertical 
supporting stresses arising from a super- 
posed load on the actual lithosphere, that 
some rather detailed calculations must be 
carried out. The deflection produced in 
the lithosphere by a single large vertical 
force may be computed with relative ease 
and is strictly analogous to the case of 
the small boy standing on a sheet of ice. 
There are practical mathematical diffi- 
culties involved in computing the down- 
ward deflection of the lithosphere for 
loads distributed over a considerable 
area. An investigation has shown, how- 
ever, that if the distributed load is linear 
in character and extends for some hun- 
dreds of kilometers in such a way that 
the load distribution can be well repre- 
sented by a single section taken at right 
angles to the linear trend, then the down- 
ward deformation for a number of dis- 
tributed loadings may be easily worked 
out. Now, because many very important 
earth features like deeps, mountain 
chains, and continental borders are es- 
sentially linear in character and extend 
for a thousand kilometers or so, this 
mathematical simplification is adequate 
for most geological interpretations and 
enables one to work out in certain areas 
the exact deformed figure of the litho- 
sphere. 

In order to determine the figure for 
any given loading, it is convenient to 
select a typical section of unit width at 
right angles to the linear trend and con- 
sider this slice of the lithosphere as a 
strong elastic beam resting on, and ulti- 
mately supported by, the weak underly- 
ing magma. It is assumed that Young’s 
modulus of elasticity and the sectional 
moment of inertia per unit width is sub- 
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stantially constant throughout the ex- 
plored region. The reader is referred to 
the original papers for the details of such 
calculations. It may be stated, however, 
that, once the load applied to the original 
undeformed lithosphere is specified or 
known, it is practicable by the solution of 
the basic equations to determine at once 
(a) the downward or upward displace- 
ment of the lithosphere, relative to its 
position prior to the application of the 
load, for all points along the section; 
(b) the distribution of the internal 
bending moment and the vertical shear; 
(c) the vertical restoring force per unit 
area, S,, due to the deformation of the 
lithosphere, which is, in fact, proportion- 
al to the space gradient of the vertical 
shear; (d) the fiber stress at every point 
in the equivalent beam, and (e) the prin- 
cipal term of the gravity anomaly result- 
ing from the internal supporting stresses. 
It is evident that all the foregoing items 
are of considerable importance in struc- 
tural studies of the lithosphere; and be- 
cause these may be determined, once the 
exact figure of the lithosphere is known, 
practical methods for its evaluation seem 
important. 

The deformation of the lithosphere 
considered as an elastically supported 
beam was first carried out by M. Smolu- 
chowski (1909), who particularly con- 
sidered solutions for the deformation on 
the assumption that the crust was com- 
posed of successive layers of rock and 
that its flexural rigidity was very much 
lower than is now believed to be the case. 
The roots of the basic equations for this 
situation lead to periodic solutions that 
this author thinks could have been de- 
scriptive of the lithosphere as a whole 
only in pre-Cambrian time, when the 
crust was very thin and weak. The roots 
of the basic equations employed in the 
author’s investigations are those corre- 
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sponding to a homogeneous and elasti- 
cally strong lithosphere, and these lead 
to logarithmically attenuated waves, 
which well describe actual geological for- 
mations. It can be shown (Gunn, 1937- 
1947) that the solution to the basic equa- 
tions having important geological appli- 
cations leads to vertical deformations, 
y, of the form 


y = Pe~* sin bx + Qe~™ cos bx , (6) 


where, in most practical cases, a is indis- 
tinguishable from the characteristic con- 
stant b, while P and Q are amplitudes de- 
pending on the imposed boundary condi- 
tions and the mechanical characteristics 
of the lithosphere. The vertical deforma- 
tion y given by equation (6) may be rep- 
resented for any assumed vertical force 
by a simple curve, as shown in figure 1. 
In figure 1, A, the downward deforma- 
tion of the lithosphere due to a down- 
ward force is shown for an unbroken sec- 
tion and in figure 1, B, for the same force 
applied to one end of a broken section. 
The abscissae are expressed in terms of 
an angle measured in radians, which is 
the product of the characteristic constant 
b and the distance x along the beam from 
the applied force. In the figure actual 
distances are also plotted for the earth. 

It is important to note from this figure 
that in the actual earth the load applied 
at a specific point appreciably deforms 
the lithosphere for a distance of some 250 
km., quite irrespective of how large that 
force may be. As the force increases, the 
downward deformation increases propor- 
tionately, and the resulting deformation 
is increased at every position along the 
beam by the same factor. The distance 
from the force to the point on the litho- 
sphere where the downward deformation 
is zero is independent of the force and is 
measured only by the flexural rigidity of 
the beam and the density of the underly- 
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ing magma. These, in turn, determine the 
value of the characteristic constant b that 
measures the ease of deformation of the 
lithosphere. 

The actual deformation produced by 
any type of load can be evaluated by re- 
placing the force, F, in figure 1, A, by an 
elementary load and summing up all the 
downward displacements caused by each 
element. Thus, if a selected load of any 
configuration is divided up into a series 
of narrowly defined loads and the down- 
ward displacement of each is calculated, 
the total downward deformation is evi- 
dently the sum of all the deformations, 
due allowance always being made for the 
exact position of each element. Certain 
geometrical forms are easily calculated, 
and a number of these have been pre- 
sented in graphical form for reference and 
use (Gunn, 1937-1947). Figure 2 shows 
the deformation of the lithosphere that 
results from a long block load of width 
2 bm, height #, and density d, as calcu- 
lated by the above method. The figure 
brings out the fact that, as the width of 
the load increases, the block forces the 
lithosphere deeper and deeper into the 
underlying magma until at a width of the 
block equal to something like 4 bm the 
block is essentially submerged to a level 
that approaches the isostatic condition. 
The degree of support provided by the 
lithosphere and its influence on isostasy 
and gravity anomalies as expressed in 
terms of the size of the block will be dis- 
cussed in a later section. In order that the 
plotted data may be of maximum useful- 
hess, in figure 2 and in the following fig- 
ures distances along the beam are meas- 
ured in radians (6x), so that the curves 
may be generally applied to any section 
of the lithosphere corresponding to dif- 
ferent values of }, its characteristic con- 
stant. The normal value of b in the 
earth’s crust is such that the distance for 
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one radian as plotted in this series of 
curves closely approximates 120 km., so 
that an approximate distance scale may 
be simultaneously plotted on each curve. 
For example, in figure 2, the semi-width 
bm = } corresponds to a width equal to 
one-fourth of 120 km., or 30 km.; and we 
see that such a block, of total width 60 
km., is deformed downward into the un- 
derlying magma by a distance equal to 
0.24 of the equivalent height of the block. 
Thus, by knowing the equivalent height 
of a geological load having the form of a 
long block and the appropriate numerical 
factor as determined from figure 2, the 
exact downward deformation in kilome- 
ters may be immediately estimated. 
The author’s papers give a large num- 
ber of figures showing the downward de- 
formation of the crust produced by vari- 
ous loads. It is worth while to note that 
almost any load can be simulated with 
sufficient accuracy by the addition of 
these standard deformations. These de- 
formations superpose, and it follows from 
the mathematical form of the equations 
that the sum of the lithospheric deforma- 
tions due to each load is the same as the 
deformation for the sum of the loads. 


THE SIGNIFICANCE AND MAGNITUDE OF 
THE CHARACTERISTIC CONSTANT 0 


The flexural rigidity of the lithosphere 
and the density of the underlying weak 
magma, together with the actual value of 
acceleration due to gravity, determine 
the magnitude of the characteristic con- 
stant of the lithosphere in such a way 
that, when its value 0 is raultiplied by the 
distance x along the effective beam, the 
resulting product is the angle bx ex- 
pressed in radians. The quantity 6, there- 
fore, is the. reciprocal of the relaxation 
distance of the lithosphere, and this may 
be defined as the approximate distance 
required for gross deformations to be ap- 
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preciably attenuated. This characteristic 
distance (1/6) approximates 120 km., 
and 6 is (8.4 + 2) 10-8 cm. 

The numerical value of b may be cal- 
culated in a straightforward manner in 
terms of the roots of the basic equations 
which define 6 if the lithosphere is reason- 
ably homogeneous. More specifically, 6 is 
equal to the fourth root of the product of 
the density of the underlying magma 
(diminished by the density of any sedi- 
ments or water that washes into the de- 
pression) multiplied by the acceleration 
due to gravity and divided by four times 
the flexural rigidity of the lithosphere, or 


es 
4 (flexural rigidity) 
7 _ 4) & | 
. £F F 


(7) 


where £ is Young’s modulus of elasticity, 
T is the effective thickness of the strong 
lithosphere, d, the density of the under- 
lying weak magma, d the density of 
water or sediments that may completely 
fill any downward deformation, and g the 
acceleration due to gravity. Thus 6 is 
large for a very weak and yielding crust 
or one in which plastic flow occurs, but it 
approaches zero for a crust that is per- 
fectly rigid. The quantity 6, therefore, 
measures the ease of deformation of the 
crust by a given load. The correct values 
for the above quantities are not exactly 
known. For example, the effective value 
of Young’s modulus in the earth’s crust 
according to competent investigators has 
been taken as low as 2 X 10" dynes/cm? 
to as high as 2 X 10” dynes/cm’. The 
thickness of the crust has similarly been 
given as from 30 to go km. Taking the 
outside limits above and assuming d, = 
3.1 gm/cm3, d = 0, and g = 980 cm 

sec’, it is easy to show from equation (7) 
that the smallest probable value for 0 is 
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5.45 X 10° cm and the highest j 
20.2 X 10-* cm~. For the present, on 
may take E as close to 1.2 X 10” dynes 
cm? and JT = 50 km. Substitution 9 
these in equation (8) leads to b = 8.4 x 
10-* cm™". 

Because the above direct calculation of 
the characteristic constant 0 is open ty 
uncertainty and because its value is oj 
first-rate importance in all discussions of 
the stress and deformation of the litho. 
sphere, it is of importance to examine 
various known geological structures to 
see how observed b changes under widely 
varying terrestrial conditions. In order 
to secure reliable values, it is necessary 
to select for analysis simple geological 
structures in which the acting forces are 
both known in position and localized in 
nature. As an example, consider the 
downward deformation produced by a 
localized load, as shown in figure 1. The 
sediments initially deposited at the 
single mouth of a great river might well 
represent such a concentrated load, and, 
instead of these sinking at the point oi 
deposition, thus establishing _ isostatic 
equilibrium, the strength of the litho- 
sphere supports the sediments locally and 
passes the load to adjacent regions. Now, 
as more sediments are deposited by the 
river, they are spread out in a great delta, 
producing a downward deformation like 
that of figure 1, A, and having a total 
width of bx = 37/2. This basin of widely 
deformed ocean bottom encourages the 
deposition of sediments at a considerable 
distance from the original mouth, anda 
great, symmetrical delta is built up. The 
mechanical behavior of the crust suggests 
that sedimentation will proceed laterally 
until the downward deflection is filled, at 
least thinly, to the edge of the original 
depression. When such a delta is well de 
veloped, as in the case of the Nile and 
Niger rivers (Barrell, 1914), the extent 
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and figure of the delta describes the ex- 
' tent and approximate figure of the under- 
) \ying downwardly deformed lithosphere. 
Barrell’s data show that the total width 
for both the above-mentioned deltas is 
about 600 km., from which it is easy to 
calculate that b is equal to 7.85 X 107° 
cm, or the relaxation distance 1/b = 
127 km. 

As is shown in the basic papers listed, 
it is possible to calculate the principal 
term of the gravity anomaly when the 
deformations of the lithosphere are 
known. The detailed characteristics of a 
highly compressed segment of the litho- 
sphere that has failed in shear along some 
zone of weakness have been worked out, 
and the expected distribution of gravity 
anomaly therein has been calculated. 
These calculations (Gunn, 1937-1947) 
show that the overriding and under- 
thrusting of the sheared ends produce 
positive anomalies in the overriding sec- 
tor and negative anomalies in the under- 
thrusting sector that extend for a dis- 
tance away from the point of failure 


equal to an angular distance, bx = 7/2. 
Beyond this negative zone in the under- 


thrusting section is a positive area ex- 
tending for a distance of bx = 7. By 
the use of observed data obtained by 
Vening Meinesz in the Javanese area, } 
may be immediately calculated, and it is 
found that 6 for this area approximates 
86 X 10°* cm. Gravity data from a 
number of points on the earth, when they 
are summarized and allowance is made 
lor geometrical variation and uncertain- 
ty, lead to a value for b = (8.4 + 2) X 
io* cm. Although the suggested values 
ib do not appear to change greatly over 
the surface of the earth, this matter de- 
erves special investigation and may 
thed considerable new light on the 
mechanical characteristics of the crust. 
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VERTICAL SUPPORTING STRESS PROVIDED 
BY THE STRONG LITHOSPHERE—ISO- 
BARIC STRESS FUNCTION 


Because the underlying weak magma 
always approaches the hydrostatic state, 
the equilibrium of the lithosphere under 
an applied load is consummated when 
the pressure along a given gravitational 
equipotential surface is constant. Thus, 
on comparing the total pressure in a typi- 
cal undeformed area with that where the 
lithosphere is notably distorted, one 
needs only to equate the expressions for 
the pressure at the equipotential surface 
in the two areas. The pressure in general 
is made up of two components, (a) the 
weight of the superincumbent rocks, in- 
cluding both those of the lithosphere and 
the load (where it exists), and (6) the 
vertical stress supplied by the strong and 
elastic lithosphere. If the lithosphere is 
bent downward by a given load, the pres- 
sure under the load at the originally se- 
lected reference level is exactly equal to 
the pressure some distance away in an 
undeformed area. It may be shown that, 
under this condition, 


— S, = hdg ont ydog = dog (H —¥) ; (8) 


where // is the “equivalent height” at 
any position, y is the downward deforma- 
tion immediately below, and S, the verti- 
cal supporting stress contributed by the 
lithosphere. The equivalent height, H, is 
the height that any given load would 
have if it had the same density as the 
underlying magma. It is defined by 
na-% 


9 
dy : ' ) 


where d and / are the density and actual 
height of the load, respectively, and d, is 
the density of the underlying magma. 
The relation of equation (8) states, sim- 
ply, that the downward force per unit 
area due to the applied load, /dg (or, bet- 











ter, Hd.g), under conditions of equilibri- 
um is exactly equal to the upward force 
per unit area due to the downwardly dis- 
placed magma, yd,.g, added to the upward 
force per unit area (—S,) due to the 
strength of the lithosphere. Since 5S, is 
taken positive in the direction of g (the 
acceleration due to gravity), S, directly 
under a load is usually negative. The re- 
lation expressed by equation (8) is in- 
valuable in all discussions of stresses in 
the lithosphere and gravity anomalies 
and is called the “isobaric stress func- 
tion.’”” By methods beyond the scope of 
this review, it is possible to express the 
vertical stress, S,, in terms of a high-order 
curvature of the lithosphere rather than 
in terms of the load and displacement. It 
was shown in the earlier sections leading 
to equation (3) how the principal term of 
the gravity anomaly can be expressed in 
terms of the vertical stress so that one 
may finally write the multiple basic rela- 
tions derived in the series of papers here 
outlined (Gunn, 1937-1947) as 


EI d‘ty_ Ag * 


H 7 y iz dog dx4 “ 2rydy Les dog 


(10) 
In these relations, one must not confuse 
the vertical stress, S., with the fiber 
stresses within the lithosphere that deter- 
mine the internal elastic equilibrium but 
should note, rather, that —S, is the sup- 
porting force per unit area contributed by 
the strength of the lithosphere itself. 
The above principles may be applied 
to give specific answers to perplexing 
problems concerning isostasy and the oc- 
currence of gravity anomalies. lor ex- 
ample, it has heretofore been impossible 
to state just how big a block load carried 
by the crust must be, in order for it to 
approach isostatic equilibrium. It should 
be clear from figure 2 that the downward 
deformation by a long block load ex- 
pressed as a fraction of its equivalent 
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thickness increases almost linearly with 
increasing width of the block up to about 
2 bm = 1 (or of total width about 129 
km.), while it is essentially submerged for 
widths greater than 2 bm = 4 (or of a 
total width of 480 km.). Thus block loads 
of limited extent are carried almost com- 
pletely by the strong lithosphere, and the 
ultimate support is passed downward to 
the weak underlying magma over a large 
area, whereas extensive loads, because of 
their considerable downward deforma- 
tion, are largely carried by the hydro- 
static forces provided by displaced un- 
derlying magma. By the use of curves 
like that plotted in figure 2, it is possible 
to determine the fraction of the load car- 
ried by the strong lithosphere itself as 
well as that carried by the immediately 
underlying weak magma. The ratio of the 
supporting forces contributed by the 
weak magma immediately below to that 
provided by the strong lithosphere is im- 
portant because it measures the approach 
of any given load to ideal isostatic ad- 
justment. Thus, if the weak substratum 
carries all the load, the ratio is unity and 
ideal isostasy prevails, whereas if the 
load is carried entirely by the lithosphere 
and distributed to the substratum over a 
large area, the ratio is small. A plot of the 
fraction of the total load carried by the 
asthenosphere immediately below is giv- 
en in figure 3 as a function of the semi- 
dimension of the block. It follows from 
this curve that, for a long block load to be 
even in approximate (96 per cent) iso- 
static equilibrium, its width must be not 
less than 330 km.—a value notably larger 
than the values frequently suggested in 
the literature of isostasy. 

FIBER STRESS AND MECHANICAL FAILURE 

WITHIN THE LITHOSPHERE 


The foregoing outline has presented 
curves showing how the downward de- 
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formation of the lithosphere under vari- 
pus types of loading can be determined, 
and one may, therefore, state the ap- 
proximate position of the neutral axis 
relative to the original surface of the 
earth at every point. Thus, by using the 
standard methods of engineering me- 
chanics, the fiber stress at every point in- 
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when subject to sufficient deformation. 
Thus, for example, in a study of the Cape 
Henry area (Gunn, 1944) it was possible 
to calculate with a fair degree of accuracy 
the approximate position where the 
break along the eastern coast of the 
United States was to be expected. This 
coincided remarkably well with the ob- 
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I'1G. 3.—The fraction of a long block load carried by the asthenosphere immediately below 


side the equivalent beam can be formally 
calculated. Actually, however, the graph- 
ical methods suggested above are not par- 
ticularly adapted to this type of calcula- 
tion, and analytical solutions are much 
more accurate. Because this type of 
analysis is beyond the scope of this sim- 
plified review, it is only remarked that, 
since the fiber stress within the litho- 
sphere can be calculated for every point 
as a function of its distortion, then, with 
the knowledge of the strength of rocks 
that constitute the lithosphere, one may 
immediately estimate where it will break 


served Fall Line and emphasizes the use- 
fulness of the techniques outlined above. 


COMPARISON OF VARIOUS TYPES OF 
LITHOSPHERIC EQUILIBRIUM 

The significance of the foregoing work, 
summarizing the mechanical behavior of 
a strong elastic lithosphere, is best under- 
stood by comparing the new assumptions 
with earlier ones. 

Three fundamental hypotheses may be 
made concerning the lithosphere. It may 
be assumed to be (a) perfectly rigid, 
(b) plastic, or (c) elastic. The implica- 















tions of conditions a and 6 are readily 
formulated in quantitative terms and 
have been widely employed in the calcu- 
lation of the Bouguer and the isostatic 
anomalies, while condition c is the as- 
sumption made in this communication. 
Comparison of the three assumptions is 
summarized in table 1. It should be noted 
that the principal term of the gravity 
anomaly for an elastic lithosphere re- 
duces to zero when dh = d,y, correspond- 
ing to the condition for ideal isostasy, 
and to 2rydh when the crust is perfectly 
rigid and y = o. Thus the elastic hy- 
pothesis leading to the principle of iso- 





278 ROSS GUNN 





gravity-anomaly data. Investigations out- 
lined in this paper show that gravity 
anomalies do not, in general, imply an 
unstable geological condition but indi- 
cate, rather, that part of the stabilizing 
system of forces arises from the strength 
of the lithosphere. The outlined prin- 
ciples are of the utmost importance in 
understanding the basic mechanisms re- 
sponsible for mountain-building. In fact, 
this investigation was an outgrowth of 
the two dilemmas of mountain-building 
that result from applying the classical 
isostatic theory. These two dilemmas 
were (a) the impossibility, according to 


TABLE 1 


COMPARISON OF VARIOUS TYPES OF LITHOSPHERIC EQUILIBRIUM 


| 
| . ess ° e 
Assumed Sufficient Approximation 
Character of 


Lithosphere 


| Resulting Character 
of Equilibrium 


for Actual Earth 
when Load Is: 


Rigid 


Plastic Uniform and larger than| Hydrostatic 
i 330km. 
Elastic Unrestricted Mechanical and 


hydrostatic 


* In these quantities, / is the height of an impressed load having a density d, y is the downward deformation of the 


Small and concentrated | Mechanical 


. _ Principal Term of 
Lithospheric Type of - wh: 
: wv ty 
Constant Anomaly 7 
Anomaly* 
Zero Bouguer | 2rydii 
Infinite Isostatic | Zero 


8.4X10%cm Isobaric 


(Gunn) 


2my(dh—d, y) 


eutral axis 


of the lithosphere, do is the density of the underlying magma, and y is the gravitational! constant. 


baric equilibrium, as here outlined, em- 
braces the others as special cases. 


REMARKS 


The foregoing outline of the basic prin- 
ciples of isobaric equilibrium emphasizes 
the limited applicability of the principle 
of isostasy and attempts to achieve a gen- 
eral understanding of the structure of the 
crust and its behavior under geological 
loadings of various types. The investiga- 
tion is more than a study of a special type 
of regional isostasy based on well-estab- 
lished mechanical principles because it 
permits one to forecast, for the first time, 
the magnitude and distribution of gravi- 
ty anomalies from structural data and, 
conversely, to deduce the structure from 





isostatic theory, of laying down prisms 
of sediments of adequate thickness and 
(6) accounting for the very rapid rise 
of geosynclines after an extremely long 
period of development and_ equilibri- 
um. A strong elastic crust having the 
properties outlined above is necessary 
to resolve these inconsistencies (Gunn, 
1937). 

The foregoing conclusions also modify 
in an important manner some of the 
basic ideas concerning the “roots of 
mountains.” If the conditions of ideal 
isostasy prevailed, actual mountainous 
“roots” would always be formed, and 
they would extend downward for a con- 
siderable distance. But, actually, because 
of the great strength of the lithosphere 
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and its elastic characteristics, an appreci- 
able share of the weight of most moun- 
tains of medium size is carried by the 
lithosphere and is spread out over a dis- 
tance approximately 250 km. each side 
of the load. Thus the root and the com- 
pensation is to a considerable extent not 
under the mountain itself but spread out 
to the sides. The resultant roots under 
the mountains consequently are relative- 
ly much shallower but much more exten- 
sive than one would estimate, were the 
region in isostatic equilibrium. 

The effect of a strong lithosphere in 
spreading the ultimate support for loads 
having a limited area over a considerably 
greater area emphasizes the artificiality 
of the isostatic hypothesis. In the reduc- 
tion of isostatic data, the assumption is 
made that the support or compensation 
isimmediately below the load. Although 
this assumption simplifies the computa- 
tions, it may be seen from figure 2 that it 
is not justified and is hardly descriptive 
of the actual geometrical configuration- 
It is evident, therefore, that gravity 


anomalies derived from isostatic reduc- 
tions are commonly unreliable and mis- 
leading except when applied to relatively 
large areas showing little deformation. 
The free-air gravity anomalies are much 
to be preferred in discussing geological 
structure because they more nearly rep- 
resent appropriately corrected observed 
data and are readily computed. Any 
proper topographical corrections must be 
made with full allowance for the fact that 
the support for any given earth feature 
is distributed over a considerable area. 
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It is recognized that the introduction 


of a strong elastic lithosphere between a 
given load and its ultimate hydrostatic 
support spreads the load and complicates 
the laws of crustal equilibrium. However, 
the principles of “isobaric equilibrium” 
well describe the actual mechanical bal- 
ance of the lithosphere and are unusually 
fruitful in describing various geological 
phenomena in terms of measurable gravi- 
tational variations. This extension of 
isostasy will, therefore, add notably to 
the interpretations of structural geology. 
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A REVIEW OF TECTONIC RELIEF IN AUSTRALIA! 


Cc. A. COTTON 
Victoria University College, Wellington, New Zealand 


ABSTRACT 


Siissmilch and others have maintained for many years that tectonic scarps in eastern Australia, mainly 
fault scarps, are prominent in New South Wales, and more recently they have described the eastern ranges of 
Queensland as made up of elongated fault-block features. Alternative explanations recognize a prevalence 
of multicycle relief forms in New South Wales and attribute the features of the Queensland ranges to subse- 


quent erosion on an ancient structure. 


In Victoria, South Australia, and parts of Western Australia, large block-faulted relief forms are recog- 
nized; but in the first-named state some scarps formerly regarded as tectonic have been more recently rele- 


gated to the fault-line scarp category. 


In Tasmania relief is dominated by fault scarps which are in a sense also structural escarpments; and 
survival of these in a single cycle from a Lower Miocene tectogenesis is claimed on the strength of paleo- 


botanical evidence. 


INTRODUCTION 

Lack of agreement as to what are and 
what are not fault scarps in Australia 
and the unsettled condition of opinion 
regarding tectonic relief in general in the 
eastern part of the continent point to a 
need for further investigation of some of 
the major features of Australian relief. 
Various methods of investigation of ma- 
jor relief have been employed, ranging 
from one which confines attention to 
underground structure to that based on 
recognition of the dislocation of former 
land surfaces tentatively restored for the 
purpose; but the features of few land- 
scape scarps themselves have been fully 
investigated or described. 

A theory of prevalent tectonic relief in 
eastern Australia was advanced long ago 
(Andrews, 1910), and it is now generally 
agreed (though possibly with some dis- 
sentients) that some of the major fea- 
tures had such an origin. With regard, 
however, to the extent to which tectonic 
relief prevails, there are distinct schools 
of thought. 


* Manuscript received August 16, 1948. 
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NEW SOUTH WALES 

Siissmilch (1910-1937) has maintained 
for many years that prominent tectonic 
scarps, mainly fault scarps, are numer- 
ous in New South Wales. On the other 
hand, however, some high scarps former- 
ly thought to be tectonic have been ex- 
plained as resulting entirely from long- 
continued selective erosion of a hetero- 
geneous terrain. It is therefore suspected 
that such an explanation might be found 
applicable to other scarps also, with the 
result that all scarps for which a tectonic 
origin has been suggested are regarded 
in some quarters with the strongest sus- 
picion. 

It is indeed desirable that ‘the sure 
methods of stratigraphy . . . and of de- 
tailed geological mapping” recommended 
by some Australian geologists should be 
applied to the problem, and all relevant 
geological data must be considered. It is 
a fact, however, that there is a consider- 
able lag in geomorphic description and, in 
the absence of this, geological methods of 
investigation that do not amount to ex- 
posure of faults and reliable inferences 
as to date of faulting lead in many 
cases to inconclusive results. Some work- 























ers, notably Siissmilch, have placed con- 
fidence in doubtful evidence of disloca- 
tion of the basalt sheets (of early Plio- 
cene age, according to Siissmilch [1937]) 
of subaerial origin, which are widespread 
in eastern Australia. Inherent in this 
view of large-scale landscape faulting is 
the conception of nearly perfect pene- 
planation over a very large area, making 
the Tertiary “great Australian,” or 
“Fastralian,” peneplain. 

Faith in this conception has been 
shaken, however, by some investigations 
(such as that of Voisey [1942]). These 
show (1) that the land surface on which 
basalt lies is, in some districts, one of 
rather strong relief and (2) that this re- 
lief was ancient when buried, being of 
a residual nature above the general level 
of the Tertiary peneplain in some ex- 
tensive areas characterized by the pres- 
ence of large bodies of resistant rock in 
the terrain. These areas include not only 
the New England Plateau of northern 
New South Wales, in which ancient relief 
has long been recognized by Andrews, 
but also the Southern Tablelands (Craft, 
1932, 19334, 6). Here, it is pointed out, 
basalt capping is found at very different 
altitudes, but this may not be interpreted 
as a result of dislocation of a once con- 
tinuous sheet covering a peneplain. “On 
the contrary they [the basalts of Monaro 
and Kiandra| are associated with sur- 
faces which bore the same relation to one 
another as they do now.” 

Throughout southern New South 
Wales and bordering the New England 
Plateau of New South Wales more or 
less prominent scarps have been diag- 
nosed as of tectonic origin by some au- 
thors. Undoubtedly, the early to middle 
Tertiary cycle of erosion was preceded in 
eastern Australia by a long period of 
earlier subaerial denudation. The region 
has been above the sea since late Paleo- 
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zoic times, and it seems that in the Ter- 
tiary cycle only the final touches were 
given to the eastern Australian pene- 
plain, where it became widely developed 
over areas underlain by weak rocks, be- 
ing there remade from a more ancient 
peneplain. It seems impossible to reject 
the belief that for some reason the tempo 
of denudation was very slow in this re- 
gion, thus allowing the survival of much 
relief in hard-rock districts. This was not 
mountainous relief, which is necessarily 
short-lived, but there were very ex- 
tensive residual reliefs on large bodies of 
resistant rock, which have been only very 
gradually subdued in successive cycles 
during eons of erosion. Such an explana- 
tion of the landforms of the New Eng- 
land district has been generally accepted 
since it was offered many years ago by 
Andrews; but its application to the ex- 
tensive plateaus of southern New South 
Wales has been made gradually and cau- 
tiously by various workers, notably 
Browne (1929, 1944) and Craft (1932, 
1933a, b). More recently, Maze (1944) 
also has recognized down-stepping pene- 
plain benches west of the Blue Moun- 
tains. 

One important conclusion reached by 
Craft is that much of the relief of eastern 
New South Wales has survived a denu- 
dation following post-Mesozoic and even 
pre-Mesozoic updomings. “The high- 
lands thus formed have been subjected 
to continuous invasion,” however, “and 
incomplete peneplains have been eroded, 

. . which lead to higher cores of gently 
concave aspect.” 

Notwithstanding the discovery that 
large parts of the area of New South 
Wales were diversified by considerable 
relief prior to the upheavals which are 
attributed to late Tertiary crustal un- 
rest—the almost-orogeny of Andrews’ 
Kosciusko epoch—the formation of tec- 
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tonic surface features as a result of di- 
verse upheavals of this date cannot be 
doubted. However, the magnitude and, 
in particular, the abruptness of the tec- 
tonic relief have been perhaps over- 
estimated by Andrews (1910, 1934), Tay- 
lor (1910, 1911, 1918, 1923), and Siiss- 
milch (1910, 1911, 1914), who have been 
convinced of the tectonic origin of most 
of the large landscape features. Craft, on 
the other hand, relegates tectonic relief 
to comparative unimportance. He rec- 
ognizes no conspicuous scarps of this 
origin, though he allows for differential 
upheaval in some upwarped areas. 
“Viewed as a whole,”’ he writes, “‘the re- 
gion has been subjected to a gentle post- 
basaltic uplift of 2000 feet on the median 
line of the highlands, with relatively de- 
pressed portions.” 

That the true explanation lies between 
the extreme views that have been ad- 
vanced—that certain scarps, as well as 
extensive broad domings, have been de- 
veloped and still remain features of the 
landscape—is shown by a description re- 
cently published of a stepped or splin- 
tered fault scarp descending eastward 
from the highest (or Kosciusko) plateau 
remnant in southern New South Wales 
(Browne, Dulhunty, and Maze, 1944). 
There is also the unchallenged authority 
of David (1908) for the occurrence of 
fault scarps in that locality. It appears, 
however, that some major scarps which 
have been diagnosed as tectonic in south- 
ern New South Wales must now be con- 
sidered the unfinished work of secular 
subaerial denudation on a heterogeneous 
terrain of rather ancient rocks. 

Perhaps the straight Cullarin scarp, or 
rather the diagnosis of this scarp by Tay- 
lor (1907) as a fault scarp, will stand the 
test of a long-delayed full discussion 
which it deserves. This scarp, about 40 
miles long but of moderate height (about 
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500 feet), is the western boundary of the 
remarkable Lake George Basin of in- 
terior drainage on the broad crest (near 
Canberra) of the main divide between 
eastern and Murray Basin drainage to 
the ocean. This scarp defines the line of 
one of “two great faults” postulated by 
Taylor (1910) to account for prominent 
lineaments of the landscape of the Feder- 
al Territory at Canberra, though on the 
second of these lines, that followed by the 
Murrumbidgee River for a long distance, 
faulting either modern or ancient seems 
to be a very remote possibility according 
to the geological work of Browne (1944). 
Discovery of the Cullarin scarp in- 
spired Taylor to produce a geomorphic 
description more detailed perhaps than 
any other to be found in Australian lit- 
erature; and probably this is one of the 
best authenticated of Australian fault 
scarps. After describing the sprawling 
spurs that descend between broad valleys 
into the alluviated Lake George plain from 
the range east of it, Taylor continues: 


Reverting to the western shore, we are struck 
by the absolute dissimilarity. Standing at the 
level of the lake we seem to be confronted by a 
great wall extending northwards for over 20 
miles from the Molonglo Plain. No broad valley 
breaks its continuity. Indeed . . . it seems un- 
broken by any definite gap, while for a large 
part it presents a steep face 500 feet high to the 
lake. From the eastern shore, however, one is 
able to see a well-defined gap about halfway 
along the western shore, where the old Southern 
road crossed the Cullarin Range... Geary’s 
Gap... . Less than two miles south of Geary’s 
Gap a stream, Grove Creek, rushes down to the 
lake. Contrast its course with that of Taylor’s 
Creek, on the eastern shore. The latter stream 
flows through a broad valley a mile or two wide, 
and shows the even grade of a long established 
river. The Grove Creek is barely a mile long, 
yet descends nearly 300 feet. Its course is inter- 
rupted by falls 25 feet high, and finally it 
emerges from a gorge or miniature canyon with 
steep sides 200 feet high. Evidently it is a stream 
which has barely reached the youthful stage. 
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He finds also that the Molonglo River 
enters the scarp farther south and crosses 
an upheaved block west of it by way of 
an antecedent gorge. Geary’s Gap may 
be a similar gorge which has become a 
wind gap as a result of defeat of the 
stream which had followed a course 
through it until it was ponded in the 
Lake George Basin by acceleration of the 
upheaval. Similar features furnish good 
evidence of block faulting in other regions. 

Acceptance of Taylor’s description 
and diagnosis of the Cullarin scarp as a 
fresh, or very youthful, fault scarp per- 
haps distinguishes it from other Aus- 
tralian tectonic scarps, most of which are 
rather the worse for wear. Either the dis- 
location on the Cullarin line took place 
more recently than most of the others, or 
circumstances have combined to slow 
down the erosional degradation of this 
scarp. The latter may well be the cor- 
rect explanation, for no erosion by a 
through river seems to have taken place 
in the Lake George Basin which could 
affect the form of the scarp. The en- 
deiric drainage of the basin is unique in 
eastern New South Wales, and the pecu- 
liarity is attributable to a combination of 
climatic conditions with accidental dis- 
position of blocks and the pattern of the 
land surface consequent on deformation. 
Remoteness from rivers and immunity 
from dissection except by short, conse- 
quent streams may allow a low scarp to 
retain its primitive form almost indefi- 
nitely. Arid erosion, moreover, favors the 
perpetuation of steep slopes, though 
these may retreat without appreciable 
change of form. 

It is of some interest to turn from the 
Cullarin scarp to examine reasons which 
led both Taylor (1910) and Siissmilch 
(1910) to regard as tectonic another long 
scarp in the Canberra district, that which 
descends to the meridional reach of the 
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Murrumbidgee River from the highland 
to the west of it, although Browne’s geo- 
logical work in the district does not sup- 
port the diagnosis and Craft (1933a) re- 
gards the scarp as a feature surviving 
from Tertiary and pre-Tertiary relief. 
For identification of a fault scarp on this 
Murrumbidgee line Siissmilch has relied 
on a postulated dislocation of a formerly 
continuous peneplain, which he claims 
was very widespread and _ low-lying 
throughout eastern Australia in mid- 
Tertiary times. The dislocation, assumed 
to take place at the time of the generally 
recognized late Tertiary or end-Tertiary 
upheavals, would make a scarp 2,500 feet 
high. Taylor’s diagnosis of faulting is also 
based in part on an assumption that it is 
necessary to postulate such dislocation of 
the landscape to account for a rectilinear 
scarp along the left bank of the Murrum- 
bidgee River. However, it must be ad- 
mitted that, apart from a regional break 
in the topography, such details of land- 
scape as are shown by Taylor’s maps and 
the photograph of his relief map of the 
Federal Capital Territory do not suggest 
a fault in the same vivid way as on the 
Cullarin line. 

In explanation of the scarped valley- 
side profile, the fact is significant that 
streams tributary to the Murrumbidgee 
River flow, probably because of adjust- 
ment to some ancient structure, parallel 
and close to the main river for long dis- 
tances and thus rob it of the accustomed 
quota of valley-side drainage. Where any 
linear cliff is deprived in this way of the 
water necessary to feed even rivulets, its 
dissection is necessarily greatly retarded. 
A valley side thus assumes and long re- 
tains a wall-like form resembling that of 
a young fault scarp. 


Some generalized accounts of the 
physical features of New South Wales 
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imply that, apart from those which may 
be attributed to marine cliffing, there are 
distinct scarps along the coastal margins 
of high plateau areas which must be tec- 
tonic. According to Craft (19330), how- 
ever, who bases his opinion on a meas- 
ured gentleness of regional slopes, it is 
unnecessary to postulate that there was 
initially any abrupt seaward descent or 
variation from very gentle regional 
monoclinal warping, which has, more- 
over, merely steepened a regional slope 
inherited from distant prebasalt times. 
“The plateau edges are primarily ero- 
sional features,” and “‘when the rocks of 
the landscape are more homogeneous the 
streams take on a symmetrical charac- 
ter” on the seaward marginal slopes. As 
Craft points out, however, “portions of 
relatively undissected plateau come near 
the coast at intervals throughout its 
length.” At such places it is possible that 
the initial seaward descent has merited 
description as tectonic scarps. This is 
suggested by some statements of An- 
drews (1914), such as: “At Illawarra and 
still farther southwards the plateau ad- 
vances almost to the shoreline’; but, 
actually, where the seaward face of the 
inland plateau becomes a high and 
bold scarp, this is a structural escarp- 
ment. That of the Illawarra coast, which 
Andrews cites, has resistant Triassic 
sandstone as cap rock, and the coastal 
profile of which it forms a part may have 
been formed by the processes of erosion 
from a gently inclined monoclinal slope 
just as well as from a fault scarp. 

A very gentle seaward tilting of a sur- 
face already sloping seaward at the end 
of the Tertiary era has been sufficient, 
according to Craft, for the rejuvenation 
of the valleys of coastal rivers in both 
northern and southern New South Wales, 
thus causing formation of falls and the 
excavation of gorges, commonly 1,500 
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feet deep and in some cases deeper. In 
northern New South Wales (New Eng. 
land) the presence of a basalt capping on 
the plateau at the headwaters of some of 
the coastal rivers has checked their head- 
ward erosion and steepened their gradi- 
ents at a zone of escarpments; and it is 
questionable whether any of the basalt 
escarpments at the margin of the high 
plateau, however imposing, are derived 
from true tectonic scarps. An example of 
such is the escarpment at Guy Iawkes, 
overlooking deeply cut headwater val- 
leys of the Bellingen and Nambucca 
rivers, “the view from which is probably 
unexcelled in New South Wales, or even 
in Australia” (Andrews, 1914). Accord- 
ing to Andrews, this feature may “be ex- 
plained either by erosive or by tectonic 
forces.”’ He adds: “Similarly for the great 
eastern scarp overlooking the South Arm 
of the Clarence.” Though (at another 
place) he views the marginal feature 
emphasized by these escarpments as 
“strongly suggestive of an important 
fold or fault,” it appears that no posi- 
tive evidence of the presence of a tec- 
tonic scarp has come to light. It is sig- 
nificant, on the other hand, that in a 
part of the seaward slope from New 
England which Craft selects as typical, 
the basin of the Macleay River, he finds 
no strong scarp in the present-day form 
of the surface, not even an escarpment 
at the margin of the basalt capping. 
Siissmilch has not attempted to rec- 
ognize tectonic features in either the 
seaward or the inland margins o/ the 
plateaus, merely remarking (1937) that 
“since their uplift the tablelands have 
suffered considerable dissection, particu- 
larly along their eastern and western 
margins.’ Benson (1917) has noted, how- 
ever, that, although “the plateaux of 
New England sink to the west and pass 
out into the plains,” one high salient or 
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buttress forming the Nandewar Moun- 
tains is bounded on the west (on the au- 
thority of Jensen [1907]) by a tectonic 
scarp. 


One high and prominent scarp in New 
South Wales has been given some recent 
attention. This scarp, on a line recog- 
nized even by Craft as a fault (‘‘Hunter 
boundary fault’), forms part of the 
southern margin of the New England 
plateaus. There is here a steep descent 
to the low-lying Hunter River Valley, a 
tectonic sag, much enlarged by erosion, 
which is continued westward in the 
Goulburn Valley and the major geo- 
graphical feature called by Taylor the 
Cassilis ‘‘geocol.” This sag separates the 
northern and middle members of the 
three broad upswellings comprising the 
New South Wales portion of the eastern 
“cordillera,”” or Main Divide, of Aus- 
tralia, which are described by Craft 
(19330) as “three low domes,” although, 
according to him, the Main Divide as a 
whole is “of essentially non-tectonic ori- 
gin.”’ Siissmilch (1940) has described 
parts of this scarped descent from the 
Barrington salient of the New England 
plateaus westward and southward to the 
Goulburn-Hunter Valley, where fault 
scarps had been earlier diagnosed by 
Andrews (1910). 

According to the description given by 
Siissmilch, these slopes contrast strongly 
not only with the subdued landforms of 
an extensive bench at their foot, which 
is regarded as part of a maturely dissect- 
ed peneplain, but also with the small re- 
lief of a plateau above, though exact cor- 
relation in the plateau and bench of dis- 
located parts of a formerly continuous 
land surface has not been made. About 
60 square miles of highland plateau sur- 
vive near the western edge of the Bar- 
rington salient at an altitude of about 


5,000 feet. The plateau consists super- 
ficially of basalt, which is found to be a 
capping 500-800 feet thick over pene- 
planed granite and other rocks. West of 
a marginal scarp, presumably tectonic, 
the base of the basalt is 1,200-1,500 feet 
lower, indicating a lowering in that di- 
rection along a “‘post-basaltic fault or a 
monoclinal fold.” Topographically the 
scarp form is not very clearly marked, 
but rivers which presumably originated 
in courses consequent on its slope are 
superposed in gorges across strike ridges 
of Devonian and Carboniferous bedrock. 
Farther west at a lower level a second 
west-facing scarp of similar height is re- 
garded as another fault scarp in a 
stepped descent. This one is related to 
the known Wingen fault, which brings 
Permian and Triassic strata on the west 
side against the folded Carboniferous. 
The displacement here, as measured on 
dislocated basalt capping, is 550 feet; 
but the scarp on this line remains a 
strongly marked topographic feature 
only where parts of it take on the charac- 
ter of a structural escarpment, owing to 
the capping of basalt. The well-dissected 
condition of these scarps indicates that 
the Kosciusko uplifts to which they are 
attributed were not very recent. 

Along the southern margin of the Bar- 
rington salient, or southwesterly pro- 
jecting buttress of the New England 
plateaus, there is a scarp feature which 
remains more strongly marked topo- 
graphically and is described by Siissmilch 
as “‘a magnificent scarp trending almost 
due east and west and giving a drop of 
more than 3000 feet from the higher to 
the lower tableland.” It “displays all the 
characteristic features [not, however, 
enumerated] of a true fault scarp, and 
was diagnosed as such by E. C. Andrews 
in 1910.” The fact that the scarp ap- 
pears to be entirely in basalt is explained 
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by Siissmilch as due to a misleading 
effect produced by the presence of much 
basalt talus on it, which entirely masks a 
bedrock consisting presumably of ancient 
rocks. No alternative hypothesis of ori- 
gin of a structural escarpment (capped 
by the basalt) at this line without fault- 
ing seems to offer a satisfactory explana- 
tion of the feature. 


The middle member of the eastern 
Australian group of high plateaus, the 
Blue Mountains—or Blue Mountain 
Plateau or Tableland, as it is better 
termed—presents the least ambiguous 
tectonic feature in the whole region, an 
enormous monoclinal flexure of the land 
surface, which descends from high pla- 
teau summit levels (in places about 4,000 
feet) nearly to sea level a short distance 
west of Sydney. This monocline and its 
local steepening northward into a con- 
spicuous scarp were noted by David 
(1g02) and Andrews (1903). The steeper 
scarp portion is described by Taylor 
(1918, 1923) as monoclinal and by others 
as a fault scarp. Such features, when they 
have undergone dissection, are not readi- 
ly differentiated. 

A brief general account of the land- 
scape of the Blue Mountains asa whole has 
been given by Taylor (1918) as follows: 

The Blue Mountains consist of an ancient 
peneplain, which isa relicof the much greaterland 
surface of early Tertiary times. The eastern 
portion has sunk beneath the sea, while the cen- 
tral portion has been buckled up to the extent 
of several thousand feet....The western 
. merge gradually into the Western 
Plains, and the abrupt scarps of the eastern 
side are not present there. 


slopes . . 


Taylor considers the peneplain to have 
been similar, prior to its upheaval, to 
that which survives far down the west- 
ern slope of the arch, where the “Cobar 
peneplain” is surmounted by scattered 
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low monadnocks with maximum relief 
of 500 feet. 

According to Craft (19330) also: 

A general relief of 500 feet may be ascribed 
to the Blue Mountains portion of the pre-basalt 
surface, and a possible altitude of 1000 feet 
above sea-level, with such ancient residuals as 
Mount Lambie rising almost to 2000 feet. ... 
The surface has been deformed. . . . There has 
been tilting or faulting towards the low coun- 
3 Ae 
At the time of the basaltic extrusion the 
Central Plateau or Blue Mountains formed a 
low region between the higher masses. Its 
general elevation probably did not exceed 
1000 feet, with a maximum of 2000 feet, whilst 
the higher mass of New England rose to 3500 
feet, and the Monaro-Kosciusko region to 
4500 or 5000 feet as a maximum. There appears 
to have been little subsequent deformation of 
the higher regions, but the Central Plateau 
may have been warped considerably. 

... The most complete peneplanation record- 
ed in the history of existing landforms is in the 
region of the Blue Mountains, where monad- 
nocks are comparatively insignificant points. 


The Central Tableland as defined by 
Siissmilch (1914) includes, besides the 
Blue Mountains, other considerable areas 
from 2,000 to 3,000 feet in altitude, 
which he has regarded as distinct blocks 
bounded and separated from one another 
by tectonic scarps. These scarps he has 
not described, however, except as “abrupt 
changes in elevation,” which are likened 
to the steeper part of the Blue Moun- 
tains monocline. The southern boundary 
of the Blue Mountains unit he has termed 
a “magnificent fault escarpment,” but 
the correctness of this identification of a 
fault scarp is perhaps questionable; for 
Craft (19336) mentions no appearance 
of faulting of the land surface at or south 
of this line and, indeed, describes the 
southerly descent from the Blue Moun- 
tains as a warped surface. 

From the foregoing remarks regarding 
Kosciusko upheavals and dislocations 
and post-Kosciusko dissection in New 
South Wales, it will be gathered that the 
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tectonic features of eastern Australia are 
of rather ancient origin and are for the 
most part well dissected. A post-Pliocene 
date seems to be out of the question for 
these upheavals; but there is no reason, 
on the other hand, for assuming that the 
origin of the surviving tectonic features 
may be dated earlier than Pliocene, not- 
withstanding the apparently necessary 
postulate of a slow tempo of denudation 
in this region. 


QUEENSLAND 


A recently published summary by 
Siissmilch (1938) of geomorphic evidence 
strongly favors recognition of numerous 
tectonic scarps, and fault scarps in par- 
ticular, in the eastern coastal districts of 
Queensland. He writes: 

In most parts of eastern Queensland very 
little geological mapping has been carried out; 
consequently definite geological proof of fault- 
ing is not available for many of the supposed 
fault scarps. These particular scarps are, how- 
ever, so similar in their physiographical features 
to those for which geological proof is available 
that it is quite reasonable to assume that all of 
them have had a similar origin. 


He lists several categories of scarps which 
he regards as fault scarps and continues: 

It may be true that some of the scarps are 
due to warping rather than faulting, but this 
would not alter the fact that they are tectonic 
scarps and not erosional scarps. If this con- 
clusion is correct, then the coast ranges and 
chains of continental islands are fault blocks, 
the corridors, both land and marine, are rift 
valleys, while the coastal plains and broad 
lowland regions are still-stand areas—i.e., 
parts of the late-Tertiary land surface which 
lagged behind during the general uplift. 


Prior to the publication of this con- 
tribution by Siissmilch and that of Stan- 
ley (1927), who has admirably described 
aportion of the Queensland coastal area, 
generally similar conclusions had been 
reached, though stated in broader terms, 
by Andrews (1910) and by David (1914). 
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Some Queensland geologists, however, 
notably Marks (1924) and Bryan (1928, 
1930), have scouted the idea of the exist- 
ence of any tectonic relief, mainly on the 
ground that the trend of elongation of 
coastal features, ridges, and corridors 
finds a sufficient explanation in the grain 
of the folded Paleozoic terrain. Thus 
Bryan in a broad generalization attrib- 
utes all elongation and alignment of 
features that lie more or less parallel to 
the coast to dominance of the north- 
northwesterly strike of the Paleozoic 
structure. The river courses, together 
with corridors which other observers 
point out are not and probably never 
have been occupied by through rivers, are 
explained as “subsequent’’—i.e., of an 
origin attributable solely to selective 
erosion; and what is true of valleys must 
be true likewise of interfluvial highland 
ridges. 

One of the main arguments employed 
by Siissmilch in Queensland, as in New 
South Wales, is based on recognition of 
dislocated parts of a formerly continuous 
peneplain under the offshore shelf, on 
the floors of low-lying coastal strips and 
corridors of small relief, and on the 
crests of highland blocks in the coastal 
ranges. No plateau (peneplain) remnants 
are recognized on the coastal ranges 
by Bryan, however, although Andrews 
(1910), Danes (1911), and others have 
there thought to identify parts of a well- 
worn surface which has commonly been 
regarded as a peneplain of Tertiary de- 
velopment, presumably continuous with 
that of New South Wales. On such a 
surface in Queensland, according to 
Stanley (1927), “dissection had reached 
maturity, ... but peneplanation was not 
nearly complete.” 

Regarded as an important “‘datum”’ 
by Bryan (1939) are red earths and lat- 
erite of a series of which the origin is 
attributed to a moist Pliocene epoch im- 
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mediately preceding the Kosciusko up- 
lifts of New South Wales. Though found 
at many places in southeastern Queens- 
land on a level surface, and on basalt cov- 
ering such a surface, this soil formation 
has been reported only at low altitudes 
on straths of an imperfect partial pene- 
plain which has failed to truncate the 
more resistant rocks of the region. The 
discovery of the red-earth datum has not 
yet afforded a key to the geomorphic his- 
tory of Queensland; but, if the date at- 
tributed to it is acceptable, it seems to 
push back the date of any recognizable 
higher-level surface of advanced denuda- 
tion there may be to a remote past. This 
is in conformity with the theory of slow 
tempo of Australian geomorphic de- 
velopment. Yet the uncertainty of the 
date of origin of such a surface, if it is 
real, does not affect the value of infer- 
ences based on its presence at apparently 
dislocated levels. 

Among the most striking features 
known as major landscape forms in 
Queensland are the so-called ‘‘corridors”’ 
(an empirical term first employed by 
Jardine [1925] for them). According to 
Bryan’s interpretation, they are no more 
than valleys formed at some time by sub- 
sequent rivers; but, according to the tec- 
tonic interpretation, they are grabens. 
They are elongated parallel to the coast, 
with either slightly emergent or sub- 
merged or buried floors. Intercorridor 
horsts, corridors, belts of outlying 
islands, and a hypothetical bedrock 
bench generally believed to underlie the 
very broad continental shelf and Great 
Barrier Reef at a shallow depth are 
elongate features aligned parallel to a 
coast whose initial form, according to the 
tectonic explanation, was caused by 
rather recent earth movements. 

The evidence of marginal and block 
faulting is found for the most part in 
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steep regional slopes that separate “table. 
lands,” which surmount the highest 
ridges in Queensland (all in the coastal 
belt), either from benches which are ap. 
parently step-faulted or from nearly level 
corridor floors, ‘coastal plains,”’ and the 
undersea surface of the continental shelf, 
which is regarded as down-faulted. Siiss- 
milch (1938) has drawn attention to 
these by sketching generalized profiles of 
the landscape, which are designed also to 
show plateau relief on the ranges. 

Beyond this, the geomorphic evidence 
as to the nature of the scarps is only 
vaguely known. Proponents of the theory 
of tectonic relief have not, as a rule, 
given much attention to description and 
discussion of the features or to the stage 
of dissection of the fault scarps which 
they have diagnosed. That this dissection 
is fully mature may be inferred from such 
statements as “both groups of scarps 
[there are scarps transverse as well as 
parallel to the longitudinal topographic 
elongation or trend] display a similar 
amount of dissection” (Siissmilch, 1938). 
A similar inference may be drawn from 
the statement that ‘differential erosion 
has played an important part in the dis- 
section of the tablelands and horsts since 
their uplift,” for this allows time for, and 
implies, advanced dissection of the scarps 
that bound these features. If these tec- 
tonic scarps were still young, there could 
be but little dissection as yet of upland 
surfaces behind them. Yet differential 
erosion is reported to have opened valleys 
right through some horsts along belts of 
weaker formations which separate granite 
massifs. 

In addition to the longitudinal scarps 
attributed to faults bounding horsts, 
grabens, etc., some are mentioned which 
are “transverse to these trends.”’ As these 
latter are transverse also to the struc- 
tural grain of the terrain, they are per- 
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haps more plausibly related to faults, al- 
though whether these be ancient or rela- 
tively recent dislocations is another mat- 
ter. It cannot be said with certainty as yet 
that differentiation has been clearly 
made in Queensland between fault scarps 
and fault-line scarps; but the balance of 
probability seems in favor of the tectonic 
origin that has been confidently ascribed 
by various authors to scarps which ap- 
pear to bound horsts and grabens in the 
coastal belt. 

Prior to the publication of Siissmilch’s 
generalization, Stanley (1927) vigorously 
advocated the tectonic theory, which he 
applied to a typical portion of the coastal 
belt, describing it as the “flexed and 
faulted-down margin of the Pliocene con- 
tinent.”’ Describing Whitsunday Island, 
one of a group rising from the submerged 
shelf which David (1911) long ago re- 
ferred to as “almost certainly horsts 
among a network of faults,” Stanley 
wrote: 

Observe that there is here no question of 
differential erosion. The physiographic units 
are constituted of many different types of rocks, 
and they cut across the geological boundaries. 
In the island festoons both granites and Paleo- 
zoic volcanics figure equally, and the same is 
true of the coastal ranges. Important as rock 
structure and relative resistance to weathering 
may be in determining topographic forms, it 
does not seem feasible to consider that it has 
acted as the broad control which is required by 
the systematic arrangement of the units. 


Some Queensland rivers are very 
strongly rejuvenated in their lower 
courses, and this fact has often been 
cited as proof of the presence of either a 
coastal fault or a,strong flexure of the 
surface. Taylor (1911) and Dane§ (1912) 
have each described the Barron River 
and its falls as an example. This river, 
Taylor wrote, 
leaves the highlands . . . by falls 600 feet high. 
This gorge cannot be due to ordinary uplift, for 





the coral reefs in the immediate vicinity show 
practically no sign of elevation. But it is quite 
probable that the conditions which obtain in 
the south-east of the continent—namely, blocks 
of the earth’s crust bounded by north-south 
faults—have also determined the topography 
here. We may imagine (as Dr. Dane’ sug- 
gested) a Tertiary extension to the east for 
several hundred miles with a normal consequent 
drainage. The greater portion of this is faulted 
and subsides gradually, and is covered by the 
coral reefs of the Great Barrier. Only the head- 
waters of these rivers remain, and they are now 
in a much more favourable position than before 
the subsidence as regards rapid erosion... . 
Faulting occurs. . .. The river is betrunked in 
a very drastic fashion. The short coastal Barron 
has now to drop 1ooo feet in ten miles instead 
of in several hundred miles, and hence cuts out 
a gorge very rapidly. 


VICTORIA 


In eastern Victoria (the Eastern High- 
lands) the high-plateau country of the 
Kosciusko district of southern New 
South Wales has its southerly continua- 
tion. Nothing is definitely known in this 
region of any tectonic scarps of Kos- 
ciusko or earlier upheaval, although 
Hills (1934) notes that differential move- 
ments have undoubtedly taken place. 
Conspicuous scarps, however, both here 
and in the Grampians (Western High- 
lands of Victoria) have been shown to be 
due to selective erosion working on a 
heterogeneous terrain (Hills, 1936). In 
the southern part of Victoria, on the 
other hand, there is ample evidence of 
differential movements which may have 
been contemporaneous with Kosciusko 
uplifts, and fault scarps are recognized 
(Hills, 1940). 

According to the findings of Hills 
(1934, 1940), the geological history of 
Victoria is different in some important 
respects from that generally accepted for 
New South Wales. The most ancient land 
surface in Victoria is preserved in some 
small residuals of high dissected plateau 











on areas of very resistant terrain; and 
this oldest surviving peneplain is at- 
tributed to a cycle interrupted before 
the end of the Cretaceous period. A later, 
lower-level peneplain was extensively 
developed at the expense of the dwindling 
Cretaceous surface prior to the extrusion 
of extensive cappings of basalt, pre- 
sumably in the Oligocene period. The 
Oligocene land surface—a_peneplain 
with not inconsiderable monadnocks still 
standing above it—is known largely as 
a fossil surface because of its preserva- 
tion under the basalt. Following this, 
middle Tertiary erosion planed off con- 
siderable areas which were later sub- 
merged in marginal districts and buried 
under marine strata; and it is subsequent 
to this (in the Pliocene period) that the 
upheavals are dated which have made 
the more prominent features of Victorian 
relief (the Eastern Highlands and the 
Central Highlands). The major relief 
forms which resulted from these middle 
to late Pliocene movements still survive; 
but the movements were, according to 
Hills, of a pre-Kosciusko date. In Gipps- 
land (southeastern Victoria) the land 
surface was dislocated by movements of 
probably later, i.e., very late Pliocene, 
date (Kosciusko?). 

The lack of agreement between events 
and their sequence in New South Wales 
and Victoria is perhaps only apparent, for 
it is possible that dates attributed to 
major earth movements in Victoria will 
be found to be applicable over a wider 
region. The date assigned by Hills to 
Kosciusko uplifts (in the most restricted 
sense) is late Pliocene to early Pleisto- 
cene, at which time he believes that 
movement was renewed in the Central 
Highlands, while at the same time new 
ranges made their appearance in south- 
ernmost Victoria. There, especially in the 
South Gippsland highlands and in the 
Otway Ranges farther west, a tectonic 
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relief of very late origin is recognized, 
One of the best-marked scarps in South 
Gippsland is the maturely dissected Yar- 
ragon scarp, which faces north, overlook- 
ing an embayment of the East Gippsland 
Plains (a tectonic basin) ; but, in addition 
to this, several others are mapped by 
Hills along the northern and southern 
margins of the upheaved area. 

Hills has attempted to distinguish 
with care between fault-line scarps (which 
are generally two-cycle features) and 
true fault scarps. Many supposed fault 
scarps in Victoria have been relegated to 
the former category; but (he states) “the 
South Gippsland highlands offer an ex- 
cellent example of block mountains 
formed by Cainozoic faulting, and they 
constitute a well-defined physiographic 
region.”” At another place: “In south- 
central Victoria the district between the 
Rowsley fault and Selwyn’s fault has 
been depressed, forming the Port Phillip 
sunkland. ...’’ One of these faults, the 
Rowsley, regarded as of Pleistocene age, 
is marked in the southern part by a well- 
defined scarp, which Hills (1940) has se- 
lected as an example of a fault scarp dis- 
sected to maturity. It forms the eastern 
boundary (with 800-foot relief) of the 
Brisbane Range, a back-tilted block. l’en- 
ner (1918) in his description of this scarp 
has drawn attention to the fact that its 
topographic form complies with Davis’ 
ideal fault scarp of the Wasatch variety 
in that it truncates indifferently the 
structures in an upheaved block and has 
a bahada in front of it. It is dissected by 
a number of short streams which are 
presumably consequent. ‘‘Deep gorges 
score the face of the high resistant scarp, 
and the streams continue on the low 
block in wandering shallow channels 
across the alluvial apron” (Ienner, 
1918). 

Hills has also figured the profile of the 
Otway Ranges uplift, farther to the 
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southwest, as a tectonic feature limited 
by post-Pliocene monoclinal flexures 
steeper on the southeast (coastal) flank. 





TASMANIA 


Fault scarps have been pointed out by 
Carey (1947) in the Launceston district 
of Tasmania and are interpreted by him 
as ancient features, though in a good 
state of preservation. In his account of 
this district, which is accompanied by 
aset of comprehensive block diagrams in 
sequence showing paleogeography, strong 
tectonic relief is described and shown as 
having its origin in the Lower Miocene 
period, and no relief of any kind in- 
herited from earlier times or developed 
later by earth movements is recognized. 
No indication has been found of even 
minor renewal of faulting subsequent to 
the Miocene paroxysm; and only very 
small movements of base-level unac- 
companied by appreciable warping are 
believed to have taken place in later 
times. A very simple picture of tectonic 
and erosional history is thus presented. 

The predeformational surface, dis- 
located in Lower Miocene times, was, ac- 
cording to this interpretation, a land 
surface in an advanced stage of senility 
after a resting period of ninety million 
years. “This peneplain reached a high 
degree of perfection. No protuberances 
suggestive of old monadnocks . . . have 
been recognized.” The senile surface was 
crusted with laterite, especially where it 
crosses a very extensive sheet of Jurassic 
dolerite intrusive between Lower and 
Upper Permian strata. 

The Miocene dislocation of this sur- 
face developed gentle folds broken by 
strong strike faults, so that the resulting 
relief was expressed mainly as a number 
of elongated, subparallel fault blocks dis- 
posed as steps on the flanks of a few brok- 
en anticlines. 

The strong relief of early Miocene ori- 
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gin is depicted as surviving very exten- 
sively into the modern landscape. Some 
fault scarps which had rock of the same 
kind on each side, even granite, are 
shown as entirely obliterated, naturally 
enough, by thirty million years of ero- 
sion; but on the sites of ancient com- 
pound horst blocks with originally thick 
dolerite capping there are still ranges; 
and the scarps that bound these, though 
they might, if an alternative multicycle 
hypothesis were available to account for 
the relief, conform at least in some parts 
to the definition of fault-line scarps, must 
be interpreted as fault scarps which still 
survive with much of their initial relief, 
if the postulate of only one cycle of ero- 
sion can be maintained. 

Until some other interpretation of the 
geological and geomorphic history is 
offered, the story told in Carey’s dia- 
grams of the extremely slow wastage of 
escarpments bounding dolerite caps on 
some of the ranges may be tentatively 
accepted as a working hypothesis; and 
these are shown by him as true fault 
scarps. 

Assuming that no question can be 
raised as to the date of faulting, it is diffi- 
cult, however, to accept unconditionally 
the one-cycle explanation of this land- 
scape. For one thing, thirty million years 
have elapsed; and, for another, the fi- 
delity to nature of certain parts of the 
land surface as depicted, or else its in- 
ferred history, may pardonably be doubt- 
ed. Though some episodes of the history 
of the district have not been written, 
the landscape has been so graphically 
delineated in the diagrams (fig. 1) that 
the story thus unfolded is none the less 
open to criticism. 

At the foot of Ben Nevis, a broad 
bench that has apparently been leveled 
by Miocene and post-Miocene erosion 
and is cut across granite at about 3,000 
feet above base-level remains unac- 
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counted for. It is at this high level of 
unexplained planation that denudation 
has obliterated a high fault scarp of 
granite. In contrast with this scarp 
degradation, in the southwestern part of 
the rectangular area figured, another 
fault scarp, which bounds the large pre- 
Miocene peneplain remnant forming the 
Great Western Tiers and is 4,000 feet 
high, retains its relief and primitive 
steepness, though exposed for an equal 
time to erosion, and has been worn back 
but a short distance. It is difficult to 
believe in the secular survival here and 
in Ben Lomond (5,000 feet high) of 
mesas of strong relief consisting of rem- 
nants of a dolerite sheet, which, though 
very thick, is weakened by its overlying 
the Permian mudstone at a high altitude. 
It would seem that the explanation given 
of these Tasmanian scarps must be con- 
sidered as only tentative, or else that the 
date of the block faulting may be later 
than has been supposed. 

The latter suggestion finds support in 
an account given by Lewis in 1939 of 
studies then in progress, which indi- 
cated his belief that the block move- 
ments responsible for producing the ma- 
jor features were of late Pliocene date, 
until which time the Tertiary peneplana- 
tion had continued, and that even in the 
Pleistocene Tasmania was still being 
broken into blocks along the lines which 
are now followed by consequent rivers. 
However, Lewis veered later to the 
opinion that the main landscape faulting 
was early Miocene, from which date, he 
wrote (1946), “a normal sequence of ero- 
sion has proceeded. This [current cycle] 
has only been broken by minor and local 
disturbances.”” He has assigned to early 
Miocene faulting the initiation of a very 
impressive scarp descending toward Ho- 
bart from Mount Wellington (4,000 feet), 
which lies immediately behind the city 
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and is described as a compound faulted 
block thrust up in compression. Conse- 
quent dissection of the major scarp has 
produced “deep narrow gullies often en- 
circled with low cliffs and ending in steep 
walls,” and “the streams” are said to 
“have cut down rapidly to base-level, 
but have not yet had time to widen the 
valleys, thus indicating a fairly recent 
uplift” (Lewis, 1946). Lewis has, how- 
ever, shown the high, bold front of 
Mount Wellington to be, like the great 
scarps of the Launceston district, in part 
a structural escarpment of the thick in- 
trusive dolerite sheet overlying Permian 
strata which dominates the geology and 
landscape of Tasmania. Thus he invali- 
dates the inference of “recent uplift,” 
which is, in any case, inconsistent with 
the rest of the story of the geomorphic 
history. 

The interpretation of at least the ma- 
jor scarps as tectonic rather than of fault- 
line erosional origin, being based on gen- 
eral geological relations, can scarcely be 
questioned; but when it comes to dating 
the earth movements responsible for 
making these features, the earlier guess 
of an end-Tertiary orogeny seems more 
acceptable to a geomorphologist than a 
middle Tertiary guess, unless much of 
the geomorphic story has been misunder- 
stood or still remains untold. Adoption 
of the early Miocene date has made it 
possible to account for the presence of 
plant remains believed to be Miocene in 
lacustrine strata believed to have ac- 
cumulated in hollows consequent on the 
great deformation. 


SOUTH AUSTRALIA 


In an account by Fenner (1930) of the 
landscape features of South Australia, 
he tells that “the Gulf and Torrens areas 
...represent major tectonic features 
that form the southern portion of ...a 





great broken zone extending northwards 
across Australia.”’ The boundaries of the 
main tectonic blocks (horsts and grabens) 
in the southern part of the region he 
shows on a map. He regards these as 
fault scarps and zones of step faulting 
and has described the stepped blocks, 
“particularly in the Flinders Ranges and 
in the eastern Mount Lofty Ranges,” as 
“almost diagrammatic in appearance.” 
The scarps bounding the Mount Lofty 
horst, especially those on the western 
side in the vicinity of Adelaide, have 
been described (Benson, 1911; Sprigg, 
1945) and are well established by evalua- 
tion of both geological and topographic 
evidence. 

In South Australia, unlike eastern 
Australia, large-scale tectonic (conse- 
quent) control of the existing drainage is 
recognized. Rivers of the past, so far as re- 
stored by paleogeographic study, followed 
different lines from the rivers of today; 
but, as a sequel to a late Tertiary fault 
dislocation of the landscape (as well as 
of a coastal strip of compound terrain 
with a cover of Tertiary strata), systems 
of streams have come into existence with 
centripetal arrangement in the Lake 
Kyre, Lake Torrens, and Lake Frome 
basins. These “are purely consequent 
upon and closely reflect the facts of the 
... ranges” (Fenner, 1930). 

In an account published recently of 
the western slopes of the Mount Lofty 
Ranges, Sprigg (1945) has applied a 
statistical method of study based on 
topographic maps, which has enabled 
him to reconstruct generalized contours 
of a deformed surface, and has thus fixed 
the positions of scarps marginal to down- 
stepping faulted benches. Sprigg has re- 
vived Benson’s (1911) early fault-scarp 
explanation of these breaks in the con- 
tinuity of an upland surface interpreted 
as a resurrected fossil peneplain. He 
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does not favor a rather complex alterna- 
tive hypothesis of faulting in two widely 
separated stages with an intervening 
peneplanation, as proposed by Fenner 
(1930) following an oral suggestion made 
by Douglas Johnson. According to this 
hypothesis the upland surface of the 
range consists in part of facets of the 
peneplain developed in the interval be- 
tween two faulting upheavals. According 
to Sprigg’s conclusion, on the other hand, 
the geomorphic cycle that was initiated 
by late Tertiary block faulting still con- 
tinues. It “has removed the bulk of the 
[Miocene] overmass sediments as well as 
much of the early Tertiary erosion sur- 
face’’—a senile surface truncating an- 
cient rocks—which it has resurrected. The 
drainage of this complex block originated 
in consequent fashion on a dislocated and 
deformed surface of the overmass. 


WESTERN AUSTRALIA 

Only brief reference need be made to 
Teichert’s recent suggestion (1947a) that 
all the deformation, including strong dis- 
location, that has been suffered not only 
by Tertiary but also by Mesozoic and 
late Paleozoic strata in Western Austral- 
ia, has taken place in Pliocene times. 
Probably this is based on a misconcep- 
tion of the geomorphic history of the re- 
gion. The suggestion of a major Pliocene 
orogeny preceding peneplanation is en- 
tirely at variance with the sequence of 
events noted by other observers and 
summarized by Jutson (1934). 

The faulting which, Teichert notes 
(19476), breaks Permian and Cretaceous 
strata “into a number of westward-tilted 
blocks” must be of a date anterior to the 
peneplanation of regional extent, which 
Jutson regards as reaching perfection in 
the Miocene period. Upheaval occurred 
after this, probably in Pliocene times, 
developing broad undulation (also of re- 
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gional extent) of aplateau surface; but this 
is broken in a few places by block faulting. 

Woolnough (1920) regards the Stirling 
Range, in the south of Western Australia, 
as a major tectonic landform bounded 
by fault scarps 2,000 feet high; and 
block-faulted plateaus and sunklands 
separated and bounded by long and 
prominent scarps regarded as tectonic 
have been reported in the Northwest 
Division and in the Napier-Geikie Range 
in the far north. 

The most remarkable scarp, however, 
which has been generally attributed with 
confidence to recent dislocation on a large 
scale is that of the Darling fault. This is a 
marginal feature of the land, not orogenic 
in the sense of bounding a horst or other 
form of block mountain uplift; but, if the 
topographic form has been produced by 
dislocation of the land surface, it is, 
nonetheless, a tectonic feature. The fault 
is roughly parallel to the west coast, 
toward which the scarp descends, at a 
minimum distance from it (near Perth) 
of 20 miles. For more than 200 miles, ac- 
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cording to the account by Jutson (1934), 
this line is marked by a prominent scarp 
800-1,000 feet high, of resistant ancient 
rocks, while at the base of this lies a so- 
called “coastal plain,” consisting of thin 
blown sands and alluvia overlying Meso- 
zoic strata. 

Successive rejuvenations of rivers de- 
scending across the scarp have been de- 
scribed (Clarke and Williams, 1926). 
These rejuvenations are attributed to in- 
termittent upheaval of the hinterland as 
the Darling fault scarp has been progres- 
sively developed, and geomorphic evi- 
dence of this kind has led to a presump- 
tion that a tectonic theory of origin is 
correct. Geologists are, however, still 
searching for evidence of dislocation 
sufficiently recent to satisfy them that 
the landscape feature has not been mere- 
ly exposed by fault-line differential ero- 
sion and accentuated by marine erosion. 
Laterite on the foreground surface has, 
it seems, not been detached by faulting 
from the “duricrust” of the inland pla- 
teau above it (Prider, 1948). 
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LEACHING OF CARBONATES IN GLACIAL DRIFT AND 


LOESS AS A BASIS FOR AGE CORRELATION’ 


RICHARD FOSTER FLINT 
Yale University 


ABSTRACT 


The rate of leaching of carbonates in Pleistocene deposits is affected by several factors, each of which 
is variable. Therefore, it is considered that correlation of two deposits of less-than-stage value stratigraphi- 
cally, based on comparable mean depth of leaching in the two, is not likely to be reliable. 


THE PROBLEM 


It has long been recognized by geolo- 
gists as well as by soils scientists that one 
of the conspicuous chemical processes in 
the formation of soil from raw material- 
such as glacial drift and loess—rich in 
calcium carbonate, is the gradual solu- 
tion of calcium carbonate by downward- 
percolating water and its reprecipitation 
(especially under subhumid and semiarid 
climates) at greater depth below the sur- 
face. Both the leaching and the reprecipi- 
tation ordinarily take place above the 
water table. The dual process is particu- 
larly evident in the drift and loess in the 
Central Lowland region, for there the 
Pleistocene deposits include a_ large 
quantity of limestone and dolomite 
ground up by the glacial mill, and the 
climate is sufficiently moist to promote 
the active circulation of subsurface wa- 
ter. East of the Lowland the drift is gen- 
erally less rich in carbonates; west of it 
the climate is too dry to permit rapid 
leaching. 

Movement of calcium carbonate is one 
of many processes involved in the forma- 
tion of soil from parent material. The 
relative intensities of these processes are 
controlled, in turn, by a number of basic 
factors. These controls have been 


‘Published by permission of the Director, U.S. 
Geological Survey. Manuscript received September 
23, 1948. 


97 


grouped in various ways. According to 

Thorp (1941) they are grouped thus: 

1. Parent-material 

2. Climate 

3. Biologic activity (the successive organisms 
that operated on and in the soil during its 
formation) 

4. Relief (topographic position) 

5. Time 

To this list might be added minor factors 

in the geologic environment such as posi- 

tion of water table. 

Modern opinion assigns a greater im- 
portance to climate, relative to time, 
than was formerly thought probable, and 
some soils are now believed to record the 
influence of climatic fluctuation during 
their formation (Bryan and Albritton, 
1943). These fluctuations have common- 
ly involved changes in vegetation covers 
and in the positions of water tables. 
Probably most of them are phenomena 
of the glacial and interglacial ages. 

The belief has been held by many 
geologists that the thickness of the 
leached zone increases, with the passage 
of time, to a maximum imposed by aver- 
age depth of water penetration. In this 
relationship, it has been reasoned, there 
may be a means of determining the rela- 
tive ages of two or more drift sheets. 
Specifically, within one climatic zone and 
in drifts of similar composition, the dif- 
ferences in thickness of the leached zone 














in two different sheets of drift should in- 
dicate the relative ages of the two drifts. 
This view stresses the importance of the 
time factor and gives little consideration 
to the other factors that influence soil 
making. 

The evaluation of depth-of-leaching 
measurements demands an account of the 
leaching process and of the attempts 
made to use depth of leaching as a guide 
to stratigraphic correlation. 


THE PROCESS 


As rain water soaks into the ground, 
it is far from saturated with calcium bi- 
carbonate (the substance present in the 
dissolved state) and it contains carbon 
dioxide derived from the air and from the 
humus zone at the surface. Calcite is 
readily soluble in water charged with 
carbon dioxide, and solution proceeds 
from the upper surface downward, until 
a state of saturation is reached. At the 
level at which saturation occurs, cal- 
cium carbonate is precipitated as a pow- 
dery filling of soil interstices, as concre- 
tions, and in some places as thick, rather 
continuous beds. Secondary carbonate is 
thereby added to the primary carbonate 
already present at this level. 

The factors that control precipitation 
of secondary carbonates are not well un- 
derstood. However, the chief process in- 
volved appears to be evaporation. In sub- 
humid regions much downward-perco- 
lating water evaporates before it reaches 
the water table. Evaporation seems to 
occur in two distinct ways: (1) by ordi- 
nary evaporation within permeable soil, 
aided under suitable conditions by capil- 
larity and in the drier regions by desic- 
cation cracks; and (2) by transpiration of 
soil moisture by plants, the roots of 
which, even in subhumid regions, may 
penetrate several feet below the surface. 
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VARIABLE FACTORS 


If we apply the factors known to in. tot: 


, : ) occ 
specific process of leaching of carbonate | ° 


fluence soil formation in general to the 


we can list a number of specific variabk 
factors. 

Parent-material: texture and stratifica. 
tion.—The conspicuous textural differ. 
ences most commonly met with are the 
differences between clay-rich till and 
sandy till, and between till and loess. In 
the Central Lowland the tills commonly 
are rich in clay and very compact; hence 
their permeability to percolating water is 
relatively small. In contrast, the loess 
sheets consist predominantly of silt and 
are less compact than the tills; hence 
their permeability is greater than that of 
till. It has been estimated that in north- 
western Illinois, loess undergoes leaching 
at about twice the rate applicable to till 
(Leighton, 1923, p. 280). In a stratified 
section, in which a body of till is overlain 
by a layer of surficial loess thinner than 
the leached zone, the thickness of the 
loess strongly influences the depth of 
leaching. In post-Iowan drift in south- 
eastern South Dakota this situation oc- 
curs commonly, with 2-4 feet of loess 
overlying clay-rich till. Where surface 
drainage is poor, the loess is almost in- 
variably leached from top to base; and 
the depth of leaching, if any, in the u- 
derlying till increases with decreasing 
thickness of the loess. In materials of 
comparable permeability, solution pro- 
ceeds more rapidly in small particles of a 
given carbonate substance than in large 
particles. 

Parent-material: —composition—espe- 
cially its content of primary carbonate.— 
It is evident that, other factors being 
equal, increase in thickness of the leached 
zone will take place at a slower rate in 
material rich in primary carbonates than 
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in material poor in these compounds. 
The richer the material, the greater the 


occur per unit of thickness; in rich mate- 
rial the percolating water becomes satu- 
rated with calcium bicarbonate quickly; 
this reduces the rate of increase in the 
thickness of the leached zone (though it 
does not necessarily reduce the rate of 
the leaching process itself). The solution 
of dolomite, of course, takes place more 
slowly than that of calcite because of its 
lesser solubility. 

Climate.—Significant climatic changes 
are known to have affected wide regions 
during the glacial and interglacial ages, 
including the time since the climax of the 
latest glacial age. Such changes may be 
expected to have influenced the rate of 
leaching directly through changes in 
precipitation and temperature and indi- 
rectly through changes in vegetation and 
in the position of the water table. Of 
these components, probably the most im- 
portant is precipitation. 

Other factors being equal, large pre- 
cipitation, distributed more or less uni- 
formly throughout the year, should result 
in a greater rate of increase in thickness 
of the leached zone than should small 
precipitation sharply concentrated in 
time. It is difficult to find comparative 
data because areas in the Central Low- 
land with markedly different precipita- 
tion are normally so far apart that the in- 
fluence of other factors on depth of 
leaching cannot be minimized. Another 
complicating factor is the capillary move- 
ment of water upward, resulting in evap- 
oration and in deposition of carbonates 
in the soils of areas having a mean an- 
nual precipitation of less than about 28 
inches of water and with a well-marked 
dry season each year (Howard, 1946, 
1947). 
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Character of vegetation—It has been 
learned empirically that in some regions 
at least, and with other factors equal, 
thickness of the leached zone is greater 
beneath a forest cover than beneath a 
grass cover. This confirms the theoretical 
consideration that downward-percolat- 
ing subsurface water beneath a forest is 
more acid than that beneath grassland. 
For example, in three counties in central 
Illinois, loess samples taken from _be- 
neath forest vegetation show a mean 
depth of leaching more than 50 per cent 
greater than that beneath grass vegeta- 
tion (Smith, 1942, pp. 170-171). 
Topography—especially surface slope. 
It has been noted repeatedly that with 
other factors equal, the thickness of the 
leached zone becomes thinner with in- 
crease in the surface slope directly above. 
In fact, beneath thoroughly dissected 
topography possessing no nondrained or 
poorly drained areas, till and loess com- 
monly exhibit no leached zone at all. 
This is true even in old drifts such as 
the Kansan. Presumably this condition 
is the result of the large proportion of 
local precipitation that runs off the sur- 
face. The large proportion of runoff not 
only reduces percolation into the soil and 
consequent solvent action but also per- 
mits mechanical erosion of the soil to 
take place, thereby reducing the thick- 
ness of the leached zone or inhibiting it 
altogether. In some places deflation also 
may have played a part in removing an 
unknown upper portion of the soil. 
Position of water table.—Since leaching 
generally occurs only above the water 
table,? the water table ordinarily sets a 
lower limit to leaching. Where the water 
table is very close to the surface, leaching 





21It is possible for leaching to occur below the 
water table if the ground water is not saturated with 
calcium bicarbonate. 
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is not apparent. On the other hand, in 
subhumid regions the water table gen- 
erally lies far below the zone of precipita- 
tion of secondary carbonates, and there 
is no apparent connection between the 
two features. From these statements it 
appears that in some situations the posi- 
tion of the water table may influence the 
depth to which leaching penetrates. 

Time.—It is not known with certainty 
whether increasing depth of leaching is a 
straight-line function of time or whether 
its relation to time is represented by a 
curve. However, we may deduce that the 
rate of increase of depth of leaching prob- 
ably diminishes as time progresses, for 
two reasons. The first is the necessity 
of redissolving ever increasing amounts 
of secondary carbonate precipitated ear- 
lier. The second is the loss of downward- 
moving water abstracted through the 
roots of plants and transpired into the 
atmosphere, as the water passes through 
an ever increasing thickness of soil ai- 
ready leached. Whether or not this de- 
duction is valid, the quantitative rela- 
tion between depth of leaching and time 
is not known; hence the results of meas- 
urement of the leached zone yield values 
that at best are only relative. 


SUMMARY OF VARIABLE FACTORS 

It was stated near the beginning of 
this paper that attempts had been made 
to derive relative time figures from meas- 
urements of depth of leaching. However, 
consideration of the variable factors in- 
volved makes it clear that the results of 
such calculations must include a wide 
margin of error. 


MINIMIZATION OF VARIABLES 
Realizing the presence of variable fac- 
tors, investigators have attempted to 
reduce them to a practical minimum. In 
consequence they have achieved results 
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which, though not accurate, indicate dis. 
tinct age differences between pairs of 
drift sheets known independently to haye 
conspicuously different ages. Reduction 
of variable factors consists mainly of (1) 
selecting drift sheets of comparable com- 
position and general texture, (2) selecting 
drift sheets in the same district with re- 
spect to rainfall, (3) selecting test points 
only on broad, flat uplands with negligi- 
ble surface slope, (4) avoiding the com- 
parison of data obtained from beneath a 
forest cover with those obtained from be- 
neath a grass cover, and (5) using auger 
borings rather than surface exposures for 
sampling, so as to minimize the possi- 
bility of contamination. Of these points 
1 and 3 are believed to be the most im- 
portant and 5 seems the least important. 
The writer’s samplings in the Wisconsin 
drift in South Dakota suggest that there, 
at least, it is almost negligible. 
NECESSITY FOR MANY SAMPLINGS 
Despite the exercise of care in reduc- 
ing the variability of contributing fac- 
tors, the resulting figures, within any dis- 
trict, show a wide range of variation 
among individual test points. Hence sta- 
tistical means derived from many indi- 
vidual tests have had to be resorted to in 
order to show consistent differences be- 
tween two drift sheets. In Indiana 
(Thornbury, 1940) 51 samplings of the 
thickness of the leached zone in Illinoian 
drift close to the border of Kay and 
Leighton’s Tazewell drift showed a mean 
of 158 inches, but with an extreme range 
of 126 inches (102-228 inches) or about 
80 per cent of the mean. Possibly this 
great range and that in the Tazewell 
drift cited below are in part the result of 
leaching under a variable thickness of 
loess. Corresponding figures for 64 sam- 
plings of the Tazewell drift, close to the 
same border, are: mean = 58 inches, ex- 
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treme range = 54 inches (36-90 inches), 
or about 93 per cent of the mean. Al- 
though there is a gap of only 12 inches 
between the minimum figure for the IIli- 
noian and the maximum figure for the 
Tazewell, the frequency distribution of 
the samplings’ supports the validity of 
the mean figures given. 

On the basis of the foregoing figures, 
thickness of the leached zone would seem 
to be a useful factor in distinguishing be- 
tween stratigraphic units, such as the 
Illinoian and the Wisconsin drifts, hav- 
ing the value of stages. In particular, 
depth of leaching should be of value in 
fixing the position of the border of the 
younger drift in detail where, as is often 
the case, no end moraine is present at the 
border. At any rate, probably it would 
be of value only within a geographically 
restricted region. 

The mean value of the Illinoian drift 
samplings in Indiana (158 inches) con- 
trasts with a mean value of 150 samplings 
of the same drift sheet in western Illinois, 
of only 97 inches (Leighton, 1923). The 
most obvious variable involved here is 
mean annual precipitation, which is 
about 40 inches in the Indiana district 
sampled and only about 34 inches in the 
Illinois district sampled. Whether or not 
the marked difference in depth of leach- 
ing is to be correlated with the difference 
in precipitation or with other factors such 
as thickness of overlying loess, the vari- 
ability of results from district to district 
is noteworthy as a barrier to the useful 
employment of leaching-depth data in 
the determination of age differences. 

Within a single drift sheet, the Taze- 
well, large groups of samplings were 
made between pairs of successive end 
moraines presumably distinguished by 


3 Individual samplings are not given in the origi- 
nal reference; W. D. Thornbury kindly supplied two 
sets of unpublished individual values. 
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minor age differences. In Indiana the re- 
sults (Thornbury, 1940, p. 454) were as 
shown in table 1. The frequency distribu- 
tion of the individual samplings cannot 
be assessed because only the means are 
given. The extreme range of these means 
is 20 inches (38-58 inches), a sizable pro- 
portion of the mean thickness derived 
from the four figures given in table 1. 
The data suggest, as was pointed out by 
the investigator who obtained them, the 
reasonable inference of gradual decrease 
in the age of the Tazewell drift from 
south to north. But they suggest also 
that, to be useful in correlation, such 


TABLE | 


Leached 


Areal and Stratigraphic Position Zone* 


Immediately behind the Shelbyville 


moraine (mean of 64 samplings)... 58 
Shelbyville drift (mean of 181 sam- 
a NAAR Pos a 51 
Champaign drift (mean of 111 sam- 
it Ri ee ety 42 
Bloomington drift (mean of 100 sam- 
Er Ad AS sivesich Satie WAR Ca 38 


* Mean thickness in inches. 


data could be drawn only from a drift 
sheet so thoroughly studied that the in- 
dividual minor layers of which it is com- 
posed, corresponding to individual end 
moraines, are known in detail. 


WISCONSIN DRIFT IN ILLINOIS AND IOWA 


The results obtained in northwestern 
Illinois (Leighton, 1923) and eastern 
Iowa (Leighton, 1923; Kay and Graham, 
1943, pp. 107-108, 195-196; Alden and 
Leighton, 1917, pp. 81-824) illustrate the 
difficulty encountered in attempting to 
correlate units having the small value of 
substages as contrasted with stages. The 

4 The data from Alden and Leighton are not com- 
parable with the data from the other sources because 
they are derived from leaching depths measured not 


from the surface but from the top of the till beneath 
a capping of loess. 
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data are summarized in table 2. The 
Iowan and Tazewell substages are fairly 
well separated from each other on a basis 
of mean thickness. On the other hand, 
the drift of the Green River lobe, on the 
basis of these figures alone, would cer- 
tainly be classed as Iowan. Yet Leighton 
(1923, p. 280) found the loess overlying 
the till of the Green River lobe to be 
notably thicker than the loess overlying 
either the Iowan or the Tazewell drifts. 
This difference, largely the result of the 
geographic position of the Green River 
lobe, influenced the leaching figures sig- 
nificantly because loess, being more per- 
meable than till, permits more rapid 


TABLE 2 
Substage and Areal Position — 
Iowan drift in eastern Iowa: 
146 samplings (Leighton)......... 62 
37 samplings (Kay and Graham).. 66 
Tazewell drift in northwestern Illinois: 


Belvidere lobe, 120 samplings..... 48 

Bloomington drift, 56 samplings... 44 
Drift of undetermined age: 

Green River lobe, 119 samplings... 68 

* Mean thickness in inches. 


leaching than does till. Taking this factor 
into account, Leighton judged the corre- 
lation of the drift of the Green River lobe 
to be indecisive; later he (Flint and 
others, 1945) dated it as Tazewell. 

In an excellent report on a detailed ex- 
amination of Tazewell drift in northern 
Illinois, Willman and Payne (1942, pp. 
176-177) presented a table of data on 
the thickness of the leached zone but 
avoided correlation of these data with 
specific drift sheets, apparently because 
of the variable factors involved. The 
table clearly shows that depth of leach- 
ing is greater where surficial silt (presum- 
ably loess) overlies sand or gravel than 
where it overlies till, because silt is much 
more permeable than till. 

Samplings of two drifts of post-Iowan 








age in southeastern South Dakota ma 
by the writer in 1946 and 1947 yielded 


the results given below. The two drifts | 
are of no more than substage value. The | 
older drift is covered with 2-5 feet of F 


loess, whereas the younger is not; other. 
wise the factors seem comparable. 
Twenty-five samplings of the older oj 
the two drifts (Cary?) near its border in 
Lincoln, Minnehaha, Turner, and Bon 
Homme counties gave a mean depth of 
leaching of 29 inches with an extreme 
range of 22 inches (18-40 inches). Fifty- 
eight samplings of the younger of the 


two drifts (Mankato?) near its border in 


Hutchinson, Yankton, and Bon Homme 


counties gave a mean value of 23 inches, f 


with an extreme range of 48 inches (0-48 
inches). The frequency-distribution plots 
made from the two sets of data fail to 
show a close grouping near the mean; ir 
consequence it seems doubtful that the 
arithmetic means are reliable indicators 
of age difference between the two drifts, 

Eighty-four samplings of the older 
drift (Cary?) in Charles Mix County, 
just west of the counties named above, 
gave a mean of 15 inches with an extreme 
range of 36 inches (o—36 inches). The fre- 
quency-distribution curve is very irregu- 
lar and unsatisfactory; hence the meat 
figure is believed to have little value. The 
irregularity is thought to be the result of 
variable composition of the till, which 
contains Niobrara chalk sharply varying 
in amount from place to place. The 
drifts sampled in the counties to the east 
contain relatively little chalk and are 
more nearly uniform in composition. 10 
this case, therefore, variable lithology is 
the chief disadvantageous factor. 


CONCLUSION 
The data given from South Dakota, 


added to the more extensive data cited 
from Indiana, Illinois, and Iowa (which 
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include nearly all the detailed published stage—although inaccurate, are reflected 
references), are believed to support the in depth-of-leaching figures. For time dif- 
ifts © conclusion deduced in the body of this ferences of lesser value the inaccuracies 
The J paper, namely, that variable factors pre- created by variable factors make for re- 
off) cude the accurate determination, sults of doubtful reliability. 
yer through depth of leaching, of the relative ACKNOWLEDGMENTS.—Helpful suggestions 
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1 of Fences corresponding in stratigraphic val- 
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GEOLOGICAL NOTES 


LIQUID INCLUSIONS IN GRANITE THERMOMETRY‘ 


STURGES W. BAILEY 


University of Wisconsin 


INTRODUCTION 

Microscopic liquid inclusions have been 
used by a number of investigators to indi- 
cate temperatures of crystallization. This 
method of thermometry consists of expand- 
ing the liquid by heat until the contraction 
vacuole normally present in each inclusion 
disappears. The temperature at which the 
vacuole disappears is the original tempera- 
ture of inclusion, neglecting pressure con- 
siderations. 

The method assumes (1) that the liquid 
completely filled its cavity at the time of 
crystallization; (2) that the liquid is an 
aqueous solution uniform in concentration 
from inclusion to inclusion within a single 
crystal and containing no carbon dioxide; 
(3) that the pressure on the liquid at the 
time of inclusion was small or that its mag- 
nitude can be estimated; (4) that primary 
liquid inclusions can be distinguished from 
secondary inclusions; (5) that there has been 
no significant change in the volume of the 
cavity itself due to pressure, solution, or 
precipitation; (6) that there has been no 
addition or loss of liquid from the cavity; 
and (7) that the samples used are represen- 
tative of the deposit under study. 

These assumptions have been discussed 
in recent works of Newhouse (1933), Inger- 
son (1947), and Twenhofel (1947). The as- 
sumptions appear to be justified in certain 
cases where moderate temperatures and 
pressures are involved. This paper reports 
an attempt to extend the method to the 
quartz of several selected granites, in the 
formation of which higher temperatures 
and pressures may be involved. 


* Manuscript received October 6, 1948. 
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APPARATUS 


The essential apparatus for investigating 
liquid inclusions consists of a heating stage 
for the microscope, a heating medium, and 
a temperature indicator. Air was found to 
be the most satisfactory heating medium 
for studying high-temperature inclusions, 
Its only disadvantage is that the rock sec- 
tion must be highly polished to see into the 
interior. In the apparatus employed, a slow- 
moving stream of dehydrated air under con- 
stant pressure passes through an alundum 
tube containing a nichrome wire-resistance 
coil. The heated air then passes through a 
thin, insulated cell of sheet tin, which rests 
on the microscope stage. The cell has an 
inside diameter of 0.6 inch and an inside 
depth of 0.2 inch. The section to be studied , 
rests on two horizontal rungs in the center 
of the cell, so that the air flows equally over 
the top and the bottom. Transmitted light 
enters the cell through a glass plate in the 
base. A thick cover glass fits over an open- 
ing in the top. When necessary, air jets set 
in the side of the insulating container may 
be played on the microscope objective to 
keep it cool. 

During operation, each rock section is 
heated gradually to the temperature of 
vacuole disappearance. The controlled-cur- 
rent flow in the heating coil regulates the 
temperature of the air flow. Once equi- 
librium has been established, the tempera- 
ture of the air flow remains constant for a 
given amperage. Ten minutes are allowed 
for the rock section to reach the temperature 
of the air stream. 

Temperatures are measured by a thermo- 
couple, which may be calibrated with re- 
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spect to steam and to the melting-points of 
standard tin, bismuth, lead, and zinc metals. 
Temperature readings are taken by insert- 
ing the thermocouple through a small hole in 
the side of the cell and placing the tip in 
contact with that portion of the section be- 
ing studied. Readings of both the base and 
the top of the section are made, and the 
average of the two readings is taken as the 
temperature of the section. With a section 
of optimum thickness the difference is not 
over a degree or two at 400° C. 

It has been found advantageous to use a 
fairly thick rock section rather than a stand- 
ard 0.03-mm. thin section. Besides providing 
more inclusions in one section, the thicker 
plate allows better determination of the rela- 
tionships of individual inclusions to planes 
and of planes to cracks or crystallographic 
surfaces. On the other hand, a thick section 
is less transparent and requires longer to 
reach equilibrium with the air current. The 
optimum thickness for work with quartz 
was found to be about 0.2-0.4 mm. 

A magnification of three to four hundred 
times is required for work with the small in- 
clusions in granites. A fairly large working 
distance is recommended, however, in order 
to avoid overheating the microscope objec- 
tive. The writer used a 21-power objective 
and a 20-power hyperplane ocular. 

DISCUSSION OF RESULTS 

Sections suitable for liquid-inclusion 
analysis were prepared from fifty-two differ- 
ent granites from the University of Wiscon- 
sin petrographic collection. Nine of these 
granites were selected for detailed work be- 
cause of the abundance of inclusions, the 
clarity of the quartz grains, or a favorable 
geologic history. For example, the Wausau 
(Wis.) granite was selected because of the 
extreme abundance of liquid inclusions. The 
Barre (Vt.) granite has few inclusions but 
is probably free of later metamorphism. The 
Mellen (Wis.) sample came from an inch- 
wide granite veinlet in gabbro. The veinlet 
must have been liquid and cannot be re- 
garded as a product of granitization. 

The liquid inclusions present are smaller 
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in diameter than those in most hydrothermal 
or pegmatite quartz bodies. Under low mag- 
nification they appear as cloudy layers and 
speckled bands. Under higher magnification 
these specks resolve into individual inclu- 
sions which lie in or close to definite planes 
and curving surfaces. 

Differentiation between primary- and 
secondary-type inclusions proved extremely 
difficult. Primary inclusions are those which 
form at the same time as the enclosing 
crystal. They are normally considered to be 
related to growth or crystallographic planes 
and directions; but this relationship may not 
be evident unless the crystallographic planes 
and directions are marked by zonal lines, 
color bands, solid inclusions, or other evi- 
dence. Because quartz is the last constituent 
of granite to form, it fills the available 
interstices and seldom shows crystal form or 
zonal growth lines. Secondary inclusions 
are those introduced subsequent to con- 
solidation of the portion of the crystal in 
which they occur. Laemmlein (1929) has 
shown that they may form by the sealing-off 
of fractures containing a liquid. Dale (1923) 
suggests that the rift and grain of granites 
are due to closely spaced, strain-induced 
sheets of secondary liquid inclusions in the 
quartz. It is also possible that secondary in- 
clusions are due to later re-working or re- 
crystallization. 

No true criteria exist for differentiating 
primary and secondary inclusions in the ab- 
sence of crystal outlines and zonal growth 
lines. Planes of primary inclusions never 
cross one another and do not cross grain 
boundaries. Planes of secondary inclusions, 
however, may cross one another and pri- 
mary planes as well. They may or may not 
cross grain boundaries. Planes of secondary 
inclusions may be superimposed upon planes 
of primary inclusions. Similar orientation of 
planes of inclusions over an entire section 
and gradation into fractures point to a 
secondary origin. Inclusions with branching, 
tubular shapes are secondary; but inclu- 
sions with other shapes, including those of 
negative crystals, are indeterminate. Opti- 
cal orientation is not readily applicable to 








planes in quartz, yielding, at best, only their 
positions relative to traces of the prism 
faces, which are parallel to the c-axis. Planes 
of inclusions, therefore, can be related to 
crystallographic directions but not to defi- 
nite crystallographic planes. 

By use of these criteria, at least half the 
liquid inclusions examined could be classed 
as definitely secondary. The remainder 
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low, but they are primarily important for 
showing that there is considerable variation 
in the temperatures at which the vacuoles 
disappear, even within the same plane of 
inclusions. In all the fifty-two granites ex- 
amined, no planes could be found whose in- 
clusions showed any visual uniformity in the 
size of vacuole relative to cavity. No zoning 
in temperature within or between planes 


TABLE 1* 


RESULTS OF LIQUID-INCLUSION THERMOMETRY: TEMPERATURE OF VACUOLE DISAPPEARANCE 
































Loca.ity 
| | = ae 
Rhyolite Grouse 
Wausau, Mellen, Barre, | Knorre, Ane, Mtn., Mtn., Petershead, Rockport 
Wisconsin Wisconsin | Vermont | Meissen, Saxony Cripple Cripple Scotland Mass. 
. | ‘ | | Saxony — Creek, Creek, . — 
| } | Colorado Colorado 
170° C. (1) 222°C. (1) | 208°C. (1) | 192°C. (1) | 195°C. (1) | 205°C. (2) | 192°C. (1) | 245°C. (2 216° C. (1) 
240 (1) 248 (1) | 216 (1) | 200 (2) | 215 (2) | 220 (1) 231 (1) | 262 (1) | 233 (1) 
300 (1) 260 (2) | 232 (2) | 208 (2) | 221 (1) | 235 (1) | 252 (1) | 283 (1) | 250 (2) 
305 (1) | 2098 (1) 215 (1) 250 (2) 286 (1) 
226 (2) 300 (1) 241 (1) 278 (1) 218 (1) 160 (1) 330 (1) 228 (1) 
243 (1) | 235 (1) 338 (2) | 260 (2) | 280 (1) | 235 (1) | 200 (1) 241 (1) 
267 (2) | 276 (1) | 340 (1) 310 (1) | 292 (1) | 250 (2) | 220 (1) | 105 (1) | 262 (1) 
298 (1) 368 (1) 330 (1) | 278 (2) | 223 (1) | 218 (1) | 265 (1) 
216 (2) | 350 (1) | 384 (1) | 223 (1) | 358 (1) | 280 (1) | 250 (1) | 265 (1) 
225 (1) | 251 (1) 335 (1) | 298 (1) | 278 (1) 
260 (3) | 198 (1) 2 (2) | 227 (1) 
270 (2) | 230 (2) 243 (1) | 180 (1) | | 240 (1) 
340 (1) 258 (1) 200 (1) 260 (2) 200 (1) | | 245 (1) 
289 (2) | 206 (1) | 285 (1) | 230 (2) | | 260 (1) 
200 (4) | 310 (1) | 215 (2) 280 (1) | | 275 (1) 
220 (4) 320 1) | 223 (1) | 200 (1) 
240 (2) | 360 (1) | 225 (2) | 242 (1) } 
260 (1) | 240 (3) | 245 (1) 
260 (1) | 276 (1) | 
216 (2) 270 (2) | 304 (1) } 
225 (2) | 291 (1) 310 (1) 
248 (1) | 300 (1) 
320 (1) | 
172 (1) | 
196 (1) 225 (1) 
330 (1) 230 (1) 
| 256 (1) 
| 260 (1) 
267 (2) 
280 (2) | 
300 (1) 
317 (2) 


* Inclusions are grouped according to the planes in which they occur. Numbers in parentheses refer to the number of inclusions 


could not be classified. None were definitely 
primary. 

Because no primary forms could be dis- 
tinguished, a number of both indeterminate 
and secondary-type inclusions were heated 
until the vacuoles just disappeared. The re- 
sults are listed in table 1, in which the in- 
clusions are grouped according to the planes 
in which they occur. Temperatures are accu- 
rate within the range +3°C. No pressure 
corrections have been made. 


The resulting figures are all relatively 








was noted; instead, vacuoles in adjacent in- 
clusions differed by as much as 163° C. in 
their temperature of disappearance. No 
carbon dioxide content was noted in any 
inclusion. 

Many inclusions containing huge vacu- 
oles, indicating leakage of liquid, were not 
run. Their temperatures would be consider- 
ably higher than any of those listed. Leak- 
age during heating was also common. In- 
clusions showing leakage were detected and 
eliminated by measuring the size of every 


























vacuole at room temperature both before 
and after application of heat to the section. 

No conclusions as to temperatures of 
crystallization can be drawn from table tr. 
The irregularity of observed temperatures 
within individual planes of inclusions indi- 
cates that at least one of the assumptions 
of the method is not justified in the present 
case. Variations in temperature may be due 
to leakage, to inclusion of a gas phase either 
above or below the critical temperature, or 
to effects of the little-understood processes 
which cause secondary inclusions. Perhaps 
primary liquid inclusions are rare or lacking 
in granites. If they are present, better cri- 
teria are needed for their recognition. Sec- 
ondary inclusions may form within the 
range of primary crystallization tempera- 
tures, but this cannot be proved until sec- 
ondary and primary inclusions can be dis- 
tinguished. 


The present investigation emphasizes 
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the difficulties attendant upon the method 
of liquid-inclusion thermometry in granites. 
The method may yet succeed, but excep- 
tional care will be required in selection of 
specimens, in evaluation of inclusions, and 
in technique. It is also possible that the 
prevalent concept of thecrystallization proc- 
ess of quartz must be modified. If crystal- 
lizing igneous quartz passes through a 
desiccating gel stage, irregular inclusions 
might be expected. Liquid inclusions form- 
ing within a viscous mass may shrink or 
swell irregularly, especially under conditions 
of nonuniform stress or of resurgent solu- 
tions. 
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SOURCE AREA OF GREAT PLAINS PLEISTOCENE VOLCANIC ASH' 


ADA SWINEFORD 


State Geological Survey, University of Kansas, Lawrence 


INTRODUCTION 

Known occurrences of the Pearlette vol- 
canic ash zone of Pleistocene age (early 
Yarmouthian) are scattered from west-cen- 
tral Texas to southeastern South Dakota 
and from western Kansas as far east as 
south-central Iowa (Frye et al., 1948). The 
volcano and silicic lava field associated with 
this immense outburst of acid glass must 
have been of giant size. The prevailing wind 
direction indicates that it may have been 
‘ Manuscript received October 4, 1948. 





situated somewhere to the southwest of the 
area of Pearlette ash. The Capulin group of 
volcanoes has been suggested as a possible 
source (Landes, 1928), but the cones are 
comparatively small, and the lavas of the 
Capulin region consist of basalts, andesites, 
trachytes, soda trachytes, and phonolites 
(Collins and Stobbe, 1942). However, these 
volcanoes are closer to the areas of Pearlette 
ash than are any other Pleistocene vol- 
canoes, and hypothetical acid tuff and pum- 
ice from the Capulin group may conceivably 
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have been completely buried under later 
flows. Other possibilities which should not 
be arbitrarily dismissed are the large areas 
of extrusive rocks in northwestern United 
States, although the prevailing winds are 
not from the northwest. A more favorable 
source is the Valle Grande volcanic area of 
north-central New Mexico. The distances 
from this area to the known localities of 
Pearlette ash range from 220 miles (Hartley 





County, Texas) to 740 miles (Minnehaha 
County, South Dakota). Their locations are 
shown in figure 1. 


VALLE GRANDE 


The Valle Grande region lies in Sandoval 
County about 30 miles west of Santa Fe and 
is part of the Jemez Mountains, between the 
Rio Grande on the east and the Nacimiento 
Mountains on the west. It is shown on the 
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Fic. 1.—Map showing Valle Grande crater and localities from which Pearlette volcanic ash has 


been examined. Not more than one ash locality is plotted from each county. 
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Jemez, Jemez Springs, and Santa Clara 
topographic sheets and is illustrated on the 
Raisz (1939) landform map as the Sierra de 
los Valles. The feature is described by Ross 
(1938) as follows: “This is an explosive vent 
or rather group of vents, that have built up a 
crater of elliptical shape, about 13 by 17 
miles in extent. This makes it the largest 
crater ever reported. The materials ejected 
are almost exclusively tuffs which reach a 
thickness of at least 1,000 feet on the east 
flank.” Williams (1941, pp. 251-252) regards 
the volcano as a caldera of collapse origin. 

The rim is surrounded by an area about 
30 miles in diameter, consisting of rhyolite, 
tuff, and pumice (Talmage and Wootton, 
1937, Pp. 41). The tuff beds on the eastern 
flank form the Pajarito Plateau, which is cut 
by radiating canyons that reach a depth of 
1,000 feet in tuffs (Ross, 1931). At least one 
of these canyons, Frijole, in Bandelier Na- 
tional Monument, is the site of ancient cliff 
dwellings, and the tuff of this canyon and its 
topographic expression have been described 
by Henderson (1913). On the west side of 
the caldera is a related broad sloping pla- 
teau of tuff and pumice which extends to the 
east flank of the Sierra Nacimiento (Renick, 
1931, p. 70). Workable pumice deposits are 
also found as far north as Rio Arriba County 
(Clippinger, 1946, pl. 1). A generalized ge- 
ologic map of the area has been compiled by 
Church and Hack (1939, p. 617). 

The tuffs, which overlie the late Tertiary 
Santa Fe formation, include a dozen distinct 
beds (Williams, 1941). They are judged by 
Ross (1931) to be comparatively recent in 
age and perhaps as late as early Pleistocene. 
Earlier workers (Lindgren et al., 1910, p. 
151) assigned a Miocene age to the material, 
and, still earlier, Reagan (1903, p. 104) con- 
sidered it to be pre-Cretaceous. 


PETROGRAPHY OF THE TUFF AND PUMICE 

The pumice and tuff are cream colored to 
light gray. The tuff forms vertical walls of 
canyons and weathers to a hard buff or 
brown surface. Iddings (1890, pp. 10-12) 
first described the petrography of the tufts, 
from samples collected by J. W. Powell in 
1887. He also described two varieties of 
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rhyolite from the area: a normal porphyritic 
variety and a mica-bearing rhyolite; and he 
concluded that the tuff was composed of the 
former rhyolite. He described the tuff as car- 
rying porphyritic crystals of quartz and 
sanidine and small fragmentary crystals of 
augite. In some of the tuffs the quartzes are 
rounded, and in others they are dihexahedral 
pyramids with short prism faces. Many of 
the sanidines exhibit a blue iridescence 
(Iddings, 1890, p. 10). 

Iddings describes pumice related to the 
mica-bearing rhyolite, which includes quartz 
and sanidine phenocrysts and also oligoclase. 
The mica consists of small flakes of biotite, 
and the quartz is in crackled grains and has 
an amethyst color. 

For the present study tuff was collected 
from the cliff on the north side of Frijole 
Canyon near the Tyuonyi ruin (samples 1 
and 2) and from a road cut between the 
museum and the checking station in Ban- 
delier National Monument (sample 3). Pum- 
ice was collected from a road cut on State 
Highway 4 on the northeast side of the Valle 
Grande caldera (sample 4) and from the 
same highway on the southeast side of the 
caldera (sample 5). 

Petrographic examination of the glassy 
part of the tuff and pumice was made ac- 
cording to methods previously described 
(Swineford and Frye, 1946). The glass of 
tuff sample 1 has the same pale-buff color 
(Ridgway 17’’f) and refractive index 
(1.500 + 0,001) as the Pearlette. The shapes 
of the shards are similar to Pearlette shard 
shapes, but sample 1 contains many more 
fibrous shards. It also has phenocrysts, up to 
I or 2 mm. in diameter, of subangular to 
rounded quartz and sanidine. There are a 
few crystals of other minerals, particularly 
magnetite; no biotite was observed. 

Sample 2 is somewhat weathered, so that 
the original color is destroyed. The shards 
are similar in shape to those of sample 1, but 
the refractive index of the glass is 1.502 + 
9.001. There are also amber-colored euhedral 
quartz phenocrysts of the habit described by 
Iddings and somewhat weathered sanidine 
phenocrysts. The sample contains a small 
quantity of tiny magnetite grains. 
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Sample 3 is an almost neutral gray 
(Ridgway 15’”’’f), has an index of about 
1.52, and contains many large quartz and 
sanidine phenocrysts. Much of the sanidine 
shows the blue iridescence or sheen noted by 
Iddings, and many of the quartz crystals 
have sharp interfacial angles. Crystals of 
other minerals are rare. 

Pumice sample 4 is fibrous and frothy in 
appearance, and its color and index are the 


GEOLOGICAL NOTES 









mr 


(Ridgway 17’’g), somewhat lighter than 
typical Pearlette. Its refractive index is 
1.502 + 0.001, and the glass is similar to 
sample 4 in general appearance. It contains 
large anhedral grains of clear colorless 
quartz and sanidine, pseudo-hexagonal bio- 
tite books, dark greenish-brown hornblende, 
and many small particles of magnetite, all of 
which can be recognized under the binocular 
microscope. 


TABLE 1 
CHEMICAL COMPOSITION OF PUMICE AND TUFF FROM VALLE GRANDE REGION, AND 
PLEISTOCENE VOLCANIC ASH FROM THE GREAT PLAINS 


(Analyses by Russell Runnels) 


Location | SiO, | ALO, 


| 
| | 
Valle Grande: | 
Ras | Frijole Canyon | 66.77) 15.68 
3 





Road cut, E. side| 74.81) 1; 





| of caldera 








| 

oe: | NE. side of cal- | 70.45] 14.25] 

| dera 

4b* | NE. side of cal- | 71.32] 14.02! 

| dera | | | 

-, | SE. side of cal- | 70.50} 13.66) 

dera | 

Pearlette: | | 

‘7... | Minnehaha Co., | 71.41) 14.00) 

S.D. | 

a | Harrison Co., Ia.| 70.98) 13.01) 

1 oe | Seward Co., Neb.| 71.29} 13. 20] 

20T | Gove Co., Kans. | 71.68} 13 20) 

427 | Custer Co., Okla.) 72.53) 14.14 
48t 


| Briscoe Co., Tex.| 70 80] 12-40) 


* After removal of crystalline fraction. 


same as those of the Pearlette. It contains 
small phenocrysts of quartz and sanidine. 
In this sample, mineral grains were sepa- 
rated by the use of bromoform diluted with 
alcohol to a specific gravity of about 2.45 in 
order to remove the quartz and feldspar as 
well as the heavier minerals from the glass. 
Heavy minerals include magnetite, horn- 
blende, augite or diopside, hypersthene or 
enstatite, pale-green olivine, apatite, a trace 
of corroded biotite, zircon, and others. The 
crystals are enclosed in glass and constitute 
approximately 8 per cent of the sample, by 
weight. 

Pumice sample 5 is a pale cream color 
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t Analyses from Frye et al. (1948). 


Of the five samples examined, 1 and 4 are 
judged to be possible correlatives of the 
Pearlette. As the phenocrysts and most of 
the heavy minerals should settle out from 
the air within a few miles of the volcanic 
vent, the properties of only the glass need be 
considered. Although the tuff (1) and pum- 
ice (4) contain more fibrous material than 
does typical Pearlette, the proportion of 
fibrous shards may also be a function of se- 
lective transport. The brown color of the 
weathered and case-hardened tuff has no 
counterpart in Pearlette deposits, but it may 
be attributed to the presence and alteration 
of hornblende, magnetite, and other iron- 
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bearing minerals in the crystalline part of 
the tuff. 

For further comparison with the Pearl- 
ette, chemical analyses of some of the Valle 
Grande samples were made. These are sum- 
marized in table 1, in which are also included 
some typical analyses of Pearlette ash. The 
pumice sample 4 was analyzed as a whole 
(4a) and also after removal of the crystalline 
constituents (4b). The latter analysis should 
be comparable with those of the Pearlette. 
Except for the slightly larger percentages of 
calcium and magnesium oxides in the New 
Mexico pumice, there seems to be no sig- 
nificant difference in composition between 
the Valle Grande material (sample 4) and 
the Pearlette ash. Sample 1, which is petro- 
graphically similar to the Pearlette, is some- 
what low in silica. 


CONCLUSIONS 


Too few samples were collected at Valle 
Grande to determine the number of different 
tuffs, their stratigraphic order, and their 
relative volumetric importance. Geologists 
who are more familiar with the region are 
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better qualified to answer these questions. 
Valle Grande is judged to be a probable 
source of the Pearlette volcanic ash for the 
following reasons: (1) Only an extremely 
large caldera, emitting highly acid lavas, 
such as Valle Grande, could supply ash 
which is as extensive as the Pearlette. (2) 
The volcano is situated to the southwest of 
the known Pearlette lentils. (3) The petro- 
graphic character and chemical composition 
of the glassy part of some of the tuff and 
pumice are similar to those of the Pearlette. 
(4) The tuffs are said to be early Pleistocene 
in age. 

It is predicted that detailed work in the 
area may show a predominance of the 
Pearlette type of pumice and may definite- 
ly link it to the Pearlette ash zone of the 
central Great Plains and glaciated region. 
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DISCUSSION 


CRUSTAL MOVEMENT IN THE GREAT LAKES AREA: A DISCUSSION: 


WILLIAM HERBERT HOBBS 
University of Michigan 


Under the above title a paper has appeared 
in Bulletin of the Geological Society of America by 
Sherman Moore, engineer of the United States 
Lake Survey Office at Detroit (1948; see also 
Moore, 1922). He has examined the records of 
lake levels at all the gauging stations on the 
Great Lakes, which records cover, on the aver- 
age, a period of fifty-two years, though some ex- 
tend over a considerably longer time. This is 
much the most complete study of its kind for the 
area, and it confirms all the earlier studies from 
Gilbert to Gutenberg in showing that the Great 
Lakes area has been tilted so that its northern 
(really east-northeastern) side is higher than the 
southern one, and that this tilt is still continuing 
(Gilbert, 1898; Freeman, 1926; Gutenberg, 1933). 

Moore has, however, offered the startling and 
wholly unwarranted interpretation that this 
canting of the earth’s surface is due to a sub- 
sidence of the area—a downtilting toward the 
south, not an uptilting toward the north. This 
is contrary to all that has been learned by Gil- 
bert (1898), Goldthwait (1907), Hobbs (1911), 
Leverett and Taylor (1915), Bergquist (1936), 
Stanley (1936, 1937), and other glacialists dur- 
ing the last half-century, from measurements 
(precise levels) of the now abandoned beaches of 
the Maumee, Whittlesey, Warren, Algonquin, 
and Nipissing glacial lakes—all ancestors of the 
present Great Lakes and formerly in occupation 
of the area. 

So far as the gauge records alone are con- 
cerned, a downtilt to the south or an uptilt to 
the north would produce identical changes in 
the water levels measured at the gauging sta- 
tions. This can be illustrated by holding a half- 
filled basin of water in a horizontal position, 
raising one side through a certain distance, re- 
turning it to a horizontal position, and then 
lowering the opposite side by the same amount. 
The movement of the water surface with refer- 
ence to the sides of the basin will be exactly the 


* Manuscript received August 10, 1948. 
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same, whether the basin is uptilted, as in the 
first instance, or downtilted, as in the second, 

It will help to make clear that the crustal 
movement has been an uptilt if reference is 
made to the map of figure 1, which shows the 
outlines of glacial Lake Nipissing, the latest of 
the glacial lakes, whose outline was little differ- 
ent from those of the present Great Lakes, ex- 
cept that Lake Superior was then a part of the 
Michigan-Huron basin. In time Lake Nipissing 
followed several earlier glacial lakes, the latest 
of them Lake Algonquin, during the existence of 
which the hinge line of no vertical movement 
had been migrating northward. The present 
elevations of the isobases of the Nipissing shore 
line, measured in feet above mean tide, are 
shown on the map by the figures enclosed in 
parentheses, and the actual amounts of uplift at 
each of them during the last two thousand or 
more years are shown by the figures not so en- 
closed. It will be noted that, in the Lake Su- 
perior area, Moore’s hinge line of no vertical 
movement and his isobase of maximum sub- 
sidence very nearly correspond in position with 
the isobase of maximum uplift and that of the 
hinge line of no uplift, respectively. 

The rise of the water level at Duluth, which 
Moore has interpreted as due to a subsidence, is 
explained not by subsidence but by greater up- 
lift of 100 feet at the north shore of the lake. 
This has caused the water to be withdrawn from 
that side and has raised the level at the south- 
western shore. A similar withdrawal of water 
from the northeastern shore of Lake Erie, where 
the outlet is likewise to the northeast, and a sub- 
mergence of the southwestern shore have taken 
place for the same reason (Moseley, 1995; 
Hobbs, 1911, pp. 27-35). No such changes have 
occurred in Lake Huron because its outlet is at 
the south end of the lake, so that the displaced 
water passes out through the outlet. 

The withdrawal of the water from the rock- 
cliffed north shore of Lake Superior and the 
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submergence of its southern shore are strikingly 
displayed to the boat traveler along them. On 
the north shore the bottom of the lake is laid 
bare for some distance out from the now aban- 
doned cliffed shore (Lawson, 1893). On the 
southern shore in all places to the south of the 
isobase, which passes through the outlet at 
Sault Ste Marie, the shore is submerged 
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(Stuntz, 1870; Wright, 1905; Hobbs, t1o11, 
Pp. 27-35). 

If we consider the beaches of the earlier gla- 
cial Lake Algonquin, they show a far greater 
amount of uplift. Over the same area, the rate of 
their average uplift is 2.6-3 feet per mile, as 
against 0.5 foot for the Nipissing Lake, which 
we have been considering (Stanley, 1939). 





Fic. 1. 
lift. The hinge line of Moore and his isobar of maximum depression are shown by the heavy and dotted 
sinuous lines (map mainly after Taylor). 


-Map showing isobases of the now warped Nipissing water plane and the amounts of actual up- 
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WOLF MAYNC 
Caracas, Venezuela 


As A. J. Eardley’s paper ‘‘Ancient Arctica” 
(1948) is incomplete and in part in error with 
respect to eastern Greenland, some corrections 
and additional information are offered. 

My detailed stratigraphical studies in 1936— 
1938 of the post-Devonian sediments in eastern 
Greenland justify a note on some errors and 
omissions in Eardley’s paper. Although recent 
works on the stratigraphy of Alaska, the Arctic 
Archipelago, Spitzbergen, etc., are referred to in 
Eardley’s compilation, the geologic literature on 
Greenland from 1939 onward is, unfortunately, 
not taken into consideration. In his bibliogra- 
phy Eardley lists papers presented at the East 
Greenland Meeting in Schaffhausen, Switzer- 
land, on March 11-12, 1939; but the data con- 
tained in these contributions have not been in- 
corporated in his article. 

In order to give an idea of what has been 
achieved during the past years with regard to 
the geology of eastern Greenland, a list of recent 
papers is given at the end of the present note. 

This list of references is supposed to stimu- 
late interest in the geography of eastern Green- 
land, which, beyond all doubt, may today be 
considered the best-known part of the Arctic 
regions. 

Although the so-called “Carboniferous” de- 
posits of eastern Greenland and of other Arctic 
regions have been shown to be Middle to Upper 
Permian in age (Aldinger, 1935; Maync, 1938, 
1939, 1940, 1942), they were still designated 
“Carboniferous” by Eardley. These beds, the 
distribution of which is, moreover, wrongly 
shown on the paleogeographic map (Eardley, 
1948, fig. 7, p. 427), were proved to be in part 
Zechstein formations, lithologically and bio- 
facially identical with the “Magnesian lime- 
stone” and ‘‘Marl slate” of Great Britain, the 
German Zechstein deposits, and the Russian 


* Manuscript received November 6, 1948. 





Kungurian-Kazanian sediments of Pae-Choi, 
Kanin Peninsula, Timan, Pinega Basin, Kazan, 
etc. (Maync, 1938, 1939, 1940, 1942). Further. 
more, these reefal Zechstein deposits of easter 
Greenland were found to be stratigraphically 
associated with fossiliferous Productus lime. 
stone and ammonoid-bearing Martinia lime- 
stone (with Cyclolobus, Medlicottia) of undoubt- 
ed Upper Permian age (Maync, 1938, 10939, 
1940, 1942; Miller and Furnish, 1940). Another 
facies member of the eastern Greenland Permian 
is the Posidonomya shale. These shales carry a 
ganoid fish fauna which contains diagnostic 
Zechstein forms known from the German 
Kupferschiefer, the English ‘‘Marl slate,” and 
the Kazanian of Russia (Aldinger, 1935, 1937). 
This fish fauna, moreover, shows close affinities 
to that of the Middle Phosphoria of Wyoming, 
a formation which is proved by its cephalopod 
faunas to be equivalent to the Word (and Capi- 
tan?) formations of Texas (Middle to Upper 
Permian), and by its foraminiferal fauna to be 
of Zechstein age. 

The occurrence of these genuine Zechstein 
deposits in eastern Greenland, of course, proves 
a marine connection with England and Ger- 
many; and the wide distribution of Permian 
(formerly Carboniferous or Permo-Carbonifer- 
ous) beds all around the Arctic region suggests 
that the present-day North Polar Sea was the 
center of this late Permian Arctic Basin. A 
broad gulf (Scandic Sea) extended southward 
between Greenland and Spitzbergen across the 
British Isles and terminated in the inland Zech- 
stein Sea of Germany (Koch, 1935; Mayne, 
1940, 1942). The “Scandinavian-Scottish Con- 
tinent” (Stamp, 1947; Eardley, 1948) of the 
Lower Carboniferous had ceased to exist. Its 
late Paleozoic remnant was the Scandinavian- 
Baltic Shield, which apparently has never borne 
any marine Permo-Carboniferous deposits. We 
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cannot share Eardley’s view, therefore, that 
“the southern extent of the Carboniferous basin 
on the east coast of Greenland is a matter of 
speculation” (1948, p. 427): Nowhere in eastern 
Greenland did any Carboniferous sea encroach 
upon the continental border, and even the sup- 
posed marine ‘‘Upper Carboniferous” beds on 
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seas. The discovery of Cyclolobus, the diagnostic 
Upper Permian ammonoid genus of the Hima- 
layan geosyncline in the eastern Greenland 
Permian obviously discloses that marine con- 
nections existed between the Arctic Basin and 
the far-eastern geosyncline (Himalayas, Ussuri 
Bay, Mongolia, Siberia). 
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Paleogeographical map of the Arctic regions in Permian (‘‘Permo-Carboniferous” or “Car- 


boniferous” auctorum) times (after Eardley, 1948; Gignoux, 1943; and Maync, 1940; 1942). 


Holms- and Amdrups Land (80°-81° N. lat.) are 
younger than Lower Permian (Artinskian) and 
correspond more or less to the Kungurian stage 
of Russia. The eastern Greenland beds are de- 
posits of the Upper Permian sea, which is trace- 
able south into the heart of Europe. Another 
arm, branching from the Permian Arctic Basin, 
ended in the Russian Kungurian-Kazanian 


Figure 1 depicts the supposed distribution 
of land and sea during the Permian (‘‘Carbon- 
iferous” or Permo-Carboniferous auctorum) as 
a supplement to Eardley’s paleogeographical 
map (1948, fig. 7, p. 427). 

Information on the Triassic formations of 
eastern Greenland is furnished by the compre- 
hensive papers of A. Rosenkrantz, L. Koch, E. 









316 


Nielsen (see Koch, 1935), L. F. Spath (19352), 
and H. Stauber (1942). 

The Jurassic faunas of southeastern Green- 
land have been dealt with by L. F. Spath in a 
masterly manner (Koch, 1935; Spath, 1935), 
1936), and a detailed paper on the Jurassic 
stratigraphy of eastern Greenland and the 
paleogeographical relations with other Arctic 
regions was published recently (Maync, 1947). 
A publication by L. F. Spath (1947) is devoted 
to Jurassic and Cretaceous ammonite faunas 
from southeastern Greenland. Until now, beds 
of Bathonian-Callovian, Neo-Oxfordian, lower 
Kimmeridgian, upper Kimmeridgian, and Vol- 
gian (Portlandian) age are known to occur in 
northeastern Greenland. Facial and faunal evi- 
dence shows that the Scandic Sea also existed in 
the Jurassic, and the landmass between Green- 
land, Spitzbergen, and Scandinavia postulated 
by Eardley seems to me, therefore, improbable 
or at least very disputable. 

Erosion after the late Cimmerian block- 
faulting has removed part of the eastern Green- 
land Jurassic, and synorogenetic boulder beds 
of the Volgian locally overlie different portions 
of the older Jurassic series with a sharp angular 
unconformity. 

A monograph dealing with the Cretaceous 
stratigraphy of eastern Greenland and the Arc- 
tic regions was finished some years ago, but its 
publication was considerably delayed because 
of postwar difficulties in Denmark (Mayne, 
1948d). It is to be hoped that its contents may 
furnish the Greenland data needed for a paleo- 
geographical map of the Arctic region. The 
above-mentioned paper contains paleogeo- 
graphical maps of northeastern Greenland in 
both Valanginian and Aptian-Albian (and 
Senonian) times and attempts to correlate the 
circumarctic Cretaceous deposits. 

In northeastern Greenland the occurrence of 
fossil-bearing sediments of Infravalanginian- 
Rjasanian (Subcraspeditan), middle Valangini- 
an (Polyptychitan), Aptian-Albian, and lower 
Senonian age has been ascertained (Koch, 1935; 
Frebold, 1935; Maync, 1938, 1939, 1940, 1948c, 
d; Stauber, 1938, 1939, 1940). Farther south, 
beds of Cenomanian-Turonian age have recent- 
ly been recorded (Spath, 1946). 

While the Lower Cretaceous faunas of east- 
ern Greenland still bear a distinct boreal charac- 
ter, proving the existence of a large isolated Arc- 
tic Aucella sea, the Aptian-Albian faunas con- 
tain cosmopolitan forms which are found in 
Cretaceous beds all over the world. This fact 
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shows the establishment of marine connections 
between the Arctic sea, the Mediterranean 
Tethys, and the Pacific geosynclines. 

Eardley’s figure 8 (p. 429, lower map), show- 
ing the distribution of the Arctic coal-bearing 
“Miocene” (recte Eocene) deposits, does not 
contain the eastern Greenland type locality on 
Sabine Island, discovered in 1869-1870 by the 
German Koldewey Expedition, the flora of 
which was dealt with by O. Heer in his famous 
Flora fossilis arctica (1868-1883). The strati- 
graphical work in 1936-1938 disclosed several 
new localities with Tertiary plant- or coal-bear- 
ing beds on Sabine and Pendulum islands, Wol- 
laston Foreland, and Hold-with-Hope (Mayne, 
1938, 1939, 1940). These Tertiary sediments 
were deposited during the interval between the 
first basaltic extrusions and the overlying basalt 
flows (traps). 

Another famous locality of early Tertiary 
plant-bearing beds, viz., Atanikerdluk on 
Nugssuak Peninsula, western Greenland, is 
sought for in vain on Eardley’s map. 

Since the publication on L. Koch’s paper 
(1935), referred to by Eardley in discussing the 
tectonics of eastern Greenland (p. 431), several 
important works have been published which 
have thrown new light on the structure of east- 
ern Greenland and its history (see the papers by 
Bierther, Biitler, Cloos, Schaub, Stauber, 
Vischer, Wegmann, listed at the end of this 
article). 

Koch’s idea of strong Tertiary faulting along 
the eastern coast of Greenland is not strictly 
correct because all the major faults recorded are 
Paleozoic and Mesozoic in age and have only in 
part been later reactivated. Since the Variscan 
orogeny, no folding has taken place in eastern 
Greenland, but crustal tensions led to an intense 
block-faulting of the continental border. 

The main phase of this step-faulting falls 
within upper Kimmeridgian times: the different 
tectonic blocks faulted en échélon (antithetische 
Schollentreppe) underwent a westward tilting, 
so that tectonic grooves were formed in the 
joints of the tilted blocks (see fig. 2). These de- 
pressions and the adjacent fault scarps have 
played an important role during the late Juras- 
sic and Cretaceous sedimentation (marine 
troughs and basins bordered to the west by 
steep, unstable fault cliffs). Furthermore, intra- 
Cretaceous movements of some of these old 
blocks are recorded (post-Valanginian, post- 
Albian, pre-Tertiary). 

The few critical remarks given in the forego- 
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ing are not intended to discredit the whole of 
Dr. Eardley’s paper, the importance of which is 
duly acknowledged. I know from experience how 
difficult it is to embark upon so delicate an en- 
terprise as the compilation of paleogeographical 
maps of little-known regions. I wish only to 
point out some omissions and misinterpretations 
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of the region that I happen to have studied at 
first hand. The rapid advances in the knowledge 
of eastern Greenland’s geological history during 
the last ten years certainly deserve careful con- 
sideration. The main shortcomings in Eardley’s 
article are due chiefly to ignorance of the essen- 
tial and indispensable literature. 
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ANCIENT ARCTICA: REPLY TO REMARKS! 


A. J. EARDLEY 
University of Michigan 


It is unfortunate that I was not more thor- 
ough in searching the literature of East Green- 
land in my article on Ancient Arctica, and 
missed the publications that appeared during 
and after the war. It is evident that Mr. 
Maync’s kind remarks should be taken as an 
authoritative supplement to the article and 
should be carefully integrated into it by those 
interested in Arctic paleogeography. 

This is an opportune time to call attention to 
still another important and recent publication 
not referred to in the article or by Mr. Maync, 
namely, one by L. R. Wager, entitled “Stratig- 
raphy and Tectonics of Knud Rasmussens Land 
and the Kangerdlugssuaq Region.” It deals with 


‘ Manuscript received January 4, 1949. 


the great plateau basalts and fossiliferous inter- 
beds and the deformation they have suffered. 
The first marine transgression in the area oc- 
curred in Senonian time, and then immediately 
afterward volcanic activity broke out to form 
the ‘lower lavas and tuffs,”’ which may in part 
be lowest Eocene. The greatest igneous event 
then followed and at least 20,000 feet of basalt 
was poured out to give the “plateau basalt 
series.”’ It is all probably lower Eocene in age. 
Shallow-water sediments below and above show 
that sinking, compatible with the thickness of 
the basalts, occurred during or soon after their 
extravasation. Immediately afterward, the 
Plateau Basalts were cast into a grand flexure 
which is now along the coast with the ocean- 
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ward side down. The differential vertical move- the Arctic and north Atlantic were vastly dif. 
ment between Greenland and the Greenland ferent than now, that the changes are due large. 
Sea basin immediately adjacent is 20,000-25,000 ly to vertical movements of the crust in the 
feet. Although the basalts are not present north magnitude of 5,000-25,000 feet, and that the 
and south of the middle east Greenland area, crust of the great region is one of continental 
the topography along the coast, according to composition—not true ocean floor, such as parts 
Wager, suggests comparable movements. This of the Pacific are. The new evidence is much 
recent evidence of profound vertical movements more significant than the faults mentioned in 
in the Greenland area is strikingly in accord the article, many of which, as Maync points 
with the general theory of the article on An- out, are very ancient, and whose significance to 
cient Arctica, viz., that the past geographies of _ the theory is doubtful. 
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of Knud Rasmussens Land and the Kangerd- vol. 134, no. 5. 
lugssuaq region. Geological investigations in 
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Geological Map of the Province of Ontario: Map 


No. 1946-3. Compiled by J. SATTERLY. On- 

tario Dept. of Mines, 1946. 1:1,267,200 or I 

inch to 20 miles, 40 X 51 in. 50 cents. 

Geologists will welcome the publication of a 
new geological map of one of the largest and 
most complex of the Canadian provinces. The 
extensive mapping program of the Ontario De- 
partment of Mines and other agencies during 
the last two decades has made compilation of 
this map a necessity. The scale of the map per- 
mits delineation of more detail than is shown on 
themap of North America. The map of Ontario, 
therefore, is a useful contribution to the geologi- 
cal literature of that province. 

Dr. Satterly exercised commendable restraint 
in his compilation. No effort was made to make 
the Archean rocks of the different parts of the 
province fit into any restricted system of no- 
menclature or correlation. The Archean is di- 
vided into only volcanic rocks (principally Kee- 
watin), sedimentary rocks (Timiskaming and 
other similar series), and basic (Haileyburian) 
and acidic intrusives. 

The Proterozoic rocks are separated into 
three catagories, namely, sedimentary (Bruce 
and Cobalt), sedimentary and volcanic (known 
locally as Animikie, Whitewater, or Keweena- 
wan), and basic intrusives (mainly Keweena- 
wan). The post-Cobalt granitic intrusives (Kil- 
larney), unfortunately, are not separated from 
the Archean; and it is not clear even from the 
map legend that there is a post-Huronian 
granite. The term ‘“‘Huronian” is missing, al- 
though the type area lies within the province. 

Rocks of four Paleozoic systems (Cambrian 
through Devonian) are appropriately shown, as 
is also the Cretaceous. 

A synoptic summary of the mineral produc- 
tion of the province to the end of 1945 is includ- 
ed on this sheet. 

es; 3 





Geological Map of Saskatchewan, Map 895A. 
Canada Department of Mines and Re- 
sources, Mines and Geology Branch, Bureau 
of Geology and Topography, 1947. 

This map, on a scale of 1 inch to 20 miles, is 
another in a series of generalized provincial geo- 
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logical maps. The geology is based for the most 
part upon the work of the Geological Survey 
of Canada, with some information from Im- 
perial Oil Limited. The Drafting and Reproduc- 
ing Division is responsible for the cartography. 

Somewhat more than 50 per cent of the 
map has been colored, but an additional 10 per 
cent is indicated to be underlain mainly by 
Upper Cretaceous formations. Rocks of pre- 
Cambrian age constitute the bedrock of about 
40 per cent of the province, Mesozoic and 
Cenozoic strata about 57 per cent. Productive 
metal mines, as well as developed and undevel- 
oped water-power sites, are indicated by ap- 
propriate type. 

The legend includes 32 map units, of which 13 
are pre-Cambrian, 3 Paleozoic, 12 Mesozoic, 
and 4 Cenozoic. The colors appear to be well 
chosen and, for the most part, match very well 
those of the Manitoba sheet (Map 850A). Sev- 
eral of the colors for pre-Cambrian rocks, how- 
ever, registered somewhat lighter than did 
those of the corresponding divisions on Map 
850A. The map units are indicated by number 
as well as by color, and it is here that difficulties 
arise. For example, map unit 19 in Saskatche- 
wan is the Lower Cretaceous Clearwater forma- 
tion, while on the adjoining Manitoba sheet the 
same number, but different color, is used for the 
Upper Cretaceous Vermilion River formation. 
It seems to this reviewer that a more desirable 
system would be one based on letters, similar 
to that of the United States Geological Survey. 

Detailed work in the large region underlain 
by pre-Cambrian rocks has been confined chief- 
ly to such areas of economic importance as 
Lake Athabaska, Lac La Rouge, and Amisk 
Lake. The difficulties of correlation are apparent 
in a glance at the legend: map units 1 and 3 are 
both chiefly volcanic rocks, 2 and 4 are chiefly 
sedimentary rocks, while 5 represents volcanic 
and sedimentary rocks. The terms “Keewatin” 
and “Timiskaming” are not used in the legend. 
Archean or Proterozoic age is suggested for the 
succeeding three map units, respectively: 
6, complex of sedimentary, volcanic, and intru- 
sive rocks; 7, chiefly basic intrusive rocks; and 
8, chiefly acidic intrusive rocks. It is in this 
group that some discrepancy with the Manitoba 
sheet was noted. Along the Churchill River at 
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the Manitoba boundary, Map 895A records 
chiefly sedimentary rocks associated with acidic 
intrusions (6c—orange). The adjacent portion 
of the Manitoba sheet shows undifferentiated 
granitic intrusions and older gneissic and 
schistose sedimentary and volcanic rocks 
(7c—pink). It seems to the reviewer that the 
Saskatchewan map should have been colored 
pink (8a) rather than that the Manitoba map 
should be colored orange. This is the only place 
where the two maps do not match in color. With 
the exception of map unit 13, which designates 
late Proterozoic diabase masses of limited ex- 
tent, rocks of definite Proterozoic age are con- 
fined to the Lake Athabaska—Cree Lake region. 
Alcock (1936) has mapped in detail the area 
north of Lake Athabaska, recognizing the 
Beaverlodge series of sediments (9), basic in- 
trusives (10), and a Proterozoic granite (11). 
The Athabaska series (12), comprising coarse 
clastics with minor basalt, constitutes about 10 
per cent of the entire map sheet and among the 
pre-Cambrian rocks is second only to the acidic 
intrusives in areal extent. 

Paleozoic rocks are represented by relatively 
small areas of color, but the boundaries of the 
possible extensions of the map units are indi- 
cated by an inconspicuous dotted line. The 
three map units represent generalized groupings 
as follows: 14, Upper Ordovician; 15, Silurian; 
and 16, Devonian. No formational names are 
given. 

Mesozoic rocks are distributed among 
twelve map units, nearly all of Upper Creta- 
ceous age. Map unit 17, found as small scattered 
patches, represents the Swan River group of 
Lower Cretaceous and possibly earlier age. Map 
units 18, the McMurray formation, and 109, the 
Clearwater formation, both of Lower Creta- 
ceous age, are shown only along the Clearwater 
River near the Alberta boundary. The Ashville 
formation, unit 20, listed as “Lower (?) and 
Upper Cretaceous,” is found only near the 
Manitoba boundary. The remaining eight map 
units are all labeled ‘Upper Cretaceous” in age. 
They include the following formations in as- 
cending sequence: 21, Favel formation; 22, Ver- 
milion River formation; 23, Lea Park formation; 


Atcock, F. J. (1936) Geology of Lake Athabaska 
region, Saskatchewan: Canada Geol. Survey 
Mem. 106. 

FuRNIVAL, G. M. (1946) Cypress Lake map-area, 

Saskatchewan: Canada Geol. Survey Mem. 242. 
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24, Ribstone Creek, Grizzly Bear, and Birch 
Lake formations; 25, Oldman formation, Pale 
and Variegated beds; 27, Bearpaw formation: 
and 28, Eastend, Whitemud, Battle, and 
Frenchman formations. Map unit 26, the Rid. 
ing Mountain formation, is designated as cor. 
relative with part of the Lea Park, together with 
the succeeding formations up to and including 
the Bearpaw. It is recognized in the eastern 
half of the province, the more complicated non. 
marine and marine succession being found in the 
western half. It is to be regretted that, in listing 
the formations in both units 24 and 28, the 
legend carries the succession in inverse order, 
the youngest beds being placed at the bottom. 
It is of interest to note that the name ‘“‘Alberta 
formation” does not appear in the legend, al- 
though it is widely used in the literature and has 
been recognized by Wickenden (1941) as cor- 
relative with the Favel and part of the Ver- 
milion River formations. 

Strata of Cenozoic age are mapped in four 
units, only two of which are at all extensive. The 
first of these, 29, the Ravenscrag formation, is 
mapped as Paleocene in the Cypress Hills re- 
gion, where it has been redefined by Furnival 
(1946, p. 106). Elsewhere the same color is used 
with the symbol 29a, where some Upper Creta- 
ceous beds are included. The Swift Current beds, 
30, of Eocene age occur only as isolated patches 
southeast of the city of the same name. The 
Cypress Hills formation, 31, of Oligocene age is 
the second map unit of any extent. It is confined 
to the southwest quarter of the sheet. The sym- 
bols 32 and 32A are used for the Wood Moun- 
tain formation of Miocene age, consisting chiefly 
of gravels. These represent the youngest rocks 
in the province, other than the unconsolidated 
deposits of Pleistocene age. 

This map and the companion sheet, Map 
896A, comprise a set which has already been 
found useful for instructional purposes by the 
reviewer. It is expected that it will be widely 
used. The cartographers have been able to in- 
clude considerable detail without sacrificing 
clarity and legibility, and they have produced 
an excellent work. Additional maps of this type 
will be most welcome. 





WICKENDEN, R. T. D. (1941) Cretaceous marine 
formations penetrated in wells near Lloyd- 
minster, Saskatchewan: Royal Canadian Inst. 
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Mineral Map of Saskatchewan, Map 896A. 
Canada Department of Mines and Re- 
sources, Mines and Geology Branch, Bureau 
of Geology and Topography, 1947. 

This map is the companion sheet, on the 
same base, of the geological map of Saskatche- 
wan (Map 895A). The eight mining districts 
into which the province is divided are outlined. 
Metallic mineral properties are indicated by a 
large spot and clearly legible number, both in 
red. The numbers are distributed through an 
alphabetical index of fifty mining properties, 
corrected to July, 1946. Metallic mineral oc- 
currences are designated by a small red X with 
the name of the mineral in small red capitals. 
Nonmetallic mineral occurrences are marked in 
the same way, but in blue print. Economic de- 
posits of clay of the Cretaceous Whitemud for- 
mation are shown by solid blue patches. In the 
pre-Cambrian portion of the map numerous 
glacial striae are shown by an appropriate arrow 
symbol. 

This sheet will be found to be a useful com- 
plement of the geological map. It should be of 
particular interest to students of both economic 
geology and geography. 

H. S. ARMSTRONG 


Geological Map of Manitoba, Map 850A. Canada 
Department of Mines and Resources, Mines 
and Geology Branch, Bureau of Geology and 
Topography, published 1946. 


This map, on a scale of 1 inch to 20 miles, is 
one of several provincial maps which have been 
issued in recent years. The geology is compiled 
from the work of the Geological Survey, the 
Manitoba Department of Mines and Natural 
Resources, and the Ontario Department of 
Mines. The cartography is by the Drafting and 
Reproducing Division. 

Most of the area between latitudes 40° N. 
and 55° N. has been colored; but, with the ex- 
ception of the area west of Granville Lake, the 
remainder of the province northward to its 
boundary (60° N.) is colored only along some of 
the major waterways. The names of Robert 
Bell, J. B. Tyrrell, D. B. Dowling, and William 
McInnes are associated with this mapping, most 
of it done prior to the turn of the century. Pro- 
ductive metal mines and deposits of nonmetal- 
lics are indicated by distinctive symbols and 
type. ; 
The legend includes twenty-two map units, 
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of which eight are pre-Cambrian, six Paleozoic, 
five Mesozoic, one Cenozoic, and two of indefi- 
nite age. Most of the colors chosen to represent 
the different units are easily recognizable. If the 
geological structure were more complex, how- 
ever, difficulty might be found in distinguishing 
the somewhat similar shades of green used for 
rocks of such diverse ages as Archean, Ordovi- 
cian, and Cretaceous. 

Rocks of pre-Cambrian age, comprising the 
larger part of the map sheet, are divided into 
eight units. The first three, of Archean age, are 
as follows: (1) chiefly volcanic rocks, (2) chiefly 
sedimentary rocks, and (3) volcanic and sedi- 
mentary rocks. This threefold division of the 
Archean rocks, a lithologic rather than a strati- 
graphic grouping, reflects the lack of sufficient 
detail in some of the source material. It is note- 
worthy that the terms “‘Keewatin” and “Timis- 
kaming” do not appear in the legend, although 
Lawson’s original ‘““Keewatin” area lies just east 
of the southeast corner. Possible Archean or 
Proterozoic age is suggested for each of the suc- 
ceeding four map units: (4) chiefly sedimentary 
rocks, (5) complex of sedimentary, volcanic, and 
intrusive rocks, (6) chiefly basic intrusive rocks, 
and (7) chiefly acidic intrusive rocks. The alter- 
native ages suggested recognize the disagree- 
ment existing among Canadian geologists con- 
cerning the age of the sediments of unit 4. The 
intrusives of unit 7 are further subdivided and 
include gray sodic types as well as red potassic 
types. This subdivision is reminiscent of the 
similar twofold division of granites in northern 
Ontario (Moore and Charlewood, 1930). The 
only rocks to which a Proterozoic age is defi- 
nitely assigned are the sediments, mostly 
quartzites, of the Churchill formation, colored 
as map unit 8. 

Palaeozoic strata form a broad central band 
extending northwesterly from the International 
Boundary. Rocks of this age occur also along 
several rivers west and southwest of Hudson 
Bay. The Ordovician, Silurian, and Devonian 
systems are represented; but in only one of these 

-the Ordovician—is any further subdivision 
made. The earliest Paleozoic rocks are Upper 
Ordovician, comprising, in ascending order, the 
following formations: Winnipeg (9), Red River 
(10), and Stony Mountain (11). These three for- 
mations are mapped separately only to the 
south of latitude 52° N.; northward the map 
shows Upper Ordovician undivided (12). The 
Stonewall formation (13) forms a well-defined 
band immediately west of the Ordovician rocks. 
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The Middle Silurian age of these dolomites 
might well have been indicated in the legend. 
Devonian rocks (14) are adjacent to the Silurian 
and are suggested to be mainly Middle Devoni- 
an, with possibly some Upper Devonian beds. No 
mention is made of the formational names Elm 
Point, Winnipegosan, and Manitoban which 
have been applied to these strata (Kindle, 1912). 

The age of the Amaranth formation (15), oc- 
curring just west of the Devonian, is not defi- 
nitely known. It is younger than Devonian but 
older than Jurassic. A recent discussion of the 
correlation of this formation is given by Wicken- 
den (1945). 

Strata of undoubted Mesozoic age, mostly 
Cretaceous, are mapped in the southwest por- 
tion of the province. Map unit 16 represents 
Lower Cretaceous and earlier beds, including 
the Swan River group of Lower Cretaceous age. 
The later Cretaceous rocks form the Manitoba 
escarpment and its dip-slope. The Ashville for- 
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mation (17), of Lower (?) and Upper Cretaceous 
age, is succeeded by the Favel formation (18), 
the Vermilion River formation (19), and the 
Riding Mountain formation (20), all of which 
are Upper Cretaceous in age. The Riding Moun. 
tain formation constitutes by far the largest area 
of Mesozoic rock in the map sheet. 

The Boissevain formation (21), of Cretaceous 
or Tertiary age, and the Turtle Mountain for. 
mation (22), of Paleocene age, underlie a small 
area on the International Boundary at longitude 
100° W. The Turtle Mountain formation repre- 
sents the youngest consolidated rock in the 
province. 

The map, which maintains the high standard 
of the Drafting and Reproducing Division, is an 
important addition to the series of more general- 
ized Canadian geological maps. As such it 
should prove to be of considerable interest in the 
classroom and laboratory and will probably cir- 


culate widely. 
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The Ancient Volcanoes of Oregon. By HowE. 
WituiaMs. Eugene: Oregon State System of 
Higher Education, January, 1948. Pp. x+55; 
figs. 9; pls. 13. $1.25. 

This recent work by Professor Williams, of 
the University of California, presents in printed 
form a series of lectures delivered under the 
auspices of the Oregon State Board of Higher 
Education. The lectureship was established in 
1944 in honor of Dr. Thomas Condon (1822- 
1907), the first professor of geology at the Uni- 
versity of Oregon. 

Although written in popular style for the gen- 
eral reader, The Ancient Volcanoes of Oregon con- 
tains a concise summary of Williams’ views on 
volcanism and therefore is of interest to geolo- 
gists. The work also presents an excellent sum- 
mary of the volcanic history of Oregon since the 
Cretaceous period. In the introduction the 
genesis of magma within volcanic regions is dis- 
cussed. The magma below volcanic vents of the 
central type is arrested in its upward ascent 
within a few miles of the earth’s surface and 
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forms a reservoir by “displacing the surround- 
ing rocks.” The magma reservoir then becomes 
the magma-feeding chamber of the volcano. If 
eruptions of a volcano are frequent and the 
magma is continually supplied from the deeper- 
seated basaltic layer of the earth’s crust, the 
products of the eruption will be dominantly 
basaltic. If, however, the magma-feeding cham- 
ber is cut off from replenishment by basaltic 
magma from below, it may and frequently does 
begin to crystallize prior to extrusion, with the 
consequent formation of a gas-rich, highly sili- 
ceous liquid phase, according to the crystalliza- 
tion-differentiation theory of N. L. Bowen. The 
composition of the extruded products will then 
be determined both by the extent to which 
crystallization differentiation has progressed 
and also by the depth or level at which the mag- 
ma-feeding chamber is tapped by fissures ex- 
tending to the surface. Among the minor causes 
of the diversity of lavas, Williams recognizes 
solution or assimilation of the wall rock of the 
magma-feeding chamber. 











eous 
18), 

the 
hich 


area 


eOus 
for- 
mall 
ude 
pre- 
the 


lard 
San 
ral- 








Williams discusses the types of volcanic 
eruptions and correlates these with forms of vol- 
canoes. The introductory section is concluded 
with a brief résumé of the growth of Paricutin 
Volcano, Mexico. 

The main section considers the volcanic his- 
tory of Oregon from the Late Cretaceous to the 
present. In the discussion of the Recent epoch, 
which in Oregon covers approximately the last 
twenty-five thousand years, the sequences of 
events for Oregon’s principal volcanoes are sum- 
marized. These include Mount Mazama (Crater 
Lake), The Three Sisters, Newberry Volcano, 
Diamond Craters, and others. In the discussion 
of The Three Sisters’ region, no mention is made 
of E. T. Hodges’ ancient supervolcano, Mount 
Multnomah (University of Oregon, 1925). 

Although noting eruptions of Mount Shasta 
in 1786, Mount St. Helens between 1831 and 
1854, Cinder Cone near Lassen Peak in 1851, 
and Lassen Peak itself from 1914 to 1917, Wil- 
liams believes that the eruptions of the Recent 
epoch represent a marked decline in the vol- 
canic history of Oregon and that the volcanic 
cycle is drawing to a close. The recent domes of 
rhyolitic obsidian representing the end-product 
of magmatic differentiation are suggested as 
evidence for an old-age decadent stage. Also the 
deeply eroded land forms common in the Cas- 
cade Mountains of Oregon suggest that forces 
of degradation are greatly in excess of those of 
construction. Yet some further outbursts from 
historically active volcanoes of the Cascade 
Mountains and possibly others can be expected 
in the future. 

Williams’ belief in waning volcanic activity 
in the Cascade Mountains may be regarded 
with some skepticism by critical readers. It 
should be remembered that, even though the 
accumulation of volcanic products within the 
twenty-five thousand years since the end of the 
Pleistocene epoch has been relatively minor, 
there were several million years in the Pliocene 
and Pleistocene for the volcanic pile which 
forms the Cascade Mountains to accumulate. 
Possibly more than one volcanic cycle may be 
recorded within this pile, judging from the di- 
versity of volcanic products and the many ero- 
sional unconformities contained therein. If this 
isso, the commencement of a new volcanic cycle 
at any one of several Cascade volcanoes remains 
adistinct possibility, though probably a remote 

one in terms of a human lifetime. 


F. M. Byers, JR. 
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Gesteine und Minerallagerstitten. Vol. 1: Allge- 
meine Lahre von den Gesteinen und Mineral- 
lagerstétten. By P. Niccut. (Unter besonderer 
Mitarbeit fiir Teil TV von E. Niggli.) Basel: 
Verlag Birkhauser, 1948. Pp. 540; figs. 335. 
Sw. fr. 46. 

This book, in itself a store of information 
massed into 540 pages of solid-set type, repre- 
sents only the first volume of a trilogy on rocks 
and mineral deposits. Volume 2, “Die exogenen 
Gesteine und Minerallagerstitten,” is said to be 
in manuscript, and volume 3, “Die endogenen 
Gesteine und Minerallagerstatten,” is in prepa- 
ration. 

The first volume is hardly a suitable textbook 
for students; it is more of a handbook of funda- 
mental facts essential to the understanding of 
the sciences of physical geology, petrology, and 
mineralogy. It is divided into five parts: (I) the 
substance matter (die stoffliche Grundlage): geo- 
chemistry, crystal chemistry, and structural 
mineralogy; (IT) the fabric: structural and tex- 
tural properties of rocks (petrofabrics) and their 
relation to the technical properties; (III) the 
physical chemistry of the formation of mineral 
assemblages: primarily the thermodynamics of 
silicates and silicate melts; (IV) geophysical 
data pertinent to the understanding of rocks 
and mineral deposits; temperature and pressure 
conditions in the crust of the earth and mechani- 
cal, electrical, and magnetic properties; (V) the 
elements of classification and systematology of 
rocks and mineral deposits. Although the works 
of other investigators are frequently mentioned, 
no references to the literature are given. This is 
a definite shortcoming, and there is little con- 
solation in the promise that volumes 2 and 3 will 
contain adequate references. 

Never before has a book of this scope been 
attempted, and no geologist other than Niggli 
could have completed successfully a task so 
diverse and manifold. His superior insight and 
intimate knowledge not only of the geological 
sciences but also of chemistry, physics, and 
mathematics make the book an excellent source 
of information. 

Although lack of space makes it impossible 
for the reviewer to discuss many details, it 
should be mentioned—for the benefit of those 
who know, or think they know, Niggli’s stand 
on the question of granitization—that Niggli 
now willingly accepts this process. He says, for 
example (reviewer’s translation): 

In addition to the metamorphism without mate- 
rial transport we have the metamorphism with more 
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pronounced exchange of matter. The magmatic 
melts in the depth are able to emit matter which 
locally or regionally will penetrate the rocks of the 
lithosphere. . . . Thus truly mixed rocks will form... 
[p. 501]. Rocks which look very much like granite 
may likewise form through ultra-metamorphism, or 
even through ordinary contact metamorphism (e.g., 
from an arkose). If the differences are negligible, or 
if no change of name is justified, then one may speak 
of metagranites (granitic rocks formed through ordi- 
nary metamorphism), or of ultrametagranites (espe- 
cially migmagranites or diabrochogranites). For bet- 
ter definition, and in order to avoid confusion, it may 
be expedient, when various bodies are present in the 
same region, not to use the simple word “granite” 
but, more precisely, to refer to the magmatic rock as 
magmagranite [p. 560]. 


In so far as the reviewer can judge, this 
unique book gives highly needed, correct, exact, 
and thoroughly reliable information, and it is, 
therefore, highly recommended to every physi- 
cal geologist. However, if the American reader 
wants to profit fully, he should have some ad- 
vance knowledge of the scientific jargon used 
by the Niggli school. There is a true jugglery 
with the symbols c, fm, t, alk’, alk, etc.; and a 
simple thing like the chemical composition of 
the feldspars is explained as follows (p. 72): 

Es sind Alumosilikate mit aihnlich gitterhaftem 
Tetraedergeriist wie die SiO.-Modifikationen, jedoch 
teilweisem Ersatz des Si durch Al. Allgemeine For- 
mel [B!¥ O4/2)nAm. In B! ist das Verhiltnis Si: Al = 
3:1, wenn A = K oder Na ist und 1:1, wenn A als 
Ca (eventuel Ba, Sr) auftritt. Auf 1 Al sind zur 
Kompensation der negativen Ladung der gitterhaft 
en Alumokieselsiure 1 Na, bzw. auf 2 Al 1 Ca not- 
wendig [!]. 

T. F. W. B. 


Geological Explorations in Western Tibet. By E. 
Norin. Octavo. Pp. 214; figs. 34; pls. 22; 
colored geological map in pocket. 

There are few regions of the earth compara- 
ble to the vast, desolate, desert highland of the 
Tibetan Plateau, still geologically a terra in- 
cognita. 

Professor Norin, together with Dr. N. Am- 
bolt, started a geodetic and geological recon- 
naissance of western Tibet and the adjoining 
parts of K’un-lun in 1932 and was able to carry 
out an extensive scientific program in spite of 
great hardships and tribulations in this most in- 
hospitable region, which is devoid of almost 
everything necessary for the maintenance of 
both man and beast and is fringed by hostile 
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tribes. The members of the expedition wer 
twice captured and kept imprisoned, once by 
Chinese soldiers, once by armed Tibetans. 

The carefully prepared report is notable for 
its lucidity and its accurate and quantitative 
statements. Excellent descriptions of the geol- 
ogy of the areas are given; but the author goes 
far beyond a mere description; discussions, de- 
ductions, and comparisons with other areas are 
generously incorporated. 

The last chapter is a treatise on the magma- 
tism in Karakoram and western Chang-Thang, 
based on extensive laboratory studies of a great 
number of rock samples: ten new chemical rock 
analyses are given. Furthermore, geometrical 
analyses, optical data, and extended petro 
graphical calculations based on new and old 
data are used for the construction of variation 
diagrams and other graphs illustrating petro- 
genetic relations. 

Also included is a geological map in colors of 
western Chang-Thang and adjacent parts of 
K’un-lun and Karakoram, on a scale of 1:750,- 
000, and covering the area between 34°-38° N. 
and 76°-82° E. 

The book is written in English but has the 
following appendixes printed in German: (A) 
“‘Jungpaléozoische Korallen von Sven Hedins 
Zentralasien-Expedition 1932,” von Franz 
Heritsch; (B) ‘“‘Permische Bryozoen aus Nori- 
tibet,” von Dr. Karl Metz, Leoben; and (C) 
““Mesozoische Fossilien aus dem Tschiptschak- 
tal und Lingschi-Thang,” von Hans Frebold. 


1.2. eae 


Festschrift Paul Niggli, zu seinen 60. Geburtstag. 
(Herausgegeben von seinen Schiilern, Mitt- 
arbeitern und Fachkollegen an den schweizer- 
sichen Hochschulen.) Zurich: Gebr. Leemann 
& Co., 1948. Pp. 568, with photo of Professor 
Niggli. 

This volume contains forty-seven scientific 
papers (which form Heft 1, Band 28, of the 
“Schweizerische mineralogische und _petro- 
graphische Mitteilungen”) and a complete bib- 
liography of the publications written by Pro- 
fessor Niggli (from 1917 through 1948), as well 
as the Ph.D. dissertations completed under his 
supervision. 

It is impossible within the allotted space to 
review individually any of the forty-seven pa- 
pers. Because all papers are written by Niggli’s 
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pupils and co-workers, the complete volume re- 
fects the extraordinary breadth of knowledge 
spanned by the intellect of Professor Niggli: 
mathematical and physical crystallography, oil 
geology, geochemistry, petrography, petro- 
genesis, paleontology, ore geology, crystal chem- 
istry, sedimentology, petrofabrics, soil mechan- 
ies, geophysics, and much more. Indeed, this 
volume is a shining manifestation of the high sci- 
entific standard, vigor, and dynamic activity of 
the Niggli school in Swiss geology. 


Ze Be We 


Technische Gesteinskunde (“Technologic Aspects 
of Rocks’). By A. von Moos and F. DE 
Quervaln. (“Lehrbiicher und Monographien 
aus dem Gebiete der exakten Wissenschaf- 
ten,’ No. 15, ‘““Mineralogisch-geotechnische 
Reihe,”’ Band 1.) Basel: Verlag Birkhauser, 
1948. Pp. viit221; figs. 115. Fr. 29 (un- 
bound); 33 (bound). 


This book, written by two staff members of 
the Zurich Institute of Technology, is a useful 
outline for engineering geologists, especially 
those who deal with problems of tunnellng oper- 
ations, foundations, highway construction, and 
testing laboratories. 

Throughout, those properties of minerals and 
rocks are stressed that need special attention by 
engineers. Descriptions are brief but up to date 
and complete; many determinative and testing 
methods are illustrated; facts and figures are 
compiled in numerous tables; and references to 
technical papers are given in footnotes. Com- 
prehensive works in the field of petrography, 
mineralogy, engineering geology, ceramics, and 
soil mechanics are cited in three pages at the 
end of the book and are followed by an index. 


R. B. 


The “Stratigraph.”’ Los Angeles 27, Calif.: Mo- 
have Sales, Inc., 1948. $7.50. 

The “Stratigraph” is a sheet of flexible plas- 
tic, 10X74 inches, that enables geologists to de- 
termine rapidly and accurately certain unknown 
structural data of inclined strata or tabular rock 
bodies, from those determined at the surface or 
through drill records. 

It is a square, composed of 100X100 small 
squares, and each 10X10 square is marked by 
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a heavier line. Circles concentric with the upper 
left corner connect the too intersection points 
on the upper margin with those on the left mar- 
gin. Around the upper left corner a plastic strip 
with ruled black line rotates across the sheet, 
and tick marks on the right and lower border of 
the square indicate in degrees the angle between 
the strip and the upper margin of the square. 
Thus, for any position of the strip, two series of 
right triangles exist on the sheet, the length of 
whose sides can be read off. 

On the back, diagrams explain the use of the 
grid for twelve selected problems, such as thick- 
ness of beds from dip and depth; depth from sur- 
face distance and thickness; dip of beds from 
depth and thickness. The accuracy of this in- 
strument approaches closely that of trigono- 
metric solutions and, like the latter, is more than 
adequate for most geological examples. For 
speedy answers this is a convenient tool. It is 
sold in a leather envelope. 


R. B. 


Heat Conduction, with Engineering and Geologi- 
cal Applications. By L. R. INGERSOLL, O. J. 
ZoBEL, and A. C. INGERSOLL. (“Internation- 
al Series in Pure and Applied Physics.’’) 
New York: McGraw-Hill Book Co., Inc., 
1948. Pp. xii+278; figs. 42. 


This volume is modeled after the well-known 
text, Mathematical Theory of Heat Conduction 
(Ginn & Co., 1913), by the senior authors, but 
the material has been thoroughly revised; much 
new material has been added, including many 
useful tables and an extensive bibliography; 
the nomenclature has been modernized; but the 
mathematical treatment is essentially the same 
as in the earlier work. 

The development and explanation of heat- 
conduction equations of many types is, for the 
most part, detailed enough to be followed with 


-a minimum of difficulty by readers moderately 


familiar with calculus; ellipsis is evident chiefly 
in the occasional differentiation of involved 
functions—the authors seem to make the er- 
roneous but flattering assumption that, al- 
though transformations and developments of 
equations involving integrals may be somewhat 
difficult, all readers have the art of differentia- 
tion at their fingertips. However, the use of the 
important basic equations is so well illustrated 
by fully worked-out examples that a geologist or 
engineer should be able to use the equations 
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whether the derivation is clear or obscure; but 
one would miss much of the fun of mastering a 
problem if unable to follow the explanation of 
the general case. The writers are to be congratu- 
lated for writing a book on a highly mathemati- 
cal subject that is so generally readable and that 
should be understandable by anyone with the 
equivalent of a one-year course in calculus. This 
reader, however, found the frequent use of 
“while” (properly a time term) for an adversa- 
tive conjunction, where “although” would bet- 
ter carry the meaning, to be singularly unfelici- 
tous in a book primarily concerned with time 
relations. 

Chapter i (“Introduction”) has been much 
improved; the tables of nomenclature and con- 
version factors and the expanded discussion of 
the application of conduction-equation mathe- 
matics to fields other than heat conduction are 
helpful improvements; the use of dimensional 
formulas could well have been preceded by a 
brief explanation and followed by a reference to 
their use in modern theory. 

Chapter ii (““The Fourier Conduction Equa- 
tion’’) contains the same material found in chap- 
ter ii of the earlier book. It presents a clear but 
concise derivation of the Fourier equation, to- 
gether with a brief section on the use of bound- 
ary conditions in determining particular solu- 
tions of the general equations. 

Chapter iii (“Steady State, One Dimension’’) 
illustrates the plan used in most of the remain- 
ing chapters of the book. The problem to be 
studied is first explained, the general form of the 
Fourier differential equation that applies is pre- 
sented, and the boundary conditions that must 
be satisfied are enumerated. This is followed by 
a brief, but adequately explained, derivation of 
an algebraic or integral equation that will allow 
arithmetical substitutions, satisfy the boundary 
conditions, and be a particular solution of the 
Fourier equation. After treating the general case 
of one-dimensional heat flow, a section is devot- 
ed to illustrative applications of the equations 
derived in the preceding section, and the chap- 
ter closes with a half-dozen problems, with the 
answers given, for the student to work out and 
check. The examples chosen in this chapter are 
of interest to engineers rather than to geologists 
and deal chiefly with problems of insulation 
and contact resistance. 

Chapter iv (“Steady State, More Than One 
Dimension”’) is also concerned with insulation 
and cooling problems of engineering interest, 
but section 4.10, “Subterranean Temperature 
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Sinks and Power Developments; Geysers,” cop. 
tains some material for the geologist. 

Chapter v (‘Periodic Flow of Heat in One 
Dimension”) treats several geologic problems; 
heat waves—cold waves, periodic heat flow and 
climates, and thermal stresses all suggest possi. 
ble applications to our field. 

Chapter vi (‘‘Fourier Series”) is identical 
with chapter vi of the earlier volume; it is de- 
voted to the step-by-step derivation of the 
Fourier sine and cosine series that must be used 
in many of the particular solutions of the Fov- 
rier heat-conduction equation, where unsteady. 
state heat-flow problems are to be solved. It is 
a pleasure to see such a generally readable ex- 
planation of a somewhat difficult subject, but it 
is regrettable that the occasional ellipsis present 
in the earlier work has not been eliminated. 

Chapters vii and viii (“Linear Flow of Heat, 
I and II’) take up many problems of heat flow 
in one dimension (taken as parallel to the «x-di- 
rection). Subjects of interest to geologists will 
be found under several headings: the ‘‘concrete 
wall,” mathematically identical with the prob- 
lem of a cooling dike, also covered by problem 4 
(p. 85) (intrusion of “a slab of molten lava”); 
the heating or cooling of a semi-infinite body 
whose surface is held at constant temperature; 
temperature of contact of thermally dissimilar 
bodies; the thawing of frozen soil; the cooling ofa 
lava sheet; the cooling of the earth, assuming an 
initially hot homogeneous mass—both with and 
without radioactive heat sources; the tempera- 
ture in rock decomposing exothermically; post- 
glacial climate calculations from temperature 
gradient measurements; and the heating or 
cooling of slabs suddenly immersed in a medium 
of different temperature (theory of the fireproof 
wall). 

The mathematical discussion of the cooling 
of the earth, on the assumption that a quarter 
of the heat loss is due to radioactivity, gives a 
figure of 45,000,000 years for its age, as com- 
pared with 22,000,000 years for a nonradioac- 
tive earth; the concluding paragraph of this sec- 
tion unfortunately states that these estimates 
are “of the same order of magnitude as those ar- 
rived at from geological considerations, such as 
stratigraphy, sodium denudation, etc.,... al- 
though far short of the 2 X 109 years which rep- 
resents the present trend of thought.” As no 
explanation of the discrepancy is offered, some 
engineering readers will be misled, and geologists 
may feel baffled and disturbed. In the section on 
temperature in decomposing granite, references 
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should have been cited for the statements (p. 
115) that “areas of granite undergoing decom- 
position are several degrees warmer than sur- 
rounding rock” and that granite gives out heat 
during decomposition of the order of 100 cal/gm. 

Chapter ix (“Flow of Heat in More Than 
One Dimension, Unsteady State”) includes 
such geologic applications as the cooling of a 
laccolith, heat sources and heat sinks as related 
to geysers, the general cases for the cooling of a 
rectangular parallelopiped and of a finite cylin- 
der in material of the same sort, and the use of 
conduction equations in such liquid-diffusion 
problems as the drying of porous solids. A help- 
ful addition to this part of the book is the table 
on pages 184-185, which summarizes the heat- 
conduction equations for many differently 
shaped bodies under specified conditions. 

Chapter x (‘Formation of Ice”’) is essentially 
the last chapter of the earlier text and gives 
both Neumann’s and Stefan’s solution of the 
problem of ice formation for one-dimensional 
heat flow involving latent heat of fusion in a 
semi-infinite body with its surface at constant 
temperature. 

Chapter xi (‘Auxiliary Methods of Treating 
Heat Conduction Problems’’) gives several ap- 
proximation methods of solving heat-conduc- 


tion problems, a few of which may take less 
time than the analytical methods which they re- 
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place. Some, however, can be used to advantage 
where the problem is too complicated to be han- 
dled readily by classical mathematics. The ex- 
planation of the step method, as illustrated in 
the problems of ice formation about pipes, and 
the warming of soils, will be found helpful by 
geologists in a variety of problems. The brief 
section on electrical methods of solving heat 
problems is good, but it is to be regretted that 
the use of thermal models is entirely neglected. 

Chapter xii (“Methods of Measuring Ther- 
mal Conductivity Constants”) concludes the 
body of the text with a brief description of sev- 
eral methods of determining &. It is followed by 
thirteen appendixes: thermal constants, an ex- 
tensive bibliography, and many mathematical 
tables especially useful in the solution of heat- 
conduction equations are given. A good index 
completes this excellent work. 

The book is well illustrated with diagrams 
and graphs, the print is clear and easily read, 
the paper is good, and the text seems singularly 
free from typographical errors—only two minor 
ones were found by this reviewer in a cover-to- 
cover perusal of the volume. The wealth of ma- 
terial of diverse geologic interest and the lucid 
presentation of analytical methods should make 
it a most valuable and welcome addition to any 
geologist’s library. 

T. S. LOVERING 





COMMUNICATIONS AND ANNOUNCEMENTS 


CAMBRIAN STRATA ON THE SOUTHERN SHORE OF DEVON 
ISLAND IN THE ARCTIC ARCHIPELAGO: 


The southeastern part of Devon Island con- 
sists of pre-Cambrian rocks unconformably 
overlain by almost horizontal strata heretofore 
designated Silurian. 

A geological reconnaissance of the Dundas 
Harbour area during the summer of 1948 by 
R. P. Nickelsen and D. B. Wales revealed that 
at least 700 feet of fossiliferous Cambrian lime- 
stones are present on the western side of Dundas 
Harbour. 

The harbor extends inland 3.5 miles and is 1 
mile wide at its entrance. A cliff approximately 
2,400 feet high rises from sea level along the 
west side of the harbor. The lower 1,700 feet of 
the cliff consist of pre-Cambrian gneiss intruded 
by several diabase dikes. The upper 700 feet 
consist of dolomitic, arenaceous, and argilla- 
ceous limestones. The surfaces of the limestones 
are intensely fractured. Limestone blocks form 
talus in many places and cover the upper part 
of the pre-Cambrian rocks. Brownish dolomitic 
limestones, from within 20 feet of the uncon- 
formity, contain scattered brachiopod shells 
identified by Dr. W. C. Bell as Micromitra. The 
presence of Micromitra suggests strongly that 
the basal beds of the limestone series are Middle 
Cambrian. Approximately 400 feet above the 
unconformity, a thin-bedded limestone talus 
block yielded numerous trilobite glabella of the 
genus Elrathia. These also indicate a probably 
Middle Cambrian age. 

The opportunity to study the rocks of the 
Dundas Harbour area was made possible by the 
United States Weather Bureau. Special ac- 
knowledgments are due to Professor R. L. 
Nichols, of Tufts College, and Mr. Charles Hub- 
bard, chief of the Arctic Section, United States 
Weather Bureau, for aid in the field work, and 
to Mr. Vincent E. Kurtz, of Dartmouth Col- 
ege, who aided in the identification of the fos- 
sils. 

Tentative plans for the summer of 1949 in- 
clude a detailed study of the stratigraphy of the 


* Manuscript received February 7, 1949. 


Dundas Harbour area to obtain information ¢ 
the thickness and age relationships of the lowe 
Paleozoic rocks of this interesting and practi 
cally unstudied area. 4 

DonaLp B. WA 
DARTMOUTH COLLEGE 


REPRINTING OF KRYNINE ARTICLE 


‘The Megascopic Study and Field Classifie 
tion of Sedimentary Rocks,” by Paul D. K 
nine (originally published in the Journal of 
ology, vol. 56, March, 1948, pp. 130-165; 
tables; 15 figures), the stock of which has bee 
exhausted, has now been reprinted as Technical 
Paper 130, and is available at 25 cents per copy. 
The limit will be 50 copies per institution or or? 
ganization. A 

Orders should be sent to: P 


MINERAL INDUSTRIES EXPERIMENT STATION 
THE PENNSYLVANIA STATE COLLEGE 
STATE COLLEGE, PENNSYLVANIA 


PENNSYLVANIA GEOLOGISTS’ FIELD © 
CONFERENCE : 


The Field Conference of Pennsylvania Geolo- 
gists will hold its annual meeting May 27-29 at 
Lancaster, Pennsylvania. The department of 7 
geology of Franklin and Marshall College will” 
be host. The major feature of this year’s field” 
excursion will be a trip through southern Lan- 
caster County to study the Martic overthrust = 
problems under the leadership of Professor 
Ernst Cloos, of the Johns Hopkins University. | 
Other field excursions will include the old Wood 
chrome mine, Gap nickel mine, other old metal 
mines, and the flood plain of the Susquehanna — 
River. q 

For further details, address Dr. Richard M. 
Foose, Department of Geology, Franklin and 
Marshall College, Lancaster, Pennsylvania. 
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